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ABSTRACT 


When using the Instrument Landing System (ILS) all aircraft must 
follow a single straight line approach path before landing. The 
Microwave Landing System (MLS) will allow use of differing approach 
paths and they need not be on a straight line. 

The objective of this research is to find out whether the intro- 
duction of MLS with its multiple approach path capability can bring 
an increase in runway landing capacity compared with conventional ILS. 

A model is developed which is capable of computing the ultimate 
landing runway capacity, under ILS and MLS conditions, when aircraft 
population characteristics and Air Traffic Control separation rules 
are given. This model can be applied in situations when only a hori- 
zontal separation between aircraft approaching a runway is allowed, 
as well as when both vertical and horizontal separations are possible. 

It is assumed that the system is free of errors, that is that aircraft 
arrive at specified points along the prescribed flight path precisely 
when the controllers intend for them to arrive at these points. 

Although in the real world there is no such thing as an error-free 
system, the assumption is adequate for a qualitative comparison of MLS 
with ILS. 

Results suggest that an increase in runway landing capacity, 
caused by introducing the MLS multiple approach paths, is to be expected 
only when an aircraft population consists of aircraft with significantly 
differing approach speeds and particularly in situations when vertical 
separation can be applied. Vertical separation can only be applied if 
one of the types of aircraft in the mix has a very steep descent angle 
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(i.e. 7.5 degrees) such as an STOL vehicle. When approaching aircraft 
are separated only horizontally, examples considered in this research 
show a modest capacity increase of 10 to 15 percent. When both verti- 
cal and horizontal separations are applied, capacity improvement can 
be greater depending on the proportion of steep descent aircraft in 
the mix. 

It was also found that the angles of entry to the extended runway 
centerline have a significant effect on landing capacity, and they 
should be optimized. 
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1. Introduction 

Improved navigational aids for approach and landing at airports 
have been under development for many years. One such aid is the Micro- 
wave Landing System (MLS) which provides multiple flight paths to the 
runway rather than a single path defined by the current Instrument 
Landing System (ILS) . 

There are a number of advantages cited for the MLS . Among these 
are (1) curved path approaches can reduce flight over noise sensitive 
areas, (2) the system is less sensitive to interference from terrain 
and man made objects, (3) since the system extends much farther from 
the runway than the current ILS, aircraft have precise guidance over 
their intended flight paths for a longer period of time before landing, 
and (4) flexibility in flight paths might increase the landing capacity 
of a runway. 

This research deals with item 4, landing capacity. Models for 
capacity are developed to reflect the multipath capability of MLS. The 
models are then applied to hypothetical situations, and the capacities 
obtainable with MLS and ILS are compared to determine if there are any 
significant differences . 

In this research it is assumed that both the ILS and MLS are free 
of any errors ; that is , that aircraft arrive at points in space when 
the controllers intend them to be there (e.g., at the entry gate to 
ILS). Therefore, no buffer is added to interarrival times. It is 
recognized, however, that in the real world there is no such thing as 
an error free system. 

The objective of this research is, then, to compare in a qualita- 
tive manner the ILS and MLS systems to see if there is a significant 
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difference in their respective capacities, and if so under what condi- 
tions . 

The term capacity as used in this research refers to the maximum 
number of landings that a runway can accept when there is a continuous 
demand for service and a certain specified set of conditions (i.e., air- 
craft mix, air traffic control (ATC) rules, etc.). These conditions 
can significantly effect runway landing capacity. 

This research focuses on one of the factors that influence capacity: 
flight path geometry, particularly that geometry which ensures maximum 
landing capacity for a specified set of the other factors that influence 
capacity, i.e., aircraft population and mix, length of common approach 
path, ATC rules. 

1.1. Characteristics of ILS 

The existing ILS is essentially a straight line in three-dimensional 
space (see Fig. 1 ) . This line ends on the runway. Aircraft follow the 
line and land on the runway. There are three pieces of information that 
the pilot obtains about the position of his aircraft, with respect to 
the line leading to the runway and the distance from threshold: 

1. position of the aircraft with respect to the alignment of the 
runway; namely, whether the aircraft is left or right of the 
centerline of the runway; 

2. position of the aircraft with respect to the required height 
above the runway, referred to as the glide path; namely, 
whether the aircraft is above or below the glide path; 
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All figures are placed in Chapter 7 . 
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3. distance of the aircraft from the runway threshold. Distance 
information is provided by markers (two or three) : one at 

the beginning of the ILS approach and the other (or two others) 
nearer to the runway threshold. 

With the current ILS (see Fig. 1), landing aircraft follow each other 
on a common path, ET (E being the "entry gate" and T the runway thresh- 
old). Normally there are two markers installed on the common approach 
path, except for ILS Category II (or lower) weather conditions when a 
third marker is added. The marker furthest from the runway (about 5 nm) 
is the "outer marker" (OM) ; the marker generally closest to the runway 
(about 0.6 nm) is the "middle marker" (MM); and the sometimes used 
third marker (about 0.2 nm) is known as the "inner marker" (IM) . 

1.2 . Characteristics of MLS 

The MLS provides glide path information up to about 15 degrees 

elevation and alignment information as much as 130 degrees relative to 

.... / 

the runway (65 degrees each side of the center line of the runway) . In 

J 

lieu of markers, continuous information on. distance to the runway will 
be incorporated into the MLS. The MLS is shown in Fig. 2. 

MLS gives information on the position of the aircraft relative to 
the runway in three-dimensional space. The area of coverage can be 
described as a quasi-pyramid with the runway threshold as the apex (see 
Fig. 3). The three-dimensional information is continuous and provides 
a means of describing different paths for aircraft to follow in the 
space covered by MLS. However, even if the information on the position 
of aircraft is accurate and the equipment is available to program any 
type of trajectory inside the pyramid, all described paths cannot 
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be followed by an aircraft. Several assumptions need then to be made 
about possible restrictions on the types of trajectories, as follows. 

1. It is assumed that aircraft need to fly along the prolongation 

of the centerline of the runway before landing (E „ „ to 

x— i ,Z , J 

T , Fig. 3). is the "entry gate" for this straight por- 

tion of the final approach for aircraft of type i . 

2. It is assumed that there are some restrictions to the curved 
paths because of the minimum turning radius of an aircraft. 

3. It is assumed that there are restrictions to the maximum angle 
of descent. 

4. It i > assumed that there are sufficient exit taxiways on the 
runway so that runway occupancy time is always less than the 
threshold interarrival time which ensures the ATC-required 
separation of aircraft in the air. 

These restrictions, together with air traffic control (ATC) separa- 
tion rules , limit the number of usefully considered approach paths . 

As section 2.2 indicates, within these limits rather simple flight paths 
are chosen, these representing the most desirable paths from the stand- 
point of runway capacity. 

1.3. Existing ATC Separation Rules 

According to ATC rules , aircraft can be separated vertically and 
horizontally . Horizontal separation can be expressed in time or dis- 
tance. It is here assumed that radar coverage is available, and there- 
fore separations are in terms of distance rather than time. Two air- 
craft approaching the runway should, then, never to closer to each 
other than the minimum horizontal distance prescribed by ATC. 
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Another important safety regulation, one that {“tyuld influence the 
capacity of a runway, is that two aircraft can not be on the runway at 
the same time; the first aircraft has to clear the runway before the 
second crosses the threshold. (It is assumed that runway occupancy 
time is always less than that threshold interarrival time, which ensures 
the ATC required separation of airborne aircraft: therefore, runway 

occupancy time is not a constraint.) 

1.4. Runway Landing Capacity Model Basic Structure 

Basic landing capacity models using the current ILS were developed 

1 

by Dr. Richard Harris of Mitre Corporation and refined and expanded by 

2 

Peat, Marwick, Mitchell and Co. The landing capacity of a runway is 
defined as that maximum number of landing operations that can take place 
on a runway in a unit of time (usually one hour) during which aircraft 
continually wish to land. This concept of capacity is often referred 
to as "ultimate" or "saturation" capacity. The maximum number of land- 
ing operations on the runway depends on a number of conditions, as 
follows : . 

1. minimum separation rules specified by ATC; 

2. aircraft mix (i.e., the proportion of different types of air- 
craft using a runway in a given period of time) ; 

3. location and type of exit taxiways (i.e. , if there are an in- 
sufficient number of taxiways , the runway occupancy time 
rather than air separation might be critical) ; 

4. geometry of the approach paths associated with the runway 
(i.e., multiple paths are available with MLS). 

To compute capacity, each of these conditions must be specified. 
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Runway land?«^ capacity is actually the capacity of a system, con- 
sisting of a runway and the airspace adjacent to the runway. The run- 
way is that part of the system where the flow of all aircraft opera- 
tions converge; consequently, the landing capacity can also be defined 
as that number of operations in unit time which pass through a point 
which all aircraft have to pass. In this analysis this point is taken 
as the runway threshold. To find the capacity it is then necessary to 

determine t. . , it being defined as follows: t.. = in ter arrival time 

1 3 

at the threshold between aircraft type i (lead aircraft) and aircraft 
type j (trailing aircraft); t_^_. should then be such that: 

a. aircraft i and aircraft j will not occupy the runway at 
the same time (i.e., when j passes the threshold i should 
have cleared the runway) ; 

b. in the air, aircraft i and aircraft j are never closer 
than the minimum separation specified by ATC rules, then. 


t. . 


min ( t . . , t . . ) where 
a xj r xj 


t. . = interarrival time at the threshold, dictated by ATC 

a x;j 

minimum separation rules for airborne aircraft, 


t. . = interarrival time at the threshold, dictated by ATC 

r xj 

runway occupancy rule: only one aircraft can occupy 

the runway during any interval time . 

Existing models of capacity assume independence, that is, the type of 
trailing aircraft j is not dependent on the type of leading aircraft 
1 . The proportion of aircraft in the mix over the unit of time being 
considered (usually one hour) has, however, to be preserved. This 
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implies the following: 





p.p. , where 
i J 

probability of the sequence ij , 
proportions of aircraft i and j in the mix. 


When t. . and p. . are found for all i and j , the expected 
ij 

interarrival time at the threshold can be computed as': 


t = 2 t p 

ij ^ 


The capacity, assuming independence in the sequence of aircraft but 
with the restriction that the aircraft mix will remain constant during 
the unit of time selected ( p^ = 1» where p fc = proportion of aircraft 
of type k ) , is 



where X - landing capacity 

f “ expected interarrival time between aircraft over the 
runway threshold. 


1.5. Factors Affecting Runway Landing Capacity 

If it is desired to maximize X , it is necessary to minimize t 

(i.e . , min 2 t,.p. .). If sequencing is not applied to the stream of 
ij ij 

the landing aircraft, with given p, , all p . . are fixed. Conse- 


quently, the term that requires analysis is t^ . . 

As noted, it is assumed that air separation rather than the runway 
occupancy rule is critical in all landing cases , so that t^_. = • 

Two more assumptions are made in the following analysis. 

1. All aircraft have uniform velocities when they are on the final 
approach to the runway. 

2 . All aircraft of the same type use. the same path for approach to 
the runway. 

These two assumptions are critical to an analysis of capacity. Consider 
simple case, only two aircraft types, "fast* 1 and "slow." We see immedi- 
ately that there are four interarrival times: 

t„ n = interarrival time between "fast" followed by "slow" 

"fast” 

"fast" 

"slow" 

Cases t„ and t co are simpler than the other two, because the land- 
ing trajectory is the same for both aircraft and the distance between 
the two aircraft measured along the trajectory is constant. The dis- 
tance has, however, to be such that during the entire approach the two 
aircraft never come closer to each other, measured on a straight line 
in the horizontal plane, than the minimum specified by ATC rules . 

The situation for t and t is more complex. First, the 

Fb br 

two aircraft might not have a common trajectory, in which case the 
separation can not be measured along the trajectory but rather by 
horizontal, vertical or diagonal separation. Second, if part of the 
trajectory is common to both aircraft (see aircraft 1 and 3 in Fig. 3) , 




“SF 


'FF 


'SS 


"slow" 

"fast" 

"slow" 




the distance between them measured on the common path is not constant; 

it depends on the relative speeds of the two aircraft. The distance 

and corresponding time will increase for v, < v. (aircraft type i 

J 

followed by aircraft j) and decrease for v. > v. . 

j i 

As noted, there are capacity models available for computation of 

ultimate capacity. These models are, however, primarily applicable 

to a single trajectory corresponding to current ILS procedures. The 

structure of these models was briefly described in paragraph 1.4. The 

1 2 

method for computing t„ in these models is shown in Fig. 4. ’ This 
figure represents a time-space diagram in the horizontal plane. Time 
is the abscissa and distance is the ordinate. 

Several remarks concerning the procedure for computation of t^ 
are in order to develop a better understanding of the assumptions. 

Figure 4 shows that all aircraft are fed into the entry gate (E) 
from the same path, an extension of the current single ILS alignment 
path which is, in turn, an extension of the center of the runway. An 
important point is that the state of the system is considered only when 
t > 0 (i.e. , only for that time after the instant, t = 0 , when the 

leading aircraft passes through the entry gate) . 

In the case of sequence FS (fast followed by slow) , shown in 
Fig. 4a, it can be seen that the horizontal distance between the two 
aircraft, d(t) , is equal to the ATC minimum required horizontal sepa- 
ration, 6 , when t - 0 and increases thereafter. It can also be 
seen that the horizontal separation constraint, d(t) > 6 , is violated 
prior to t = 0 . (See dotted lines in Fig. 4.) 

There are three possible explanations how one could have a slow 
aircraft only 6 behind a fast aircraft when the last aircraft is at 
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the entry gate , E . 

The first is to have the approach path unchanged in the horizontal 
projection and assume that the velocities of an individual aircraft are 
not uniform during the approach. Velocities of the two aircraft should 
be at least equal (v £ v ) for t < 0 . As it is not likely that v 

b r r 

will be smaller when t < 0 than when t > 0 , v g has to be increased 

until it at least equals v„ (see Fig. 5) . This increase in speed is 

r 

not very practical, especially in the case where two approach speeds 
are very different. 

The sequence FS is of particular interest since the interarrival 
time, t c , is critical from a capacity point of view since the time 

rb 

gap between F and S opens as they approach the runway along the 
common path (see Fig. 4a). In addition, if it is assumed that an in- 
crease in speed of the slow aircraft (vg) is possible, the problem of 
the two aircraft being separated exactly 6 when t = 0 remains only 
partially solved; i.e., it is still assumed that both aircraft are on 
the straight line (a prolongation of the runway centerline) when t < 0 
and are coming from infinity on this line, continually separated by at 
least S . This assumption does not of course fully consider the real 
world since one of the following events must have occurred before 
t « 0 . 

1. Aircraft F has overtaken aircraft S (observed in the hori- 
zontal projection) on the straight line, indicating that ver- 
tical separation was imposed (the fast aircraft went either 
under or above the slow) . 

2. Aircraft S joins the straight line path (observed in the 
horizontal projection) after aircraft F has passed through 
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the point E , this point being where the slow aircraft joins 

u 

the common path. 

A second possibility is that two aircraft can be separated verti- 
cally in the approach air space before t = 0 . Fig. 6 shows how this 
separation might look using distances which are appropriate to current 
ATC rules and procedures. It is assumed that when t = 0 , vertical 
separation is imposed in such a way that the slow aircraft is above the 
fast and both are flying level. Two conditions are necessary to ensure 
the vertical separation shown in Fig. 6. 

1 • h - x 

where h is vertical distance between the two aircraft and x 
is the minimum permissible vertical separation. At present 
ATC rules specify x as 1»000 ft (305 m) . 

2. The fast aircraft should be able to perform straight level 

approach at an altitude equal to threshold elevation plus H , 

a distance of 2 „d (see Fig. 5), before intercepting the 
I c 

glide path at the entry gate , E . 

In Fig. 7 another possible vertical separation for aircraft approaching 
the entry gate, E , is shown, when the fast aircraft is at E , and 
the slow one is a distance of S behind and x below the fast. Two 
conditions are necessary to ensure this vertical separation. 

1. (H - X)/sin 0 has to be greater than the distance necessary 
for the slow aircraft to stabilize on the glide path. If this 
condition is not satisfied, the resulting problem could theo- 
retically be solved by the slow aircraft climbing to intercept 
the glide path after t - 0 . The slow aircraft climbs — or 
climbs and levels — until it intercepts the glide path (see 
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2. The second condition is similar to the second condition of 

the previous vertical separation case (a slow aircraft above 

a fast); but it is far more difficult to satisfy: the slow 

aircraft is required to perform straight level approach at an 

altitude equal to threshold elevation plus H - X * f° r a 

distance of 2 d - 6 before reaching point F . Consider- 
E c 

ing. the values of H (shown in Fig. 6) this condition is al- 
most impossible to satisfy since X is 1>000 ft, i.e., slow 
aircraft would be required to fly level at a height above the 
terrain of only 500 to 1500 ft which is not acceptable (see 
Fig. 7). 

The third possibility of having the slow aircraft only 6 behind 
the fast one, when the fast one is at E , is briefly discussed below. 

Does, for example, a path leading to E exist, such that, if the 
slow aircraft is behind the fast on this path separated by distance 5 
(measured along the path) when the fast passes E (t = 0) , the distance 
between the two aircraft (measured as a straight line in the horizontal 
plane) is never (when t < 0 and when t > 0) less than 5 ? (See 
Fig . 8a.) 

To satisfy the condition d(t) > 6 , d(t) being the distance be- 
tween the slow and fast aircraft measured along a straight horizontal 
line, for all values of t , it is also necessary to satisfy this con- 
dition when t = 0 . Obviously, this condition cm only be satisfied 
if the path of the slow aircraft to E is on a straight line leading 
to E for at least length 8 before E is reached. This implies 
that at the moment when t = 0 the slow aircraft should be somewhere 
on the circle of radius 5 which has its center at E (see Fig. 8b). 
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However, strictly speaking, though the slow aircraft is on this circle 
when t = 0 , the distance that the slow aircraft has to fly along the 
path to E does not equal 5 ; in fact, this distance is more than <5 
with the exception of position 1 (see Fig. 8b) for the following reason 
The distance that the slow aircraft must fly to reach E is greater 
than 6 because of a change in heading which is required for positions 
2, 3, and 4. This change can only be made on a curved path; so the 
sum of the straight and curved path will always be greater than the 
radius of the circle. However, even if the slow aircraft is at point 1 
when t = 0 , it will violate the condition d(t) k 6 when t = -At 
(see Fig. 4). So we can conclude that strictly speaking the path we 
were looking for does not exist. 

In spite of the shortcomings which have been discussed in this 
section, existing runway capacity models can be used to find potential 
sources for an increase in landing runway capacity. Fig. 4 allows us 
to examine interarrival times, t g , t gF , t^ and t gg , and at- 
tempt to decrease them, as a way of decreasing mean interarrival time, 
t , and consequently increasing landing capacity, X . For ease of 
computation, v g can be expressed as v g = yv^, , 0 < y S 1 . 

From Fig. 4 

£ se ■ I; 

t = «_ 

FF v p 

fc ss = 


, and 


• • - 1 1 t - 1 l i 
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t m 1_ f S + Y _ v *| „ 1_ 6 +y(l - u ) = 6 ± Y(1 - yi) 

FS V F 1 y YJ V F y u v F 

Possible ways to decrease t. . are then as follows. 

iJ 

1. The most efficient would be to decrease 6 , thus simultane- 
ously decreasing all t F j * s • 

2. Increase both v e and v , thus ensuring a higher velocity 

b r 

of aircraft stream at the runway threshold which results in a 
higher flow. 

3. Increase ]i only, keeping v constant. (This can be achieved 
by requesting slow aircraft to maintain cruise speed as long as 
possible.) 

4. Decrease y , the common part of final approach. 

The implementation of MLS could reduce 6 (case 1) and y (case 4) . 

In summary, then, this research is concerned with whether the ap- 
proaches along the curved paths described by MLS will improve the land- 
ing capacity of a runway, and, if the answer to this question is posi- 
tive, what is the value of this increase over that capacity obtained 
with ILS? 

The following analysis of capacity will begin with the assumption 
that vertical separation within the operating area of MLS is not per- 
mitted, i.e., horizontal separation is crucial. The situations governed 
by this assumption are examined in Chapters 2 and 3. The assumption 
that vertical separation may be permitted (either horizontal or verti- 
cal separation could therefore be employed) governs the situations 
examined in Chapters 4 and 5. 
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2 . Analysis of Capacity With Horizontal Separation Only 
2.1. Introduction 

In this chapter the influence of multiple flight approach paths on 
the landing capacity of a runway is analyzed. It is assumed that all 
approaching aircraft are on the same horizontal plane, that is, no 
vertical separation exists. 

The procedure for the current ILS is shown in Fig. 9. All aircraft 
landing in IFR have to pass through a point E (entry gate) , located 
distance y (common approach path length) from the landing threshold, 

T , along the extended runway centerline, and remain on the centerline 
until they reach the threshold. As shown in Fig. 10, with MLS there is 
no need for aircraft to use a common approach path of length Y . Sim- 
ilarly, all aircraft do not need to pass through the common entry gate, 

E . Multiple paths are possible. Each path intersects the prolonga- 
tion of the runway centerline at E_^ , and, as noted in Chapter 1, 
Section 1.5, this path is used only by aircraft of type i . Each 
intersection of the runway centerline, E_^ , is located y^ from the 
runway threshold. Each intersection, E_^ , can be considered as the 
entry gate for type i aircraft. 

The objective of the analysis is to develop a model to describe 
multiple rather than common entry path geometry. The model should rep- 
resent the two approach path situations (ILS and MLS) shown in Figs. 9 
and 10 and determine the ultimate capacities of the two (under given 
conditions) . The model should also be able to determine a geometry 
(i.e., set of approach paths) which maximizes ultimate capacity for 
either the ILS or MLS situation. Once the maximum ultimate capacity 
for each of the two situations has been determined the effect of MLS on 
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runway landing capacity can be obtained by comparing the two capacities. 

Analysis of the effect of geometry on the ultimate capacity, either 
with ILS or MLS, will concentrate on an aircraft pair consisting of two 
aircraft types, a fast aircraft, F , and a slow aircraft, S . When 
analyses of cases involving more than two aircraft types in a popula- 
tion are necessary they can be treated in pairs. These pairs in all 
cases consist of fast and slow aircraft. This approach of breaking 
down an aircraft population in couples and dealing with four general 
case interarrival times, t^ , t^g , t gF , and tgg (i.e., t g - 
interarrival time at threshold for a fast aircraft followed by a slow 
one, etc.) will be used throughout the analysis since any of the inter- 
arrival times, t^_. , can be represented by one of these four times. 

From Fig. 4 we can draw some rough conclusions about aircraft 
operations with ILS. In conventional ILS all aircraft pass through an 
entry gate, E , and travel on a common approach path through a dis- 
tance of Y to the runway threshold. As Fig. 4a illustrates, in the 
case of a slow aircraft following a fast one, the interarrival time 
over the threshold, t^g , increases as Y is increased. Similarly, 
in the case of MLS for a situation involving a fast aircraft followed 
by a slow, J , the common path for the two aircraft should be con- 

r o 

sidered rather than the single common path, Y , where 

T FS - min(Y F> Yg) . 

A new assumption is therefore introduced; it will be used in fur- 
ther analysis; it is 




This assumption requires that an aircraft with a higher approach 
speed will need at least as long a final straight approach along the 
extended runway centerline as an aircraft with a lower approach speed. 
That is. 


v = y 

'fs 's 


i.e., the common approach path is equal to the length of the final 
straight approach for slower aircraft. 

Therefore, if Y c < Y » is decreased. However, even if 

b r b 

Y c = 0 (i.e., that a slow aircraft does not need any portion of straight 

b 

level wing approach before landing) the interarrival time, t^ s , re- 
mains equal to 6 /Vg (see Fig. 4a). 

Fig. 4 also illustrates that if Y is reduced the only possible 
reduction in interarrival times that can be made is a reduction of t g . 
The other three times, t , t gF , and t gg , will remain unchanged. 

Intuitively, it can be expected that an increase in landing capac- 
ity will be small if a decrease in t is the only reduction of inter- 
arrival time possible. For example, consider the case where an aircraft 
population is equally divided between fast and slow aircraft, that is, 
p = 0.5 and p = 0.5. If the sequence of the arrivals is considered 
as random, then the fraction of slow aircraft following fast aircraft 
will be 0.25, which is the result of p f x p g . Seventy-five percent 
of the arrivals will constitute pairs of aircraft with interarrival 
times, t_,_ , t , and t . These times cannot be changed. There- 
fore, a significant reduction in the mean interarrival time, t , and 
consequently, an Increase in landing capacity, A , cannot be expected. 


This- conclusion will be examined in detail in the following sections* of 
this chapter. 

As noted in Chapter 1, Section 1.5, in order to decrease t a 

r b 

slow aircraft must be brought as close as possible to the threshold 
when a fast one is over the threshold. Minimum horizontal separation 
rules must not, however, be violated. 

The following examination of the effect of MLS on runway capacity 
will involve analysis of changes of interarrival times over the thresh- 
old, t_^ . The first step in this analysis is to see how these changes 

respond to certain variables. 

The variables that affect t. . are as follows: 
v.p = velocity of the fast aircraft 
v Q = velocity of the slow aircraft 

b 

11 = V V F 

Y_, = length of the approach on extended runway centerline 

r 

of the fast aircraft 

Yg = length of the approach on extended runway centerline 
of the slow aircraft 

5 = minimum horizontal separation required between the two 

aircraft measured in the horizontal plane. 

If the values for all of these variables are given, then the effect 
on capacity of various possible aircraft paths , herein referred to as 
trajectories, is necessarily the subject of the following analysis. 


4 


2.2. Types of Trajectories to- Be Considered - - • 

The following assumptions and examples will indicate that family 
of trajectories that will be considered in the later analysis of capac- 


It is assumed that there are no obstacles in the approach area, 
nor are there any constraints due to noise; therefore, there are no 
restrictions as to the type of approach paths. 

The case of a fast aircraft followed by a slow aircraft will be 
initially considered since, as noted in Section 2.1, this is the criti- 
cal sequence. 

Assume that a turning radius of an aircraft is zero, that is, its 

heading can be changed instantaneously. Further, assume that a fast 

aircraft is at Eg , the entry gate for slow aircraft , where t = 0 

(as shown in Fig. 11). A useful question is: where should the slow 

aircraft be to ensure that the interarrival time t_„ is as small as 

FS 

possible? Obviously, the two aircraft should be separated at a minimum 
distance of 5 (i.e., the slow aircraft should be somewhere on the 

circle of radius d with center in E ) , but the condition of minimum 
separation has to be maintained until and after the moment the fast 
aircraft reaches E g . Assuming that a slow aircraft is always flying 
toward E , Fig. 11 shows that to maintain a minimum separation, 6 , 
its location at t = 0 has to be somewhere on the portion of the circle 
marked with a heavy line (e.g., point A). For any position of a slow 
aircraft on the light line portion of the circle, the condition d(t) £ 6 
will not be satisfied for t > 0 , where d(t) is equal to the hori- 
zontal distance between two aircraft. The conclusion is, then, that a 


slow aircraft should approach the point E anywhere from the heavy 


J 


1 


,. ; .:*.: ., .,,- 4 . -.- • 
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• line portion of the circle; in other words, from the right side of the 

line, perpendicular to the extended runway centerline passing through 
the point E . 

b 

Consider a slow aircraft at point M and a fast one at Eg (when 
t = 0) , as shown in Fig. 12. What kind of a path should the slow air- 
craft follow between M and Eg ? If a slow aircraft must reach Eg 
as quickly as possible, the desirable path is the straight line MEg . 

On the other hand, if it is necessary for some reason to increase the 
time that the slow aircraft takes to reach E , a curved path between 
M and E is one alternative. Another alternative is to require the 
slow aircraft at t = 0 , to be at a position, N , this position 
being further away than M on the same straight line; and then to have 
the slow aircraft continue on the straight line to E . Consequently, 
it can be concluded that straight line entries to Eg are at least as 
good as any other family of trajectories. The straight line entries 
because of their simplicity will be exclusively used in further analy- 
sis, A general straight line approach path configuration is shown in 
Fig. 13. 

Before proceeding with analysis of this configuration it is neces- 
sary to establish arbitrarily that the index of aircraft type increases 

i 

with approach speed, i.e., if i > j then v_^ > v , where i and j 
are indices of aircraft types. As noted in Section 2.1, that Y^ > Yg , 
it follows that Y.^ 5 Y 2 ••• 5 Y^ , where the subscript 1 refers to 
that aircraft type which needs the shortest straight-line part of the 
final approach, i.e. , that aircraft type that can intersect the exten- 
sion of th< runway centerline closest to the runway threshold. 

A useful question, then, in terms of understanding the character- 
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istics of the family of trajectories that is being considered is, what 

should be the angles a, , a„ , • • • a ? 

12 n 

It has been shown earlier in this section that la. I < 90° for 

' i 1 

any i , a.^ being positive when measured coni' terclockwise from the 

extended runway centerline and negative when measured clockwise . 

All approaches are assumed to be in a horizontal plane. Therefore, 
to assure that the two approach paths will not intersect, the following 
conditions must be met: 


a. < a . for 
i 3 

all i 

and 

3 

such 

that 

i > 

3 » 

a. > 0 

X 

and 

a. > 0 
3 

a. > a. for 
i 3 

all i 

and 

3 

such 

that 

i > 

j > 

a. < 0 
i 

and 

a . < 0 . 
3 


These conditions state that If i and j are a pair of aircraft 

such that v. > v, , then when both a. and a. are positive, a. 

x j i 3 x 

has to be less than a. . However, if a. and a, are both negative, 

J i 3 

then 0^ must be greater than a_. (i.e., the absolute value of 

is smaller than that of a.). 

3 

The above conditions enable us to use any possible combination of 
paths and compute t„ , where i can be F or S and j can be F 
or S , thereby preventing any two paths to intersect before they merge 
on the extended runway centerline. 

If the case of only two aircraft types, fast and slow, is consid- 
ered, then the necessary geometry of approach paths is shown in Fig. 14. 
This geometry requires that the following conditions have to be satisfied. 

0 < a g < 90° 

-90° < < a g . 
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At the beginning of this section an assumption, "zero turning 
radii," was introduced. If this assumption were not made the geometry 
of the approach paths would follow that shown in Fig. 15, rather than 
that shown in Fig. 14. We will now determine whether the geometry shown 
in Fig. 14 is a close enough approximation of the geometry shown in 
Fig. 15. 

For the purpose of this determination, even for very unusual con- 
ditions (high approach speed, high angle of interception of extended 
runway centerline, and low rate of change of heading) the straight line 
paths shown in Pig. 14 give a good approximation of the circular arc 
paths shown in Fig. 15. For example, if the approach speed of an air- 
craft is 160 kts and the rate of change of heading is 3° per second, 
the difference in the length of the approach path measured along the 
straight line paths and the circular arc is 0.09 nm. If the angle of 
interception a = 60° , and only 0.02 nm if the angle of interception 
a = 30° . These differences in the length of the approach paths are 
for the purposes of this study negligible. 

The conclusion, then, is that the geometry of approach paths that 
should be considered in further analysis of capacity is that shown in 
Fig. 14 for the case of two aircraft types, or that shown in Fig. 13 
for the more general case . 



Consider the case of the two aircraft shown in Fig. 16 . The thresh- 
old interarrival time, 


A 



I 


where . is increasing with initial separation _^_.d Q * Initial sepa- 

A 

ration . .d .is the distance between two aircraft i and j at the 
time t = 0 (leading aircraft i enters the common path on the extended 
runway centerline) , measured along the trajectory of trailing aircraft j . 
If for the given situation (aircraft types i and j and a given ap- 
proach path geometry) cne wants to minimize t , then the problem 
could be stated as follows: 


min . .d 
ij o 


subject to d..(t) - 6.. for t < 

J ij 13 v. 


where 


d..(t) = distance between aircraft i and j measured on a 

ij 

straight line in the horizontal plane 
6 „ = minimum horizontal separation for aircraft i fol- 

lowed by aircraft j , given by ATC rules. 


Condition t 1 states that the ATC horizontal separation, <5 . . , is 

v i 13 

required only while both aircraft i and j are airborne. This condi- 
tion can be released after leading aircraft, i , crosses the runway 

Y- . 

threshold, at t = — . However, the overall objective function of 

i 

the horizontal separation model is much more complex; it can be stated 
as follows : 


mxn t = 


E t . . p . . 
ij « ^ 


1 1 • 

subject to d. .(t) > <S . . , for t < — , for all i and j 

13 xj v. ’ 
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In this situation the geometry is not given, as it was in the pre- 
vious very simple example. The problem becomes, then, one of finding a 
geometry which for a given aircraft population and given ATC rules will 
provide the highest landing capacity. 

The next step is to study how the distance between two successively 
landing aircraft changes with time. 

2.4. Equations for Distance between Two Aircraft in a Plane 

If the selected geometry of approach paths is that shown in Fig. 14, 
and if two types of aircraft, fast F and slow S , are considered, 
there are four possible approach sequences to be studied: FF , SS , 

FS and SF (where FS refers to fast followed by slow, etc.). It 
has been noted in Section 2.3 that for the sequence of any two aircraft 
ij , t - 0 is the instant in time when the leading aircraft, i , 
enters the first point of the common path on the extended runway center- 
line . 

All four sequences mentioned above are shown in Fig. 17 at the 

point when t = 0 . As can be seen, for the sequence SF there are 

two cases: at the point when t = 0 and the aircraft S is ’t Eg , 

aircraft F could already be on the extended runway centerline, i.e., 

aircraft F could have passed through point E (case a) or it could 

still be on that part of the approach path leading towards E (case b) . 

/\ 

The important point here is that in all cases _d Q * s the dis- 
tance between two aircraft when t — 0 , measured along the path of 
the trailing aircraft, j . 

The next step is to analyze interarrival times by studying the 
distance between two aircraft measured along a straight line in a 
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horizontal plane, as a function of time. This distance will be denoted 
by d^(t) and will be studied for the four sequences mentioned above. 
However, as Fig. 17 illustrates, cases 1 and 2 (SS and FF sequences) 
represent essentially similar cases; they will therefore be studied as 
one case. 

2.4.1, The case of fast aircraft followed by a fast aircraft, or slow 
aircraft followed by a slow aircraft 

From Fig. 18 it can be seen that the function d(t) , the horizon- 
tal distance between two aircraft, is a piecewise function of t . 
Therefore , 

1. for t < 0 (neither of the two aircraft is on the extended 

centerline yet) 

d(t) = x d(t) = d Q ; 

A 

d 

2. for 0 < t < — (only the leading aircraft is on the extended 

runway centerline) 

d 2 (t) = 2 d 2 (t)= {vt + (d Q - vt)cos a} 2 + {(d Q - vt)sin a} 2 

= 2v 2 t 2 (l - cos a) - 2d vt(l - cos a) + d 2 ; (1) 

o o 

A 

d 

3. for t > — (both aircraft are on the extended runway 

v 

centerline) 
d(t) = 3 d(t) = d Q 

jV . . 

Only equations which are referred to later are enumerated. 
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2.4.2. The case of fast aircraft followed by slow aircraft 

From Fig. 19 it can be seen that the function d (t) , the hori- 

zontal distance between two aircraft (fast followed by slow) , is also 

a piecewise function of t . Therefore, 

1. for t < - ^ — (neither of the two aircraft is on the extended 
V F 

runway centerline yet) 
d FS (t) = l d FS^ 

= { ( FS d Q “ yv F t)cos a g - 3 + (v^t + 3) cos a p } 2 
+ {( FS d Q - ViVpt)sin a g + ( V]T t + 3) sin a p } 2 

= t 2 v 2 {y 2 - 2 cos(a g - a^) + l} 

+ t2Vp(3( (1 - cos a ? ) - |i(cos(a g - a p ) - cos a g )} 

a a 2 

+ FS d o {c ° s(a S " V " W + FS d o 

^ 2 
+ 23 FS d o (cos(a s - a ? ) - cos a g } + 23 (1 - cos a^,) 

( 2 ) 

8 ps d o 

— < t < •— — (only the leading aircraft is on the 
V F V F 

extended runway centerline) 


2. For 
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d FS^ 2 d FS^ t ^ 


a 2 A 2 

[( FS d o “ VV )co8 a * + v * t] + [( ^ d . " yVr)sin a q ] 


‘S F 


FS o H F 


t 2 v 2 (y 2 - 2ycos a g + 1) 


+ t2 FS d o V F <COS a S " ^ 


+ FSo > 


(3) 


3. For t > 


FS d o 

y 


(both aircraft are on the extended runway 
centerline) 


d FS 3 d FS^ C ^ 


FS d o + V 1 ' y) C • 


(4) 


2.4.3. The case of slow aircraft followed by fast aircraft 
Here two cases can be distinguished: 

a. at t = 0 trailing aircraft is already on the extended 
runway centerline, as shown in Fig, 20a (gp^ 0 < 3) > 

b. at t =? 0 trailing aircraft is not yet on the extended 
runway centerline, as shown in Fig. 20b (gj,^ > $) . 

In both cases, d„ (t) is the horizontal distance between the two 

DC 

aircraft and is a piecewise function of t . The two cases will be 
treated separately. 
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Case a 


1 


(See Fig. 20a, gF d Q < 3 .) 


For t < 


SF d o ' ^ 


(neither of the two aircraft is on the 
extended runway centerline yet) 


d SF (t) l d SF (t) 


= [v F t(ycos a R - cos a v ) 4 ' Q17 d n - 8) cos a r 4- 0] 


F SF o 


4 [v F t(ysin ct g - sin a p ) 4 ( gp d o - 0)sin a p ]' 


t 2 v F f.y 2 - 2ycas(a g - 0^) 4 1] 


4 t2v (0[(1 - cos a ) - y(cos(a - a ) - cos a )] 


+ SF d o [yCOs(a « - %> - H) + 


S F 


SF o 


- 28 gF d o (l - cos a f ) 4 23 (1 - cos a F > . (5) 


SF d o " ^ 

For < t < 0 (only the trailing aircraft is on the 


extended runway centerline) 


d SF^ 2 d SF^ t ‘ > 


^2 2 

- [v F t(ycos a - 1) 4 ^dj 4 [yv^t sin a g ] 


SF o J 


t 2 v 2 (] 2 -2]ioosr-,t- 1) 4 t2v , d ( 1 “| icos(x t , ) 4 d 2 . 
I' !> I* S I 1 o b or o 


( 6 ) 


29 


3. For t > 0 (both aircraft are on the extended runway 
centerline) 

d SF^ = 3 d SF < ' t ' ) 

= SF d o " V F (1 - u)t • <7) 

A 

Case b. (See Fig. 20b, „„d. > 6 -) 

' 1 ' b t O 

1. For t < 0 (neither of the two aircraft is on the extended 
runway centerline yet) 

4 (t > - l d SF (t > 

A 2 

= [v F t( ycos a g - cos a F ) + ( gF d o - 3) cos a p + 3] 

A 2 

+ [v F t(hsin a g - sin a ? ) + (g^ - 3) sin a p ] 

= t 2 v 2 [y 2 - 2vicos(a s - a F ) + l] Hr 

+ t2v F (3[(l - cos a F ) - vi(cos(a g - a f ) - cos a g )3 
+ sr d o [nc OS (a s - V - 1]) + SF d 7 
- 23 SF d o (l - cos a F ) + 23 2 (1 - cos a p ) . (8) 

Note that this ^d gF (t) is the same as in the previous 


case a. 


2. For 0 < 
d SF^ 


3. For t > 
d SF^ 
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t < 


SF d o " 3 


(only the leading aircraft is on the 
extended runway centerline) 


- 2 d Sf (t) 

= [v p t(u - cos otp) + ( S pd o - 6) COS a F + 3] 2 
+ [-v F t sin a F + ( SF d o - 3) sin a F ] 2 
= t 2 v 2 [vi 2 - 2ycos a F + 1] 


A 

+ t 2v ? [ SF d o (l - ycos a p ) - 3(1 - cos a F )(l + y)] 


SF o 


d 2 - 2 


d 3(1 - cos a_) +23 (1 - cos a ) . (9) 


SF o 



(both aircraft are on the extended 
runway centerline) 


3 d SF^ C - 

’ SF d o ' V F (1 ' v) C • (10) 
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2.5. Initial Separation 

This section will describe procedures for determining those optimal 
initial separations which satisfy model requirements. 

Initial separation here refers to that distance between two runway 
approaching aircraft, measured along the trajectory of the trailing 
aircraft, that exists when the first aircraft intersects the common ap- 
proach path (see Fig. 16). 

Section 2.3 notes that the threshold interarrival time is an in- 
creasing function of initial separation. This section also notes that 
the objective of the model is to minimize initial separation subject to 
constraints imposed by ATC rules on the required horizontal distance 
between two approaching aircraft. 

As noted in Section 2.4 the horizontal distance between two air- 
craft is a function of many variables rather than simply a function of 
time. The function could be written as follows. 


d..(t) = d. . (v. ,v. ,T. ,Y. ,0t. , a . ,t, d ) . 

3-1 11 1* 1* i 1* X* 1* *±j O' 


However, v, . v. , Y. and Y can be considered as fixed parameters 
ill j 


rather than variables , so 


d, . (t) - d . (a. ,a , t, . .d ) . 

ij il i j ij ° 


To 


maximize the capacity for a given geometry (i.e., ou , cu , 


and y. are fixed) one wants . . d to be as short as possible but still 
T ! il o 

such that the ATC rule, d,.(t) £ 6 . . , is not violated for any 

Y.. 1J „ 

t < — 1 i . For this condition and for , d to be minimized, the two air- 
y il 0 

craft have to be separated exactly when they are closest to each 


other. 


i 


I 


Using a partial derivative of d. . with respect to t , the 

xj 

moment of minimal separation t can be obtained, as follows. 


t - 


At the t moment the two aircraft should be separated by horizon- 
tal separation <$ , i.e., 

d..(t) „ = 6.. 

13 t=t* xj 


d..(t) - £(a.,a.,..d ) = 6, 

xi t=t* 1 1 in o 3 


The initial separation . .d is then found: 

ij o 


ij d o - K.a.,6^) . 


Note that (L can also be taken for a parameter with a fixed 
value; initial separation could then be written; 


. .d = £ (a. ,a.) 

1 ] o * 13 


As noted, the threshold interarrival time, t.. , is a function 

xj 

A 

of initial separation, ^_.d Q » this function will vary depending on 
the particular pair of aircraft types to be considered, so that 


f (a. ,a .) for i = j (two aircraft of same type) 

X- j , 


,a.) for 1 > i (fast followed by slow) 
I'b i 1 


(a ,a_.) for 1 < j (slow followed by fast) 

v 1* JL 1 


I 


l 


1 


) 


The functions t are included in the objective function of the model. 


min 2 t . . p . , 
13 3-3 


(p^ are given) 


s.t. d. . (t) > 6 . . 

3-3 3 .3 


for t 


< In 


for all i and j . 


The objective function, Et p^ » ^ ecom es very complex for the fol- 
lowing reasons. 

2 

a. Summation £t^_.p^. ^ as n members, when n represents 
number of aircraft types in airport population. 

b. Members are nonlinear functions of a. and a. . 

3- 3 

c. Each member belongs to one of the three types of functions 
mentioned above 


f (a. ,a.) for i - j 
3- 3 


f FS (a i’°0 for i > j 


f SF ( a i> a j) for 1 < 3 


However, function f can take two different forms, function f four 

r b 

and function f even five different forms, depending on the values 

b F 

of input parameters v^ , v j » Y^ » Yj and <5^.. ■> as well as on 

the values of variables a. and CL. themselves. The several possible 

3- 3 

forms of the three functions are the result of the fact that <$..(t) , 

. ... ' . 2-3 

being a piecewise function of t , could have a global minimum in more 


M r. ■ i - J i i » * * 
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than one of its segments. For example, the two aircraft could be mini- 
mally separated with neither being on the extended runway centerline, 
or minimal separation could occur when the first aircraft has already 
entered the extended runway centerline, or, when both aircraft are on 
the extended runway centerline. 

Because, as noted, the objective function is complex and because 
the constraints are non-linear, it was found necessary to compute initial 
separations and, consequently, interarrival times and runway landing 
capacity for fixed Ok and a_. . This computation involves the devel- 
opment of computer programs capable of determining capacities fpr dif- 
ferent discrete values of these intercept angles, 

2.5.1. The case of fast aircraft followed by a fast aircraft, or 
slow aircraft followed by a slow aircraft 

This case is shown in Fig. 18 and distance equations are given in 

Section 2.4.1. 

From equations for ^d(t) , 2^(0 gd(t) it can be seen that 

the min d(t) , i.e., the minimum separation between two aircraft, will 

A 

occur for 0 < t < d^/v , i.e., the minimum will be somewhere on the 

A 

segment 2 d(t) faction d(t) . Range (0, d Q /v) for variable t 

represents that time when the leading aircraft is already on the extended 
runway centerline and when the trailing aircraft is still approaching it. 

From equation (1) and for 

8 0 d 2 (t) 


the time when minimum separation occurs can be found by 
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At that time when they are minimally separated the two aircraft have 
to be separated distance 6 , 


>d 2 (t) 


= <5‘ 


t=t r 


i.e., if t^ is substituted for t in equation (1) and the distance 
is 6 , the following result is obtained: 



'/ 


f<i 


+ cos ) 



Subscript 02 indicates that the initial separation is corap., ited for the 
case when min d(t) occurs on the second segment, 2^(0 > ° f func- 

tion d(t) . 

The graph of function d(t) is shown in Fig. 21. As can be seen, 

the initial gap (horizontal separation) between the two aircraft begins 

to close when t > 0 and reaches its minimum, which is equal to $ , 

when t = t^ . The two aircraft are, therefore, always separated more 

* 

than is required, except at the moment t 2 . It is possible,, however, 

* 

for the first aircraft to land(e.g., very short Y ) before t = t^ 

(see Fig. 22); in this case, the separation of the two aircraft should 

be distance 6 at the moment when the first aircraft lands , if the 

conditions d. . (t) > 5. , is to be satisfied for that period when both 
ij ij 

aircraft are airborne, i.e,, 
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,d 2 (t) 


-= X 


6 2 . 


Therefore, from (1) the fallowing quadratic equation is obtained 


Ad 05 + Bd 05 + C 


where 


B ^ - 2 y(l - cos a) 

C = 2 y 2 (1 - cos a) - 6 2 . 
Consequently, the necessary initial separation is 


-B + / B 2 - 


4C 


05 


Finally, the necessary initial separation, d^ , for two aircraft 
of the same type, depending on which of the two cases shown in Fig. 21 

A A 

and Fig. 22 one deals with, must be either d^ or dg,. . 

Subroutines SDOFF and SDOSS of the program ..compute values for 

A A 

„„d and -^d . These subroutines are given in Appendix B.2. 

FF o SS o rr 


2.5.2. The case of fast aircraft followed by slow aircraft 

The above case is shown in Fig. 19, its distance equations are 
given in Section 2.4.2. From these distance equations it can be seen 
that the first two segments of the functions d^ s (t) , ^d^Ct) an d 
„d (t) are square roots of a parabola and, being distances, they are 
positive square roots of a parabola. The third segment of d^g(t) — 

_d c (t) — is a linear increasing function of t . It can be shown 

J Jc b 
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further that the slope of _d (t) for t = d /yv_ is always greater 
, Z r b r b o r 

than the slope of „d (t) . A conclusion can therefore be drawn: the 

J r b 

function d„ c (t) could resemble any of the various shapes shown in 
r b 

Fig. 23a; the following table lists all of the possible shapes of this 
function. 

l d F s(t> 2 d FS^ mln d FS^ at point 


c 

d 


a 

b 

b 


A 

A 

C 

B 

C 


Whichever of the above forms of the function d fS (t) is encoun- 
tered the same approach is taken, that is, when the two aircraft are 
separated by a minimum distance, that distance should be equal to the 
minimum horizontal separation required by ATC rules. This minimum 
separation condition is therefore used to find the necessary initial 
separation. 

However, there is another case to be considered in addition to 
those shown in the above table: the case presented in Fig. 23b. If 

the minimum distance between two aircraft occurs during the second 
segment of d Fg (t) , i.e., at 2 d Fg (t) , and if the leading aircraft 

lands before the gap between two aircraft reaches its minimum, i.e.. 


Yr 


FS 2 
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then, as the argument in Section 2.5.1 suggests, the minimum horizontal 
separation 5 should be imposed when t = y /v . In short, the 
function d Fg (t) reaches its minimum at the upper limit of its domain. 

Following are solutions of initial separation derived for the cases 
when min d (t) appears in one of the four points. A, B, C, or D 

rb 

shown in Fig. 23. 

Case A: min d-^t) at segment ..d-^t) , i.e., for t < - $/v_, 

r b x r b r 


From equation (2) and for 


3t 


= 0 , 


the time when the minimum separation appears, t , can be found: 

Fb X 


* _ FS“o l 
FS C 1 


d n [cos(a s -a F ) - y] + (1 - cos a v ) - y[cos - cos « R ] } 


S f 


{y 2 - 2ycos (a g - a p ) + l} v p 


Using the condition that the minimal distance between two aircraft is 
exactly equal to the minimum required by ATC rules, i.e., 


l d FS (t > 


* 

^FS*! 


= 6 


FS ’ 


from equation (2) the following is found 


FS1 A FS d 01 + FS1 B FS d 01 + FS1 C 


0 , 


where 


[cos(a g -a F > - y] 


A = 1 _ 

FS1 X 2 


y - 2ycos(a c , - a„) + 1 
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FST 


B 


2g{ [cos(ot g - a^,) - cos a c ] 


[cos (otg-otp) - y ][(1 - cos ot^) - y(cos(a g -a F ) - cosa g ) ] 
2 

U ~ 2ycos(a g -ap) + 1 




and 


FS1 


3 (2(1 - cos olJ 

r 


[(1 - cos a F ) - yCcosCag-ctp) - cos a g )] 1 
- - - } - 6 


y - 2]icos (a g -a F > + 1 


FS 


The necessary initial separation for the A case is therefore, 


+ A - 2 


FS d 01 


“fsi b + ■ FS1 B ~ 4 FS1 A FS1 C 
2 fsi a 


Case B: min d F g(t) at the limit of segments d Fg (t) and 2 d FS ( t ) » 

i.e. , for t = -3/v 

F 

From equations (2) or (3) and for 


1°FS ( t > 


t= -3/v 2"FS 
F 


2 d FS ^ ^ 


— 5 

t= -3/v f " FS 


the following is found 


FS4 A FS4 d 04 + FS4 B FS4 d 04 + FS4 C 0 


FS4 


A = 1 


where 


40 


fs4 B = - 2 $(cos a s - y) 


FS4 


c = 3 2 (y 2 - 2 y cos a s + 1) - fi| g 


The necessary initial separation for Case B is therefore, 


/ 


FS d 04 


__.B + / . B - 4 _.C 

FS4 FS4 F34 


Case C: min d Fg (t) at segment * i,e *» for ^ < t < 

F 

From equation (3) and for 


FS o 


3 2 d FS (t > 

St 


= 0 


the time when minimum separation appears, pg t 2 » can f° unc ^ : 


* 

FS t 2 


FS d o (cOS a S ~ y) 
(y 2 - 2ycos a g + l) 


Using again the same condition as in Case A, i.e., 


2 d FS^ 


* 

t= FS t 2 


FS 


and equation (3) , it is found that the necessary initial separation for 


Case C is : 


a > 


i 


FS d 02 


(cos ou - y)' 


U -2 cos a + 1 

D 


Case D: min d (t) at, the segment d (t) but when t = Y c /v„ , 

r b Z r b or 

i.e., when leading aircraft is landing. 


From equation (3) and for 


2 d FS (t) 


t- V v F 


it follows that 


where 


A A + n a 4. n 

F-S5 FS5 05 FS5 FS5 05 FS5 


2 Y g (cos a g - y) 


Yg (y 2 * - 2 y cos a g + 1) - 6 2 g 


The necessary initial separation for Case D is therefore. 


FS d 05 


FS5 B + ' FS5 B " 4 FS5 C 


The above solutions to cases A, B, C, and D illustrate how to compute 


and 


The question now is which of 


FS d 01 ’ FS d 02 * FS d 04 


FS d 05 * 


these four values represents ^-d 

FS o 

The rule governing this choice will be as follows 


„_d 
FS o 


m £ n FS d olc * 


such that 


d FS^ " 6 FS 


for all t and 


d FS (t) 


t= FS t 


k* 


= 6 . 


FS 


Algorithm for is given in Appendix A.l . Subroutine SDOFS 

A 

which computes i s given in Appendix B.2. 

2.5.3. The case of slow aircraft followed by fast aircraft 

This case is shewn in Figs . 20a and 20b , and corresponding dis- 
tance equations are given in Section 2.4.3. From the distance equa- 
tions it can be seen that the first two segments of d g p(t) , ^d gF (t) 
and are positive square roots of a parabola and that the 

third segment, d c _(t) , is an increasing linear function of t . 

) or 

As noted in Section 2.4.3, this SF (slow followed by fast) se- 
quence consists of cases a and b ; they will be treated here sep- 
arately. 

A 

Case a: < $ (See Fi g* 20a and Fig. 24a.) 

From equation (6) and for 
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A 

the time of minimum separation on second segment, SF t 2 » is f° un d Co 


SF d o (1 - ycos a s ) 

2 

(y - 2pcos a Q + l)v T 


It can be seen that sp t 2 > 0 ; the second segment of d gF (t) is 
therefore a decreasing function of t over its entire range (see Fig. 
24a) . 

On the other hand, the third segment of d gF (t) is a decreasing 

linear function of t . It can also be shown that the slope of 2 d gJ1 (t) 

when t = 0 (t = 0 being the limit between 2 d gF (t) and 3 d gF (t)) is 

steeper than the slope of 3 d gF (t) * 

The above characteristics of the function d gF (t) can t^ en be 

said to support the following conclusion: the shape of function d gF (t) 

could correspond to one of the alternatives shown in Fig. 24a. 

The following table lists all of the possible shapes of function 

d (t) , see also Fig. 24a. 

Si? 


ld gF (t) 


2 d SF ( ' t ' ) 


min d gF (t) 
at point 


Using, then, the argument developed in Section 2.5,2, the solutions 
for initial separations are derived for the cases when min d gF (t) 
occurs at points A or D (Fig. 24a) . 
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Case A: min d gF (t) at segment jd gF (t) , i.e., for t 


SFo 


d„ - g 


From equation (5) and for 


9 l d SF ( t> 


- 0 . 


3t 


The time when the minimum separation appears , 0 -pt , is then found to 

or i. 

be 


* _ SF 
S¥ t l 


d o [pcos(a s - a F ) - 1] + $[(1 - cos a^) - y(cos(a s -a F ) - cos a s )] 
(y 2 - 2ycos(a s - a p ) + l)v F 


Further, from condition 


l d SF (t) 


t= t* 
C SF 1 


SF 


and from equation (5) , the following is found 


SF1 A SF d 01 + SF1 B SF d 01 + SF1° ° 


where 


[y cos(a s - a F ) - 1} 


A = 1 - — 

SF1 2 


y - 2ycos(a g - a ? ) + 1 


SF1 B = - 23 {(l~cosa F > + 


[ycos(a g -a F )-l] [(l-cosa F >-y(cos(cx s -a F )-cosag)] 

' ' 2 ’ " — — - ) 

y - 2ycos (a g - a F ) + 1 
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SF1 C = 3 (2(l-cosa F ) - 


(1 - cosa F ) - ]i(cos(a g -a F ) - cosa g ) 
2 

U - 2 y cos(a g - a^) + 1 


The necessary initial separation for Case A is therefore, 


SF d 01 


SFl^ + * sfi b ~ 4 sfi a sfi c 
2 sfi a 


Case D: min cL_(t) at the upper limit of the range of segment 0 d (t) , 

b l 1 j br 

Yg 

i.e., for t t_ = (the time when leading aircraft is 


'5 yv T 


landing) 


From equation (7) and for 


t = t c 


* _ J_S_ 

SF t 3 ” yv 


and using the condition 


3 d SF^ * 

t= t ■ = t 
5 SF 3 


it follows that the necessary initial separation for case D is: 


SF d 03 SF d 05 5 + (1 “ y) y 
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Case b. ^d > 3 (See Fig. 20b and Fig. 24b, c, and d.) 

— d r O 

The only general conclusion about function d CT ,(t) is that in 

bi? 

case b the third segment, ^dg F (t) » is a decreasing linear function 
of t . 

The possible alternatives for the shape of d (t) and the loca- 

br 

tions of the min dg^Ct) are shown in Figs. 24b, c and d; these alter- 
natives are presented in the following table. 


l d SF (C> 

2 d SF (t) 

3 d SF (t) 

min dg F (t) 
at point 

Fig. 

a 

e 

f or g 

A 

24b 

b 

e 

g 

D 

11 

c 

e 

g 

D 

11 

c 

e 

f 

C 

If 

d 

e 

f 

C 

II 

d 

e 

g 

D 

ft 

a 

h 

g 

A 

24c 

b 

h 

g 

D 

If 

d 

h 

g 

D 

II 

a 

i 

f 

A 

24d 

a 

i 

g 

B 

11 

d 

i 

f 

B 

II 

d 

i 

g 

D 

M 


Comparing equations (5) and (7) or (8) and (10) it can be seen 
that they are identical, but the ranges of functions are different. 


However, the initial separations found for the cases when min dg^(t) 


i 



3 

■ $ 
f 


i 

?• 

$ 


■k 

. 4 *: 


$ 





l 


l 
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occurs at points A or D are the same , for both the a and b cases . 
Therefore, only the initial separations for the situations when min d gF (t) 
occurs at points B or C have to be found. 

Case B: min d gF (t) at the limit of segments ^^(t) and 2 d SF^ * 

i.e., when t = 0 

From equations (8) and (9) and for 


l d SF (t) t=0 


2 d SF (t) 


t-0 


SF 


the following is found 


SF4 A SF d 04 + SF4 B SF d 04 + SF4 C 


0 , 


where 


SF4 


A = 1 


sf4 e = - 2 6(1 - cos a F ) 

SF4 C = 2B 2 (1 - COS dj) - 6 2 f 


The necessary initial separation for case B is therefore. 


SF d 04 


SF4 B + ' 4 SF4 C 


2 


I 


J 


t 


Case C: min d gF (t) at segment 2 d SF (t) * i,e, » wlien 0 < t < 


SF d o”^ 


From equation (9) and for 


3 2 d SF^ 


the time when minimum separation appears, _ t. * can found, as 

br Z 

follows 


* S F d c U - ycos a p ) - 0(1 - cos « F )(1.+ y) 
SF t 2 = 71 “ , 77 

(y - 2 y cos a p + 1) v p 


Therefore from equation (9) and for 


2 d SF^ * 

t "SF t 2 


the following is found 


SF2 A SF d 02 + SF2 B SF d 02 + SF2° 


where 


(1 - ycos a p ) 

~~2 

y - 2ycos a + 1 
F 


(1 - ycos a ) (1 - cos a p ) (1 + y) 

-2{3{(1 - cos ou) - 2 ) 

y - 2ycos a +1 
F 
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SF2 C 


2 (1 - cos a_) (1 + p) 

g {2(1 - cos Op) ~2 } - <5 

p - 2 y cos a + 1 

F 


The necessary initial separation for case C is therefore, 


SF d 02 


SF2 B + >/ SF2 B " 


* SF2 A SF2 C 


SF2 


The above solutions to cases A , B , C , and D illustrate 

A A A A 

how to compute SF d Q1 , S pd Q2 , S pd Q ^ and S pd Q5 . The question now 

A 

is which of these four values represents gpd Q • 

The rule governing this choice will be as follows. 


such that 


SF d o = “k 11 SF d ok ’ 


d SF (t) - 6 SF £or a11 11 


and 


d SF^ t ^ 


t= S¥ t 


k* 


SF * 


The algorithm to find „d is shown in a flowchart in Appendix A. 2, 

SF o 

A 

Subroutine SDOSF which computes OT1 d is given in Appendix B.2. 

or o 

2.6. Equations for Interarrival Time at Threshold 

When all initial separations for the four sequences (FF, SS, FS 
and SF) are found, the threshold interarrival times for the same se- 
quences should be computed. 

Using the expression for t. . given in section 2.3, 




.1 - .. 


i a 


iiO 


c iJ 


h i. + ' iiS _ lii 


V. 

J 


V. 

X 


and the definition of y. . from Section 2.1. 

ij 


Y ±j “ min(Y 1 ,Y j ) , 


The interarrival times at the threshold for the four above sequences 
will be 


'FF 


FF d o 


SS d o 


'SS 


Y + d 
1 S FS o 


'FS 


'SF 


>Y S + SF d o 

v„ 


Y. 


Subroutine SDOT given in Appendix B.2 computes the values of 
interarrival times. 

In the case when there are more than two aircraft types in a 
population, the matrix of the interarrival times should be computed in 
such a way that the aircraft population is considered pair by pair; 
each pair always consists of a fast and a slow aircraft. 
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6 . Go to 1 * 

The restrictions for model input parameters as well as the ranges 
for variables OL and ot are given in Fig. 25. 

b r 

Program computing capacity for the case of only two aircraft types, 

CAPSF, is given in Appendix B.3. 

For the following input: 

-- aircraft approach velocities: v and v 

r b 

— necessary straight approach lengths: Y and Y c -I 

— proportion of aircraft types in the mix: p p and pg 

— matrix of horizontal separation rules: <$_ = ,6^ ’^SF’^SS^ * \ 
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the following output is obtained: 

— matrix of capacities A(a c »a ) 

b r 

A 

— matrices of initial separations sF^o^ a S ,a F^ 

FS S o ( V a F> 

FF^O ( a S ’ a F> 
SS 3 o<W 


— matrices of indices showing the shape of function d^(t) and 
location of minimum separation (see Appendix A) 


INDXSF(a g ,a F ) 
INDXFS(a g ,a F ) 
INDXFF(a s ,a F ) 
INDXSS(a g ,a F ) . 

An example of this output is given in Appendix C.l. 


2.8. Model for Arrival Runway Capacity When Aircraft Mix Consists 
of Three or Four Aircraft Types 

A procedure similar to that discussed in Section 2.7 is used to 
compute capacity for these two cases, a difference being that instead 
of preserving the values of capacity for all combinations of angles 
generated, the highest N values for capacity are stored during the 
computation and then printed, together with the values of the angles 
which provide those N maximum capacities. For the purpose of com- 
parison the lowest N values are preserved and printed as well. 

The procedure to determine N maximum and N minimum capacities 

■k 

for the case of three aircraft consists of the following steps. 

^ 

Note that the case when a population consists of four aircraft types is 
not illustrated, since this case is essentially similar to that fol- 
lowing three aircraft. 
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1. Generate angles and • 

2. Group the aircraft types in couples, i.e., 

1 and 2 

1 and 3 

2 and 3 , 

and treat each of these couples as a pair of slow and fast, 

S and F . 

3. For each of the couples compute initial separations, thus ob- 
taining a matrix of initial separations , 


ij d ° 


4. From the matrix of initial separation, d Q » compute the 

matrix of threshold interarrival times , t , . . 

il 

5. Compute mean interarrival time, 


t 


E 


t. . 
1-3 


P i 



6. Compute landing capacity, 

x =4- • 

t 

7. If the capacity found belongs to the set of N highest or 

N lowest capacity values, store it along with the correspond- 
ing angles. If the capacity does not belong to either of these 
sets, storage is unnecessary. 

8. Go to 1. 

The restrictions for model input parameters and the ranges of 


J - - - '• -- ■ ■ I I . i .... J \ 
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variables a. are given in Figs. 26 and 27, for the populations of 
three and four aircraft types. 

Programs computing capacities for the two populations, CAP 3 and 
CAP4 , are given in Appendices B.4 and B.5. 

For the following input, 

— aircraft velocities v for i = 1, 2, 3 or i - 1, 2, 3, 4 

— necessary straight 

approach lengths Y. " M 

— proportion of aircraft 

types in the mix P . " " 

matrix of horizontal 


separation rules 

d. . for i = 1, 2, 3 

or i - 

1, 2, 


and j = 1, 2, 3 

or j = 

1, 2, 


the following output is obtained: 

- N maximum values of capacity X and corresponding angles CL 
for i - 1, 2, 3, or i = 1, 2, 3, 4 

- N minimum values of capacity X and corresponding angles 
ct^ for i = 1, 2, 3 or i = 1, 2, 3, 4. 

An example of output from program CAP4 is. given in Appendix C.3. 
Value N is arbitrarily set at 40 . 
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3. Application of Capacity Models with Horizontal Separation Only 
3.1. Input Data 

Data for inputs into the models were obtained from: several ref- 

2 3 4 

erences given in the reference list ’ * , particularly those concerned 
with approach speeds and the minimum length of straight final approach; 
consultation with airlines and aircraft manufacturers; and approach 
speed studies at the TRACON facility at the Oakland International Airport. 

A review of these three data sources indicated that the following 
tabled inputs will usefully demonstrate the application of the proposed 
models , 



Aircraft Classification 

Aircraft Type 

Description 

Approach Velocity 

i 


v i (kts) 

1 

Propeller Driven 1 

100 

2 

Propeller Driven 2 

120 

3 

Nonheavy Jet 

140 

4 

Heavy Jet 

150 


Minimum Length of Straight Final Approach (y) 


Y • 

Aircraft Type x 

i (nm) 

1 2 

2 2 

3 , 4 


I 


1 


l 


l 


.ri- 


se 


Minimum horizontal separation rules for aircraft i followed by aircraft j 

a. = 5 mi when the leading aircraft is aircraft type 4 (i = 4) 

and trailing aircraft is either type 1, 2, or 3 (j = 1, 2, 3) 
= 4 nm when the leading aircraft is aircraft type 4 (i = 4) 

and trailing aircraft is aircraft type 4 (j = 4) 

= 3 nm when the leading aircraft is either aircraft type 1, 

2, or 3 (i = 1, 2, 3) 

b. 6. . = 3 nm for all aircraft types 

3-3 

c. 5.. = 2 nm for all aircraft types 


ij 


3-3 


3-3 


Aircraft Mix 

Aircraft Type percentage of aircraft type i in mix 


i 


a 

b 

c 

d 

e 

f 

s 

1 






30 

20 

10 

2 






30 

20 

10 

3 


80 

60 

40 

20 

20 

40 

60 

4 


20 

40 

60 

80 

20 

20 

20 

Using ILS, 

all 

aircraft 

have to 

pass 

through a 

common 

entry gate; 

therefore, the 

length of approach along the extended 

1 runway 

centerline 

is the same for 

all 

aircraft 

types ; 

that 

is, 

L “ 

y - 6 nm (6 nm is 


assumed to be the average distance from the entry gate to the runway 
threshold) . 

Using MLS, each aircraft of type i enters the extended runway 
centerline at its own entry gate, » at distance Y.^ from the 

runway threshold (see the table on the previous page) . 

For both ILS and MLS, the angles of entry to the common approach 
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path are optimized to provide maximum capacities. These capacities are 
then compared . 

As stated in Chapter 1, this comparison is made assuming that both 
systems are free of errors: aircraft precisely maintain assigned flight 

paths and arrive at points in space precisely when they are expected to 
be there . 

3.2 . Analysis of Results 

The first case involves aircraft mixes that contain only two air- 
craft types (mixes a, b, c and d) , namely, heavy and nonheavy jet air- 
planes. The results are shown in Fig. 28, and the output tables for 

Points A and B are given in Appendices C.l and C.2. 

The curves in Fig. 28 correspond to the optimum values of the 
angles of intersection of the entry path with the final approach path , 

a_ and ou . For point A, ot 0 = 10° and a„ = -10° (see table CAPSF, 

F o o r 

Appendix C.2) . The analysis also shows that a g and oc p can vary 

considerably and yet provide about the same capacity as long as the 

geometry of the flight paths is similar. 

There are several points that can be made about the angles a g 

and Op, which are important when landing capacity is considered (see 

Fig. 29). One is the value of the relative angle a (the angle 

K. 

between the two paths before they enter the extended runway centerline) . 

~ i a S "* 1 * 

This relative angle is important when the two aircraft have a common 
entry gate, as in ILS (Fig. 29c and d) or when they have different 
entry gates, as in MLS (Fig. 29a and b) . The optimum range of ot 
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v s 

depends on the ratio of approach speeds of two aircraft, y = — . 

V F 

The greater is y the smaller is a_ . From table CAPSF, in Appendix 

IV 

C.2, it can be seen that when v = 140 kts and when v = 150 kts , i.e., 

o F 

for a high value of y , the optimum range for a is 10-40° . 

R 

Another point about path configuration is how well centered a 

R 

is relative to the extended centerline of the runway. The configurations 
in which the approach paths are more or less symmetrical relative to the 
extended runway centerline, as shown in Figs. 29a and 29c, result in 
larger capacities than those when ot^ is not symmetrical with the ex- 
tended runway centerline (Figs . 29b and d) . This point can be seen in 
the CAPSF table, Appendix C.2. (The optimal range of angles for 
and 0L, is circled in that table.) 

Returning to Fig. 28, it can be seen that there is no essential 
difference in landing capacity between MLS and ILS when the mix consists 
of aircraft types with similar approach speeds, There is a slight in- 
crease in capacity as the proportion of faster aircraft in the mix is 
increased, when the minimum separation rules for all aircraft types are 
5 = 2 nm and 6 = 3 nm . This slight increase is to be expected be- 
cause the average speed of the stream of aircraft increases. When, 
however, the current 5 = 3, 4, 5 nm, separation rules are applied; 
the capacity decreases as the proportion of fast aircraft increases. 

This decrease occurs for the following reasons : because the fast air- 

craft are also the heaviest, therefore, when a heavy aircraft is fol- 
lowed by a light one the pair must be minimally separated a distance 
of 5 nm; the minimum separation must be 4 nm when a heavy aircraft is 
followed by another heavy aircraft. In short, the average spacing 
between approaching aircraft increases when the minimum separation rule 
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6 = 3, 4, 5 nm is applied, and this has more influence on the capacity 
than the mentioned increase of average speed of the aircraft stream. 

Fig. 30 shows the importance of optimizing entry angles. (The 

CAPSF tables in Appendices C.l and C.2 should be consulted for a more 

detailed treatment of this subject.) Note that the figure’s curves were 

plotted not only for the optimal angles of ±10° (see also Fig. 28 but 

also for values of a„ and a of ±40° , ±60° , and ±80° . De- 

b r 

creases in landing capacity of 4%, 11% and 21% were found relative to 
the optimal configuration. Similar results could be shown for cases 
when the minimum separation rules are 6 = 2 nm or 6 = 3, 4, or 5 nm. 

This implies (compare Fig. 28 and Fig. 30) that under the assumptions 
of the model and for the aircraft types considered, runway landing 
capacity is more sensitive to changes in entry angles than to other 
parameters, such as, for example, the proportion of fast and slow air- 
craft, or the length of the common approach along the extended runway 
centerline. As shown in Fig. 28, shortening of the common approach 
length from 6 nm to 4 nm does not result in any significant increase of 
landing capacity even for the optimal angles of entry, because the dif- 
ference in approach speed of the two aircraft is relatively small. 

A second case which demonstrates the application of the model 
involves four aircraft types in a mix (mixes e, f and g from Section 
3.1). The results are shown in Fig. 31. (Example output for point C in 
Fig. 31 is given in Appendix C.3.) From this figure it can be seen that 
the use of MLS increases landing capacity, and for the given example the 
increase is about 6-10% when minimum horizontal separation is 6 = 3, 4, 

5 nm, 7-11% when 5 = 3 nm, and 10-16% when 6 = 2 nm. 

Note that the increase of capacity is greater when the aircraft 
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mix consists of approximately half fast and half slow aircraft (e.g., 

P 1 = p 2 = *3 » P 3 - - . 2 ) than in the case when one type of air- 
craft prevails (e.g., = p^ = .1 , p^ - .6 , p^ = .2 which is a mix 

with 80% of fast aircraft) . This point agrees with the conclusion made 
earlier in this section that the increase of capacity obtainable with MLS 
is .rather small when a population consists of aircraft types with simi- 
lar approach speeds. This point also agrees with the discussion of ex- 
pected capacity increase given in Section 2.1. 

The capacity increase obtained in the second case (four aircraft 
types) is significantly greater than that obtained in the first (two air- 
craft types), for two reasons: first, the differences between the veloc- 

ities of the aircraft types are much greater in the second case than in 
the first; second, the decrease of the common path length is greater 
in the second case (on the average) , and particularly because the slow 
aircraft types enter the entended runway centerline only 2 nm from the 
runway threshold. 

As the model computes capacity for discrete values of angles , 
it should be noted again (as in the case when there are two aircraft 
types in the mix) that, for all practical purposes, a group of solu- 
tions (rather than a single solution) ensure maximum capacity. This 
can be seen from the output example for CAP4, Appendix €.3. 

Some of the solutions are presented graphically in Pig. 32; all 
result in about the same capacity. Differences between solutions are 
the result of; slight changes in some of the entry angles; an inter- 
change of the given trajectories for two similar aircraft types; or 
a combination of both. It should, however, also be noted that the 
flight path configurations shown on Fig. 32 are not radically different, 
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If entry angles radically different from those shown on Fig. 32 are 
used, the capacity could be reduced by as much as 20% (see the Minimum 
Capacity table in Appendix C.3). 

Some additional examples indicate that further decrease of the 
length of common straight approach (e.g., y^ = 1 nm, y^ = 2 nm, y^ = 

3 nm and y^ = 4 nm) will not bring any significant increase in runway 
landing capacity (see Fig. 31). 


1 




1 


I 


l 


7.1 t i 
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4. Analysis of Capacity Models with Horizontal and Vertical Separation 
4.1. Introduction 

Chapters 2 and 3 have considered a landing capacity model which is 
governed by the assumption that in accordance with current ATC rules only 
horizontal separation between aircraft is allowed during approach. 

In this (and the next) chapter, runway landing capacity will be 
analyzed. It will be here assumed that horizontal and vertical separa- 
tion between aircraft approaching a runway is possible. Landing capac- 
ity will be computed using, for each particular aircraft pair, the more 
efficient of the two possible aircraft separations. 

Some preparatory remarks should be made before proceeding. First, 
no two-segment approach paths will be considered. Second, as the fea- 
sible angles of descent under IFR conditions vary from 2.5° to 8-10°, 
this span of about 7° does not allow for more than two "useful" paths 
in the vertical plane, useful referring to those paths which provide 
vertical separations that are more suitable from a capacity point of 
view than the corresponding horizontal ones; consequently, there is no 
point in dividing all aircraft types into more than two categories , 
these categories being determined by the descent angle capabilities of 
the aircraft type. 

These two categories will be: (1) aircraft capable of performing 

steep descent (STOL aircraft) and (2) the remaining aircraft population 
(GTOL aircraft) . Second category aircraft are assumed to be in the 
same horizontal plane and are separated by horizontal separation rules . 

If a vertical separation of only 1000 ft (305 m) were applied (for 
descent angles of 2.5° or 3°) the resulting horizontal spacings are 
either greater than the prescribed minimum given by the horizontal 
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separation rules or are infeasible because of wake vortex separations . 
The argument in chapters 2 and 3 will apply to the situation when both 
groups of aircraft appear in the mix. Vertical separation will be ex- 
amined only for the cases when at least one of the two consecutive land- 
ings is performed by an aircraft with steep descent capability. 


4.2. Vertical Separation 

If a pair of aircraft land one after the other and one of them is 
capable of performing steep descent, it is assumed that two distinct 
approach paths in the vertical plane are followed. (This situation is 
presented in Fig. 33.) These paths are denoted as higher approach, H , 
and lower approach, L . Similar notations will be used for all air- 
craft utilizing these paths. It is assumed further than the aircraft 
using the higher path will have lower approach speeds than those using . 
the lower path, i.e.. 


V H < V L * 


This assumption implies that if only horizontal separation is 
employed, the aircraft pair landing sequence low followed by high (LH) 
will be critical from a capacity point of view because it actually rep- 
resents the fast followed by slow aircraft sequence (FS) . This sequence 
is shown in Fig. 34 at the moment when the landing aircraft, L , crosses 
the threshold. If the vertical separation, X » between the two air- 
craft produces , the horizontal spacing between the two aircraft 

at that moment — which is less than the spacing required under ATC hori- 
zontal separation rules, some gains in capacity can be expected. How- 
ever, to assure that a vertical separation between the two aircraft of 


A 


J 
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at least x is continually preserved before the leading aircraft, L, 
lands, the following condition has to be satisfied: 

v sin 0 > v sin 0 , 

ri Li L 1 j 

i.e., the vertical component of the speed of aircraft H has to be 

equal to or higher than that of aircraft L . If this condition is 

satisfied, the vertical gap between the two aircraft continually closes 
during approach until it reaches x > the minimum required vertical 
separation, when the leading aircraft lands. 

This condition can also be: 

v sin 0 T 
S Li L 

v 

H - . n 

sin 0 

H 

This variation of the condition suggests, since angles 0 and 

L 

0 TT are small, that if the angle 0 is three times as big as angle 
H n 

0 , then v has to be greater than one third of v . This is an 

L rl 

approximation, but it shows that when 0 is large enough the condi- 
tion is always satisfied. 

An exact solution for the minimum speed of an aircraft on the 

upper approach path, when v^ = 150 kts, 0^ = 3° and 0^ = 7.5° , 

is that v > 60 kts. This condition will be satisfied because the 
H 

aircraft capable of steep descent have expected minimum approach speeds 
of 70 to 80 kts. 

When the necessary vertical separation condition is satisfied the 
separation between the two aircraft in sequence LH , at the moment 
when leading aircraft lands, will be 
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When 6 =7.5° , 


The table indicates that vertical separation causes some capacity 
improvement, since the horizontal spacings of vertically separated air- 
craft are shorter than the minimum required horizontal separations, 
except when <5„ = 2 nm and y ~ 2000 ft (610 m) . 

The sequence under consideration has been LH , an aircraft on a 
lower path followed by an aircraft on a higher one. However, there are 
three other possible same pair sequences: HR , 1IL , and LL . 

Sequence HH will require a vertical separation X when the lead- 
ing aircraft lands, at which moment (see Fig. 34) the horizontal spacing 

will be & 

H X 

For sequence HL , -orizontal separation will be such that when 
the leading aircraft, H , lands, the trailing aircraft, which has, 
as noted, a higher approach speed, will be distance 6 from the thresh- 
old. 

Sequence LL requires horizontal separation only. The model de- 
scribed in Chapter 2 will be used to compute threshold interarrival 
time as it relates to horizontal separation. 


I 

H X 


X_ 

tan 0 


H 


£ has the following values. 

A 


X 


„ £ 

H X 


(ft) 

(m) 

(nm) 

1000 

305 

1.25 

1500 

457 

1.87 

2000 

610 

2,50 


I 


I 


L 


i 


r 


The requirements for each of the above four approach sequences 
suggest that the threshold interarrival times for a pair of aircraft 
using different descent paths are 


t 'LH " t HH 


"HL v T 


t = to be found using the model for horizontal 
,LL» 

separation only. 


The next section will describe how these interarrival times can be 
combined and runway landing capacity will be computed. 


4.3. Model for Arrival Runway Capacity 


The following description of this simple model will consider four 
differing aircraft types using the same runway. The types are denoted 
as 1, 2, 3 and 4; speeds ascend from type 1, the slowest, to type 4, 
the fastest. Only type 1 aircraft are capable of performing steep 
descent, i.e., they are type H ; all other types are type L (see 
Section 4.2.). The following model parameters are 


= aircraft approach speed 
— min straight final approach 
lengths 

= proportion of aircraft type, 


= 1, 2, 3, 4, 


= 1, 2, 3, 4 
= 1, 2, 3, 4 


min horizontal separation matrix i = 1, 2, 3, 4; j = 1,2, 3, 4 


0^ = angle of descent 

X = minimum vertical separation. 


0 1 > 0 2 - 0 3 - S 4 


To find runway capacity when aircraft type 1 use a steeper angle 
of descent than the other types, and when vertical separation between 
aircraft type 1 and the others is applied, the following steps should 


be made. 


1. Check whether the necessary condition for vertical separation 
(given in Section 4.2) is satisfied, i.e.. 


v sin 0 T 
> L L 

sin 0„ 


in this case it should read: 


v. sin 0, 
> _4 4 

sin 0., 


and v. > v £ v and 0 = = 0 / • 

4 3 2 2 3 4 

If the condition is satisfied vertical separation can be 
applied. 

Form the matrix of interarrival times, as follows. 


Trailing aircraft 


t l,l t l»2 t l,3 t l,4 


Leading 

aircraft 



1 


l 


1 
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^1 1 -i ^<1 1 1 t_ 


0 

HI 


X 


'1,1 2,1 3,1 4,1 HR v v tan 6 v 1 tan 8 

H n HI X 


■ 1,2 


1,2 


'1,3 


'1,4 


1,3 


1,4 


'LL 


is computed as the mean interar rival time when the population 
consists of three aircraft types: computation involves use 

of a model which assumes horizontal separation only (Program 
CAP3) , All the input values of this model are the same as 
detailed earlier in this section, except that the proportion 
of aircraft in the mix is modified in such a way that 

X p! = 1 , where p! are new proportions and the 
i=2,3,4 r 1 

relation between , P 3 » and P^ is same as IP was 
between p^ , p^ and P^ , namely. 


P i 


P i 


1=2,3, 4 


P i 


It should be noted that when obtaining t^ L the optimal 
horizontal path geometry for aircraft types 2, 3 and 4 will 
be found. 

Find the mean threshold interarrival time, t . 





4. Compute runway landing capacity, which is 
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5 . Application of Capacity Model with Horizontal and Vertical 
Separation 

5.1. Input Data 

The second of the two examples given in Chapter 3 will be here 
modified so that the influence of vertical separation on capacity can 
be determined. 

The example's classification of aircraft types is the same as 
that given in Section 3.1, page 55, excepting that type 1 aircraft are 
considered capable of performing steep descent. A corresponding set 
of descending angles are assumed as follows. 


A/C type Angle of descent 



Another exception is that if type 1 aircraft perform a descent 
where 0^ = 7.5° they will have lower approach speeds , so that v^ - 
80 kts when vertical separation is applied. Speeds for aircraft types 
2, 3, and 4 are the same as those listed in Section 3.1, page 55. 

The remaining inputs for the capacity model with horizontal and 
vertical separation are as follows. 

— Minimum lengths of straight final approach for all four air- 
craft types are the same as those given in Section 3.1, page 55 , 
— Minimum horizontal separation rules are as given in Section 
3.1,- page 56. 
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-- Minimum vertical separation rules are assumed to be the same 
as given in Section 4.2: x = 1000, 1500 and 2000 ft (305, 457 

and 610 m) . 

— Aircraft mixes used in this example are the mixes e , f and 
g given in Section 3.1, page 56. 

In summary, modification of the selected example involves the fol- 
lowing assumptions: an aircraft type 1 performs descent at a 7.5° angle 

with an approach speed of 80 kts; type 1 aircraft are separated ver- 
tically; type 1 aircraft are also separated from other aircraft types 
by vertical separation; a vertical separation between any two aircraft 
types can be either 1000, 1500 or 2000 ft. Aircraft types 2, 3 and 4 
are .separated horizontally. 

Using the model described in Section 4.3, the results plotted in 
Figs. 35, 36 and 37 were obtained. Note that in these figures, plots 
for only horizontal separation are actually those given earlier in 
Fig. 31. 

5.2. Analysis of Results 

In Section 4.2, it was noted that for a steep angle of descent, 

8 = 7.5°; minimum vertical separations of y are 1000, 1500 and 

2000 ft (305, 457 and 610 m) ; horizontal spacings between two aircraft 
are 1.25, 1.87 and 2.50 nrc, respectively, at that moment when the lead- 
ing aircraft lands and trailing aircraft is on a steep descent path. 
These spacings are much shorter than the horizontal spacings between 
two aircraft at the moment when the leading aircraft lands if current 
ATC horizontal separation rules are applied. Consequently, significant 
improvements in landing capacity should be expected. An exception is 


-I 


] 


V 
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the case when required horizontal separation is only two nautical miles 

and vertical separation is either 1500 or 2000 ft; this exception 

applies particularly when vertical separation of 2000 ft is required 

and when 0 = 7.5°; the resulting horizontal spacing, l , equals 

n u X 

2.5 nm; this horizontal spacing is larger than that required by mini- 
mum horizontal separation rules; therefore, the use of vertical separa- 
tion in this case would decrease rather than increase capacity. The 
greatest improvement should he expected in the case when minimal hori- 
zontal separations are in their high range, e.g., 3, 4, 5 nm separation. 

Figs. 35, 36, and 37 shows runway landing capacities for ILS and 
MLS procedures, the latter involving the use of: only horizontal sepa- 

ration, and horizontal and vertical separation. Input values for these 
capacities are those listed in the previous section 

The single difference in the capacities shown in Figs. 35, 36, 
and 37 is a result of differing minimum horizontal separation rules . 

The minimum horizontal separation rules used in computation of the 
capacities shown in Fig. 35 require that 6 =■ 3, 4, or 5 nm. As this 
figure indicates, MLS procedures can produce a capacity increase of as 
much as 40 to 50% when compared under the same aircraft mix conditions 
with ILS procedures. 

Fig. 36 shows the capacities obtained when the horizontal separa- 
tion rules require that 5 = 3 nm. The increase in capacity when MLS 
rather than ILS procedures are used can be as much as 30 to 40%. 

The capacities shown in Figs* 35 and 36 indicate that when any of 
the three vertical separations (X - 1000, 1500, or 2000 ft) are applied 
higher capacity increases, compared with capacities obtainable when 
MLS procedures with only horizontal separations are applied, are 


obtained. However, when the horizontal separation rules are from a 
capacity point of view improved, i.e., minimum horizontal separations 
are shorter, the additional gain in capacity, if any, is, as noted, 
smaller, see for example Fig. 37. This figure shows that a vertical 
separation of 1500 ft produces about the same result as a horizontal 
separation of 2 um, i.e., the increase in capacity obtainable in both 
MLS cases (with and without vertical separation) is about 15%, when 
compared with the capacity obtainable with ILS. The figure also shows 
that a vertical separation of 10Q0 ft produces a capacity improvement 
of 30%, when compared with that obtainable with ILS. In contrast, the 
figure shows that if type 1 aircraft are separated vertically, from 
themselves and from other aircraft, by 2000 ft, rather than horizontally, 
by 2 nm, the effect of vertical separation on the capacity is negative. 

Note that the highest obtainable increase in landing capacity when 
using MLS procedures occurs when the aircraft population consists of 
roughly equal divisions of fast and slow aircraft. (This point is dis- 
cussed in Section 3.2.) 

In Section 4.3 it was noted that when the capacity model (which 
responds to both horizontal and vertical separation rules) is applied, 
it produces optimal horizontal path geometry for aircraft types 2, 3, 
and 4. Further, in Section 3.2 it was noted that when horizontal sepa- 
ration rules are applied (to either MLS or ILS procedures) more than 
one optimal configuration results. Fig. 38 illustrates one of these 
approach paths configurations; this configuration maximizes landing 
capacity and produces the results shown in Figs. 35, 36 and 37. The 
paths are as follows. 

— Aircraft type 1 approaches on a glide slope of 7.5°, performing 
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either a straight (as shown in Fig. 38) or curved approach. 
Aircraft types 2, 3 and 4 all approach on a 3° glide slope. 
Aircraft type 4 approach directly along the extended runway 
centerline. 

Aircraft type 3 intercept the extended runway centerline from 
its left side, at an angle of 20°, 4 nm from the threshold. 
Aircraft type 4 intercept the extended runway centerline from 
its right side, at an angle of 30°, 2 nm from the threshold. 
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6 . Conclusions 

MLS allows the use of multiple approach path geometry in the final 
stage of approach, whereas ILS requires use of a common approach path. 
This research has considered runway landing capacity increases obtain- 
able when MLS rather than ILS procedures are employed. Conclusions 
are as follows. 

1. The approach path configuration shown in Fig. 13 is at least 

as suitable from a capacity point of view as any other possible 
configuration, when only horizontal separation is allowed (see 
Section 2.2). 

2. The angles of entry, ou , to the extended runway centerline 
have a significant effect on landing capacity; they should 
be optimized (see Section 3.2). 

3. When only horizontal separation is allowed and when an air- 
craft population consists of aircraft with similar approach 
speeds, MLS procedures (differing final straight approach 
path) do not produce any improvement in landing capacity, 
when compared with capacities obtainable with ILS procedures 
(a common final approach path). (See Section 3.3.) 

4. When only horizontal separation is allowed and an aircraft 
population consists of aircraft with considerably differing 
approach speeds, and for mixes of roughly half fast and half 
slow aircraft landing capacity improvements of 10-15% can be 
expected, when MLS rather than ILS procedures are employed 
(see Section 3.2). 

5. If both vertical and horizontal separations are allowed during 
the approach phase, depending highly on the aircraft mix and 
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the applied minimum vertical separation rules, very significant 
improvements in landing capacity, as much as 50%, can be achieved, 
when MLS rather than ILS procedures are employed (see Section 
5.2). These improvements are due to the vertical separation. 
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2. v s <v F 

3. Pp + p s e I ; P F >0; p s >0 

4. Angle increments 10° 

except when ctgSsdp 

then a s -d p * 5° 
e.g. 
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Figure 25 
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Figure 26 
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b. 


20° < d| < 90° 
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-90° < d 3 < - 10° 
d 3 + 10 < ^4<d2"IO 


Figure 27a 
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Figure 29 
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Figure 37 
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8. Glossary 

Subscripts 

1 - aircraft type. This subscript can take numerical values 

1, 2, ...» n, where 1 is the aircraft with the lowest 
approach speed. It can also take values F or S, 
denoting fast or slow aircraft in the case of only two 
aircraft types in the population. 

ij - aircraft pair, consisting of two aircraft landing one 

after the other, i is the leading aircraft type; j 
is the trailing aircraft type, j can take all values 
that i takes (see the preceding paragraph). 


Capital Letters 

E - entry gate, point at which aircraft joins final straight 

approach path before landing. 

E.^ - entry gate, for aircraft of type i in the case of multiple 

approach paths. 

T - runway threshold. 

Small Letters 


d(t) 


d . . (t) 
ij 


iVo 

2 V 0 


- distance between two aircraft measured on a straight line 
in horizontal plane. This distance is a function of time. 

-distance d(t) , for a pair of consecutively landing air- 
craft (i is the leading aircraft type and j is the 
trailing aircraft type) . 

- first segment of function d^ (t) . 


- second segment of function d 


ij 



< T> 
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3 d ij (t) 


xj om 


third segment of function r?.^(t). 

initial separation between aircraft i and j. This is 
the separation between the two aircraft at time, t = 0 
(the moment when leading aircraft i reaches the begin- 
ning of the common straight approach path) , measured 
along the path of trailing aircraft. 

initial separation between aircraft i and j computed 
under the condition that the minimum of function d^ (t) 
is located as indicated by subscript m 


m = 1 min d^ (t) occurs at the first segment of the 

/\ 

function. ..d n _ is found from 
xj 01 


i d ij (t) 


t= t* 

ij 1 


6 ij 


m 


min djj (t) occurs at the second segment of the 

A 

function. . ,d rt „ is found from 
xj 02 


2 d ij ( t) 


t= t* 
ij 2 


- 6 .. 

ij 


m = 3 min d..(t) occurs at the third segment of the 
ij ■' 

A 

function. . .d_.„ is found from 
xj 03 


3 d lj (t> 


* 

t=. .t_ 
ij 3 


- 


m = 4 min d..(t) occurs at the limit between the first 
xj 

A ' ' . 

and second segment of the function, ^d^ xs 
found from 


i d tj (t) 


* 2 d ij (t) 


t= t. 

ij 4 


= 6 . . 

xj 


t= . t, 

ij A 
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m = 5 min (t) occurs at the upper limit of the domain 
of the function for t = ^t^ (i.e., leading air- 

A 

craft i lands) .,d._ is found, therefore, from 
ij 05 


d H (t> 




t=. .t- 
13 5 


13 


t. . 
13 


t. . 
a ij 


t. . 
r ij 


. .t* 

13 1 

t* 
ij 2 

t* 

ij 3 


proportion of aircraft type i in the population, 
probability of the pair ij aircraft occurring in the 
vehicle stream, i.e., aircraft type j landing after air- 
craft type i . 
time 

Note: t = 0 indicates for any pair of consecutively land- 

ing aircraft the time when leading aircraft reaches the 
first point on the common straight approach path (on the 
extended runway centerline) . 

interarrival time at the runway threshold (time separation 
over the threshold) between leading aircraft of type i 
and trailing aircraft of type j. 

time separation over the threshold dictated by the minimum 
separation rules in the air. 

time separation over the threshold dictated by the runway 
occupancy time. 

expected interarrival time at the threshold. 

time at which „d. . (t) reaches minimum. 

1 ij 

time at which _d .. (t) reaches minimum. 

2 ±3 . . 

time at which „d..(t) reaches minimum. 
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1 * 


t 2 * 

ij 1 

etc. 

t 5 * 
ij 3 


in general 

k* 

. .t 
ij 


time at which d (t) reaches minimum, computed for 

■ ” V 

/s /v 

. .d = . .d.. . 
xj o xj 01 

time at which ^d^.Ct) reac ^ es m inimum, computed for 

/N ^ 

..d = . .d „ . 

xj o xj 02 

t im e at which 3^^^^ reaches minimum, computed for 

A A 

. « d = , , d n _ • 
xj o ij 05 


time at which d . . (t) reaches minimum, computed for 

A A 

..d n « . .d . . 
xj 0 ij ok 


. .t. 
xj 4 


IJ 5 


- time limit between the first and the second segment 
(jd^j (t) and function d^(t) . 

time at which first of two aircraft in an aircraft pair 
lands 


y. . 

t.. = 

xj v. 


ij fc 6 


v. 


time limit between the second and the third segment 

(_d. . (t) and „d . . (t)) of function d..(t) . 

2 xj 3 ij xj 

- approach speed of aircraft of type i. 


Greek Letters 

ol - angle of entry of the aircraft of type i to the extended 

runway centerline measured relative to that line. Positive 
anti- clockwise , negative clockwise. 

a„ - relative angle of entry for two aircraft of different types 

K. 
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a, 


R 


a. - a. 
1 3 


i * 3 


e.g. 


a. 


R 


V - a F 


Y 


Yj 


Y. • 


6 . . 

13 


0 , 

l 


X 

u- 


difference in the necessary lengths of straight final 
approaches along the extended runway centerline for 
fast (F) and slow (S) aircraft 

3 = Y f -Y s 

length of the straight common final approach along the 
extended runway centerline, for ILS case, 
necessary length of straight final approach for aircraft 
type i. 

length of the common straight final approach for aircraft 
type i and type j , for MLS case 

Y. . = min(y. >Y- ) 

13 i 3 

minimum horizontal separation between aircraft type i 
followed by aircraft type j, required by ATC rules, 
angle of descent in the final approach of aircraft of 
type i. 

runway landing capacity. 

ratio of approach speeds of slow (S) and fast (F) air- 
craft. 


X 


M = — < 1 

v — 

F 


minimum vertical separation between two consecutively 
landing aircraft. . 
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Notes 

1. It can be shown by geometrical proofs that all of the exits from 
the algorithm shown as "ERROR” cannot occur, if the input values of 
the variables and parameters are correct. However, such exits are 
preserved in the programs to facilitate location of eventual errors 
in programming or errors in the input data. 

2. Index "INDX" describes the shape of tbe function d^gCt) and shows 
the range of this function where minimum separation between the two 
aircraft occurs. The key to values of "INDX" is given in the two 
following figures. 
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SF d o ~ SF d 03 


INDX - 3011 


d = d 
SF o SF 03 


INDX = 3021 


d = d 
SF o SF 01 


INDX =1021 


















Notes : 


1. 2*, 2" and 2 ,M are exactly the same as 2 but with changed error 

indices as follows : 


2 

2' 

2” 

2 m 

1001 

1005 

1013 

1020 

1002 

1006 

1014 

1021 

1003 

1007 

1015 

1022 

1004 

1008 

1016 

1023 


2. See note 1, Appendix A.l, page 136. 

3. See note 2, Appendix A.l, page 136. 
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Function d c (t) 

o r 
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All Subroutines 


COAS 

cos a g 

COAF 

cos a F 

COASAF 

cos(a s - a p ) 

AMAS 

2 

y - 2ycos a g + 1 

AMAF 

2 

y - 2Ucos + 1 

AMASAF 

2 

y - 2ycos(a g - 

AMI 

U 

VF 

V F 

VS 

V S 

GAMAF 

y f 

GAMAS 

Y S 

BETA 

3 

IAFAR 

a in degrees 

r 

IASAR 

a_ in degrees 

u 

A 

DOSF 

SF d o 

A 

DOFS 

__d 
FS o 

A 

DOSS 

SS d o 

A 

DOFF 

FF d o 

INDXFS 

INDEX for d Fg (t) 

INDXSF 

INDEX for d gF (t) 

INDXSS 

INDEX for dg S (t) 

INDXFF 

INDEX for d (t) 

x r 

TSF 

t SF 

TFS 

t FS 

TSS 

hs 
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IFF 

ft 

DELTSS 

fi SS 

DELTFF 

6 ff 

DELTSF 

6 sf 

DELTFS 

6 fs 

■ ■■■■■ ; '* 

Subroutine SDOSF 


A 

DOS FI 

SF d 01 


A 

DOSF2 

SF d 02 


A j ■ ■ 

D0SF3 

SF d 03 


a] 

DOSF4 

SF d 64 

visa 

t 1 * 
SF 1 

V2S1 

2* 

SF fc l 

V3S1 

t 3 * 
SF 1 

V4S1 

t 4 * 
SF 1 

V1S2 

t 1 * 

SF l 2 

V2S2 

^ 2 * 

SF~2 

V3S2 

3* 
SF 2 

DV2S1 

l d SF ( 

DV3S1 

l d SF^ 

DV1S2 

2 d SF^ 

DV3S2 

2 d SF* 

DV1S3 

3 d SF^ 

DV2S3 

3 d SF^ 

* 

Subroutines 

have i 


2*v 

t 3 *) 

f 1 ' 

& 

& 

t 1 *) 

r 3 1 

A*) 


as the first letter, e.g., SDOSF is subroutine which computes DOSF . 
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DV4S3 3*SV { SVH ) 

DODO 3 2 d SF^ computed using SI> dQ 3 

Subroutine DOFS 


A 


D0FS1 

FS d 01 

A 

DOFS 2 

FS d 02 

DOFS4 

A 

FS d 04 

DOFS5 

A 

FS d 05 

T1S1 

t 1 * 

FS 1 

T2S1 

i2* 

FS^ 

T4S1 

t 4 * 

FS 1 ; 

T5S1 

t 5 * 

FS 1 

T1S2 

t 1 * 
FS 2 

T2S2 

.2* 
FS 2 

T4S2 

t 4 * 

FS 2 

DT2S1 

1 FS-FS 1 } 

DT5S1 

5* 

l d FS ( FS t l * 

DT1S2 

2 d FS^FS t 2 ^ 


C-APSF 


DOOSF 

SF d o 

(matrix) 

DOOFS 

FS d o 

(matrix) 

DOOSS 

■ 

SS d o 

(vector) 

DOOFF 

FF d o 

(vector) 

INXSF 

INDEX 

for d SF (t) 

INXFS 

INDEX 

for d FS (t) 

INXSS 

INDEX 

f° r d ss (t) 

INXFF 

INDEX 

for d FF (t) 



CAP3 (Similar is for CAP4) 


VKT 

approach speed in knots 

(vector) 

GAMA 

(vector) 


P 

p^ (vector) 


PP 

(matrix) 


DELTA 

6 . . (matrix) 


TT 

(matrix) 


IALF1 

a ^ 



IALF2 

a 2 

• values producing 

optimal N 

IALF3 

a 3 . 



CAPAC 

optimal N capacities 


CAPMIN 

the lowest N capacities 


IALFMl 

a 1 



IALFM2 

a 2 

f values producing 

the lowest 

IALFM3 

a 3 . 




‘ configurations 


N capacities 





SUBROUTINE SDOSF 
DIMENSION DELTA (5,5) 

COMMON /SSFFS/ 

1 CO AS r CQAF , CQASAF, AM AS , A MAF , AM ASAF, AMI , VF f GAMA S , GAMAF , BE TA , DELTA . 
21 AFAR, I AS AR , DOSF , DOFS , DOSS , DOFF , I NDXSF , I NDXFS , INDXSS , INDXFF, 

3T5F , TFS , TSS , TFF 
4 , OELTSS , DELTFF , DELT SF , DELTFS 
COMMON /SF/ 

2D0SF 1 , DOSF 2 , DO Sr 3 , DOSF 4 , VI SI , V251 , V3S I , V4SI , V 1 S2 , V2S2 • V3S2* 
3DV2S1 t DV3 SI ,D V 1 52 , D V3 S2 , DV 1 S3 , D V2S3 , D V4S 3 , D0D03 
CALL SD0SF3 
PT 1 2= {DOSF3-BET A )/VF 
IF (PT12) 700,700 , 100 

100 CONTINUE 

CALL SV3S2 

IF (V3S2-PT12) 200,200,105 
105 CONTINUE 
CALL SV3S1 

IF ( V3S 1 ) 115,115,110 

110 CONTINUE 

I NDXSF = 31 1 1 
DOSF=DOSF3 
GO TO 995 
115 CONTINUE 

CALL SDV3S1 

IF (DV3S1 -DELTSF ) 125,125,120 

120 CONTINUE 

INDXSF=312I 
D0SF=D0SF3 
GO TO 995 
125 CONTINUE 

CALL SDOSF1 
CALL SV1S1 

IF { VIST) 135, 135, 130 

130 CONTINUE 
I NDXSF=1 
DQSF=C. 

GO TO 995 
135 CONTINUE 

CALL SV1S2 

PT12=(D03F1-8£TA)/VF 
IF (V1S2-PTI2) 155,155,140 
140 CONTINUE 

CALL SDV1S3 

IF ( DV 1 S3 -DELTSF ) 150,150,145 
145 CONTINUE ' 

I ND XSF= 1121 
DOSF=DOSF 1 
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GO TO 995 
150 CONTINUE 
IN OX SF— 2 
OOSF=0 • 

GO TO 995 
155 CONTINUE 

IF (VI 52} 185,185,160 

160 CONTINUE 

CALL SDV1S2 

IF { DV IS2-OELTSF ) 180, 1 80 , 165 
165 CONTINUE 

CALL SDV1S3 

IF (0V1S3-0ELTSF) 175,175,170 
170 CONTINUE 

I NDXSF= 1122 
DOSF-DOSFl 
GO TO 995 
175 CONTINUE 
INDXSF=3 
OOSF=0. 

GO TO 995 
180 CONTINUE 

INDXSF=4 . . . 

DOSF=0. 

GO TO 995 
185 CONTINUE 
INDXSF=50 

OOSF=0 • ' ' 

GO TO 995 
200 CONTINUE 

IF IV3S2) 400 ,400, 203 
203 CONTINUE 

CAL! 8DV38? 

IF (DV3S2-DELTSF) 300,300,205 
205 CONTINUE 

CALL SV3S1 

IF < V3S1 ) 215,215,210 
210 CONTINUE 

IN0XSF-31 12 
D0SF=D0SF3 
GO TO 995 
215 CONTINUE 

CALL SDV3S 1 

IF (DV3S1-OELTSF) 125,125,220 
220 CONTINUE 

INDXSF=31 22 
DOSF = OOSF 3 
GO TO 995 
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300 CONTINUE 

CALL SD0SF2 
CALL SV2S2 

PTi 2=( D0SF2-BETAJ/VF 
IF (V2S2-PT12) 310,310,305 
305 CONTINUE 
INDXSF=9 
DOSF-O . 

GU TO 995 
310 CONTINUE 

IF ( V2S2 ) 315,315, 320 
315 CONTINUE 
INDXSF=10 
D0SF=0. 

GO TO 995 
320 CONTINUE 

CALL SV2S1 

IF (V2S1) 340,340,325 
325 CONTINUE 

CALL SDV2S3 

IF ( DV2S3-DELTSF > 335,335,330 
330 CONTINUE 

INOXSF=21 12 
DQSF = DOSF 2 
GO TO 995 
335 CONTINUE 
I NOXSF= 1 1 
DOSF = 0 • 

GO TO 995 
340 CONTINUE 

CALL SDV2SI 

IF ( DV2S1-DELTSF) 360,360,345 
345 CONTINUE 

CALL SOV2S3 

IF (DV2S3-OELTSF) 355,355,350 
350 CONTINUE 

INDXSF=2122 
DOSF —DQSF 2 
GO TO 995 
355 CONTINUE 
I N D X SF= 1 2 
DOSF=0 • 

GO TO 995 
360 CONTINUE 
GO TO 125 
400 CONTINUE 

CALL SV3S1 

IF ( V 3 S i } 600,600,405 
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405 CONTINUE 

CALL SD0003 

IF ( QODQ3— DELTSF > 415,415,410 
4 1 C CONTINUE 

I NDX SF= 3113 
DQSF=D0SF3 
GO TO 995 
415 CONTINUE 

CALL SD0SF2 
call 

PT12 = '<D03F2-8ETA)/VF 
IF (V2S2-PT12) 425,425,420 
420 CONTINUE 
I NOX SF= 1 7 
OOSF = 0 • 

GO TO 995 
425 CONTINUE 

IF (V2S2) 500,500,320 
500 CONTINUE 

CALL SDOSF1 
CALL SV1S1 

IF (VIS 13 540,540,505 
505 CONTINUE 

CALL SDQSF4 
CALL SV4SI 

IF (V4S1) 510,510,515 
510 CONTINUE 
INDXSF=24 
OOSF-O. 

GO TO 995 
515 CONTINUE 

IF (V4S1) 525,525,520 
520 CONTINUE 
I NOX SF = 25 
DOSF = 0 » 

GO TO 995 
525 CONTINUE 

CALL SDV4S3 

IF ( 0V4S3-DELTSF) 530,530,535 
530 CONTINUE 
I NO XSF— 26 
D0SF=0. 

GO TO 995 
535 CONTINUE 

I N0X5F=4 113 
DOSF — DOSF 4 
GO TO 995 
540 CONTINUE 
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CALL SV1S2 

P T1 2= { DOSF l —BET A ) /VF 
IF { VI S2— PT 12) 560,560,140 
560 CONTINUE 

IF (V1S2) 570,570,160 
570 CONTINUE 

INDXSF= 1 1 23 
DOSF=DOSF1 
GO TO 995 
600 CONTINUE 

CALL SDV3S1 

IF i 0V3SI-DELTSF) 610,610,605 
605 CONTINUE 

I NQX SF=31 23 
OOSF=DOSF3 
GO TO 995 
610 CONTINUE 

CALL SOOSF1 
CALL SV1S1 

IF ( V1SI ) 620,620,615 

615 CONTINUE 
I N D X SF =30 
DOSF- C • 

GO TO 995 
620 CONTINUE 

CALL SV1S2 

IF (VIS2) 625,625,640 
625 CONTINUE 

CALL SO V 1 S3 

IF ( DV 1 S3-0ELT SF ) 635,635,630 
630 CONTINUE 

l NDXSF= 1 1 23 
DOSF=DOSF 1 
GO TO 995 
635 CONTINUE 
I NOX SF =31 
DOSF=0. 

GO TO 995 
640 CONTINUE 

PT 1 2= t DOSF 1 —BET A ) / VF 
IF (V1S2-PT12) 650,650,645 
645 CONTINUE 
IN0XSF=32 
DOSF=0 • 

GO TO 995 
650 CONTINUE 

CALL SD0SF2 
CALL SV2S2 
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IF t V2S 2 ) 655,655,660 
555 CONTINUE 
I NDXSF— 33 
DOSF=0. 

GO TO 995 
660 CONTINUE 

PT12=( D0SF2-8ETA) /VF 
IF (V2S2-PT12) 670,670,665 
665 CONTINUE 
I NDXSF =34 
DOSF=0 • 

GO TO 995 
670 CONTINUE 

CALL SV2S1 

IF ( V2S I J 680,680,675 
675 CONTINUE 
INDXSF=35 
DOSF=0. 

GO TO 995 
680 CONTINUE 

CALI SDVPS1 

IF { DV2S1-DELT SF ) 685,685,690 
605 CONTINUE 
I ND X SF = 36 
DOSF=C . 

GO TO 995 
690 CONTINUE 

CALL SDV2S3 

IF ( DV2S3— DELTSF > 693,693,696 
693 CONT INUE 
IN0XSF=37 
D0SF=0 • 

GO TO 995 
696 CONTINUE 

I NDXSF=21 22 
DOSF= DQSF2 
GO TO 995 
700 CONTINUE 

CALL SV3S1 

P T 12=1 OOSF3— BETA ) /VF 
IF (V3S1-PT12) 710,710,705 
705 CONTINUE 
-V: INDXSF=30 1 1 
DOSF=DOSF3 
GO TO 995 
710 CONTINUE 

CALL S0V3S1 

IF ( DV3S 1 —DELTSF } 720,720,715 
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715 CONTINUE 

INDXSF=3021 
DOSF — 00 SF3 
GO TO 995 
720 CONTINUE 

CALL SD0SF1 
P T 1 2= <D0SF1— BETA) /VF 
IF { PT 1 2 ) 725,725,800 
725 CONTINUE 

I NDX SF= 10 21 
OOSF=OOSF1 
GO TO 995 
800 CONTINUE 

CALL SV1S1 

IF (V1S1) 810,810, 80S 
805 CONTINUE 
I NDXSF— 38 
DOSF=0-« 

GO TC 995 
810 CONTINUE 

CALL SV1S2 

IF ( V IS 2 > 815,815,320 
315 CONTINUE 

INDXSF= 1123 
DOSF^DOSFl 
GO TO 995 
820 CONTINUE 

PT 1 2= ( OOSF 1 —BET A ) /VF 
IF (V1S2-PT12) 830,830,325 
325 CONTINUE 

I NDX SF = 1121 
DOSF=DOSFl 
GO TO 995 
830 CONTINUE 

CALL SDV1S2 

IF (OV1 S2-DELTSF) 840,840,835 
835 CONTINUE 

I NDX SF = 1 122 
OOSF=OOSF1 
GO TO 995 
840 CONTINUE 

CALL SDOSF2 
CA* I 9VPS2 

IF~( V2S2 ) 845,845,850 
845 CONTINUE 
IN0XSF=39 
DOSF= 0 * 

GO TO 995 
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850 CONTINUE 

PT12= (D0SF2— BETA) /VF 
IF IV2S2-PT12) 860,850,855 
855 CONTINUE 
INQXSF=40 
DOSF-O. 

GO TO 995 
850 CONTINUE 

INDXSF=21 22 
OOSF = OGSF 2 
GO TO 995 
995 CONTINUE 
RETURN 
END 


SUBROUTINE SOOSF1 
DIMENSION DELT A < 5 , 5 ) 

COMMON /SSFFS/ 

1COAS, COAF , COASAF, AMAS, AMAF , AM ASAF, AM I , VF , GAMAS, GAMAF, BETA, DELTA , 
21 AFAR, I ASAR, DOSE, DOFS, DOSS ,DOFF, INDXSF, I NDXFS , I NQXSS , INDXFF, 

3TSF ,TFS ,TSS,TFF 
4 , DELT SS , DELTFF, DELT SF , DELTFS 
COMMON /SF/ 

2DOSF1 , D0SF2 , D0SF3, D0SF4 , VI SI , V2S I , V3S 1 , V4SI , V IS 2, V2S2, V3S2, 
3DV2SI ,DV3S1 ,DV1 S2,DV3S2 , D V 1 S3 , D V2S3 , D V4S3 , D0003 
ASF 1= 1 1 .-AM I^COASAF ) **2. /AMASAF 

8SF1=-2.*BETA*( 1 .-CQAF+ (AMI ^COASAF— 1 • )*! 1 .-COAF-AMI* 

1 < CO ASAF— COAS ) >/AMASAF) 

C 1 SF l = BET A * * 2 . 4 ( 2 ( 1 .-COAF)-( 1 • -COAF-AM I * { CO AS AF-COAS > )*#2/ 

1 AMASAF ) 

C SF 1=C 1 SF 1 — DELTSF 4^2 
DETSF l=8SFl**2-4.=!tASFl *CSF 1 
IF iDETSFl) 10,20,30 
10 CONTINUE 

PRINT lOOjIASAR, IAFAR, DETSF1 

100 FORMAT (//, ❖SUBROUTINE SDOSF 1 , DETERMI NANT NEGATIVE *,3X, 
l*ALFAS =* , I3,3X,*ALFAF =* , I 3, 3X, #DETSF1 -#,F10.S/) 

GO TO 40 
20 CONTINUE 

DOSF 1=-8SF1/(2.*ASF1 ) 

GO TO 40 
30 CONTINUE 

D01SFI=I-BSF1-S0RT(DETSF1 ) )/(2.*ASFl ) 

D02SFl= (-BSF1 +SQRT( DETSF 1 ) ) / ( 2 .*ASF 1 ) 

DOSF 1=D02SF1 
40 CONTINUE 
RETURN 
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END 


SUBROUTINE SD0SF2 
DIMENSION DELTA (5*5) 

COMMON /SSFFS/ 

ICO AS t CO AF,COASAF, AM AS , AMAF, AMASAF, AMI ,VF r GAMAS, GAM AF, BETA, DELTA, 
21 AFAR, IASAR,DQSF f DOFS#DOSS f DOFF f INDXSF, INDXFS, INDXSS, INDXFF, 

3TSF ,TFS,TS5, TFF 
4 t DELTSS fDELTFF , DELTSF , DELTFS 
COMMON /SF/ 

2D0SF1 ,OOSF2,DOSF3 , D0SF4 , VI SI , V2S1 , V3S1 , V4S1 , V 1 S2 , V2S2 , V3S2 , 
3DV2S1 ,DV3SI ,DV1S2,DV3S2 , D V 1 S3 , D V253 , D V4S 3 , D0DO3 
ASF 2=1 I AMI *CUAF}**2./ AMAF 

BSF2=-2#*BETA*( 1 ,-COAF- ( 1 .-AMI*CGAF)« ( 1 •~COAF)*{l . + AMI ) /AMAF ) 
C15F2=BETA**2,*(2.*C1 .-COAF J~( I «-CO AF I ♦ * 2 .* ( l.+AMI )**2./AMAF) 
CSF2=C 1SF2-DELTSF**2 
DETSF2=8SF2**2-4.*ASF2*CSF2 
IF (DETSF2) 10,20,30 
10 CONTINUE 

PRINT 100,1 AS AR , I AFAR , DETSF2 

ICO FORMAT (//, ^SUBROUTINE SDOSF2 , DETERM I NANT NEGATIVE *,3X, 

1*ALFAS , I3,3X,*ALFAF =*, I3,3X,*DETSF2 =*,F10.5/) 

GU TO 40 
20 CONTINUE 

DUSF2--BSF2/( 2 .4ASF2) 

GO TO 40 
30 CONTINUE 

001 SF2 = {— BSF2— SORT C DETSF2 ) ) / ( 2 • *A SF2 ) 

00 2SF2= { — BSF2 +SQRT C DETSF2 ) )/( 2 • 5 J t ASF2 ) 

D0SF2=D02SF2 
40 CONTINUE 
RETURN 
END 


SU3R0UTINE SD0SF3 
DIMENSION DELT A { 5 , 5 ) 

COMMON /SSFFS/ 

ICO AS , CO AF , CO ASAF , AM AS , AMAF, AM ASAF , AM I , VF , G AMAS , GAM AF, BETA , DELT A , 
21 AFAR, IAS AR, DOSF, DOFS, DOSS, DOFF, I NDXSF, I NOXFS , I NDXSS , I NDXFF , 

3T SF , TF S , T SS , TFF 
4 , DEL TSS , DELTFF ,DELTSF , DELTF S 
COMMON /SF/ 

2DOSF1 , DOSF2 , D05F3 , D0SF4 , Vl SI , V2S1 , V3S 1 , V4S 1 , V 1 S2, V2S2, V3S2 , 
3DV2S1 , DV3 SI , D V 1 S2 ,D V3S2 , D V l S3 , D V2S3 , D V4S3 ,D0D03 
D0SF3=DELTSF+< 1 .-AMI )^GAMAS/AMI 
RETURN 
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END 


SUBROUTINE SD0SF4 
DIMENSION DELTA<5,5) 

COMMON /SSFFS/ 

1 CO AS, COAF, CO A SAP, AM AS* AMAF, AM AS AF, AM I , VF , GAM AS , GAM AF, BETA, DELTA. 

2 I AFAR, I ASAR,DOSF, DOFS, DOSS, DOFF, INDXSF, I NDXFS , INDXSS, INDXFF, 

3TSF ,TFS,TSS, TFF 

4 , DELTSS , DELTFF , DEL T SF , DELTFS 
COMMON /SF/ 

2D0SF1 , D0SF2 , DOS F3 , DOSF4 ,V1S1,V2S1 , V3S1 ,V4S1 , V 1 S2, V2S2 , V3S2 , 

3D V2S1 ,DV3S1 ,D VI S2,DV3S2 ,DVI S3 , DV2S3 , D V4S 3 , D0D03 
ASF4-1. 

eSF4=-2.*BETA*{ 1 .-COAF) 

C 1SF4=2.*BETA**2**{ 1 .-COAF > 

CSF4=C 1 SF4-DELTSF**2 
DETSF4=BSF4**2-4.*ASF4*CSF4 
IF ( DETSF4 ) 10,20,30 

10 CONTINUE 

PRINT 1 00 , 1 ASAR, I AFAR , DETSF4 

100 FORMAT I //,* SU8R0UT I NE SD0SF4 , DETERMI NANT NEGATIVE *,3X, 

1 rALF.AS , 13 , 3X , * ALFAF =*, 13, 3X,*DETSF4 =4,F10.5/) 

GO TO 40 
20 CONTINUE 

D0SF4=-BSF4/(2.*ASF4) 

GO TO 40 
30 CONTINUE 

DQ1 SF4=<-3SF4-SQRT(0ETSF4) >/<2.*ASF4> 

002SF4= (— 8SF4+SQRT ( 0ETSF4 ) )/(2 .FASF4) 

D0SF4— D02SF4 
40 CONTINUE 
RETURN 
END 


SUBROUTINE SV1S1 
DIMENSION DELTA (5 ,5 ) 

COMMON /SSFFS/ 

1 CO AS , COAF , COAS AF , AM AS , AMAF , AM AS AF , AM I , VF , GAMAS, GAMAF , BE TA, DELTA , 
21 AFAR, I ASAR, DO SF, DOFS, DOSS, DOFF, I NDXSF , I NDXFS , INDXSS , I NDXFF, 
3TSF,TFS,TSS,TFF 
4 , DELTSS, DELTFF, DELTSF, DELTFS 
COMMON ./SF/ 

2DOSF 1 , D0SF2 , D0SF3 , D0SF4 , VI SI , V2S 1 , V3S1 , V4S1 , V 1 S2 , V2S2 , V3S2 , 
3DV2S1 ,DV3S1 , DV1S2,DV3S2 , D VI S3 , D V2S3 , D V4S3 ,00003 
V3SI=-{D0SF1*( AMIVC0ASAF-1 • )+BETA*{ 1 .-COAF-AMI* (COASAF-COAS) ) ) 
l/< AMASAFSVF) 
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RETURN 

END 


SUBROUTINE SV2S1 
DIMENSION DELTA! 5,5 ) 

COMMON /SSFFS/ 

I CO AS, COAF, CO A SAP, AM AS , A M AF , AM ASAF , AM I , VF f G AM AS, GAMAF , BETA , DELTA , 
21 AFAR , I ASAR , DO SF , DOF S , DOSS , DOFF , I NDXSF, I NDXFS , I NDXSS, I NDXFF, 
3TSF,TFS,TSS,TFF 
4 , DELTSS , DELTFF , DELTSF , DELTFS 
COMMON /SF/ 

2D0SF1 , DOSF2, DOSF3 f D0SF4, VI SI , V2S1 ,V3S1 , V4S1 , V 1 S2 , V2S2 , V3S2 t 
3DV2S1 ,D V3S1 , D VI S2 , D V3S2 , O VI S3 , 0 V2S3 , D V4S 3 f D0D03 
V2S1=-(D0SF2*( AMI*COASAF.-l • >+3ETA*< 1 . -CO AF-AM I* { COAS AF-COAS ) ) ) 

1/ ( AMASAF*VF) 

RETURN 

END 


SUBROUTINE SV3S1 
DIMENSION DELTA (5, 5) 

COMMON /SSFFS/ 

ICOAStCOAF, CO ASAF, AM AS, AMAF, AM A S AF, AM I , VF , G AM AS , GAM AF , BET A , DELTA , 
2 I AFAR, I AS AR,OOSF, DOFS, D OSS , DOFF , I NOXSF, I NDXFS , I NDXSS , I NDXFF , 
3TSF,TFS,TSS,TFF 
4 , DELTSS , DELTFF , DELTSF , DELTFS 
COMMON /SF/ 

200SF1 ,DQ3F2,D0SF3 f D0SF4 , VI SI , V2SI f V3Sl , V4S 1 , V 1S2, V2S2, V3S2 » 
3DV2S1 ,DV3SI ,DV1S2, D V3S2 ,D VI S3 ,DV2S3,D V4S 3 , DODQ3 
V3SL=-(D0SF3*{ AMI*C0A3AF-1 . )+BETA*< I , -CO AF-AM I M ( CO AS AF-CO AS > ) > 
I/( AMASAF^VF) 

RETURN 

END 


SUBROUTINE SV4SI 
DIMENSION DEL T A < 5 , 5 ) 

COMMON /SSFFS/ 

1COAS, COAF , CO ASAF, AM AS, AMAF, AMASAF, AMI , VF f GAMAS, GAMAF , BET A , DEL.TA, 
21 AFAR, I ASAR, DOSF, DOFS , DOSS , DOFF , l NDXSF, I NDXFS, I NDXSS , I NDXFF , 

3TSF , TFS,TSS,TFF . 

4, DELTSS, DELTFF , DELTSF, DELTFS 
COMMON /SF/ 

2DOSF 1 , DCSF 2 , DOSF 3 , DOSF 4 , VI SI , V2SI ,V3S1 ,V4S1 , V 1 S2 , V2S2 , V3S2 . 
3DV2SI ,DV3SI ,DV1S2,DV3S2 ,DVIS3, DV2S3,DV4S3»D0D03 
V4SI = - ( DOSF 4# ( AM I CO ASAF— 1 . J+BETA* ( 1 • — CO AF— AM I* ( COAS AF— CO AS ) } ) 
1/1 AMASAF# VF > 
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RETURN 

END 


SUBROUTINE SV1S2 
DIMENSION DELTA (5,5) 

COMMON /SSFFS/ 

1 CO AS t CO AF , CCA SAF T AM AS , AM AF , AMASAF ,AMI , VF , GAM AS , GAMAF, BET A , DELT A, 
21 AFAR, I ASAR,OOSF, DOFS, DOSS, DOFF, I NDXSF, I NDXFS , I NDXSS, INDXFF, 

3TSF , TFS , TSS , TFF 
4, DELTSS .DELTFF,DELTSF,DELTFS 
COMMON /SF/ 

2DOSF1 , DO SF2 , DO SF3 , DO SF 4 , VI SI , V2SI , V3S1 , V4 S 1 , VI S2, V2S2 , V3S2 , 
3DV2S1 ,D V3S1 ,D VI S2 , DV3S2 , DV 1 S3 , D V2S3 , D V4S 3 , DO»03 
V3S2={ DOS FI* ( I .-AMI*COAF )-BETA*< 1 .-COAr}*( I ,+AMl ) ) /( AMAF*VF) 
RETURN 
END 


SUBROUTINE SV2S2 
DIMENSION DELT A { 5 , 5 ) 

COMMON /SSFFS/ 

1 CO AS, CO AF , CO A SAF , AM AS , AM AF , AMASAF, AMI , VF, GAMAS, GAMAF ,BETA , DELTA , 

2 I AFAR, IAS AR, DOSF, OOFS, DOSS, DOFF, I NDXSF, I NDXFS , I NDXSS ,1 NDXFF , 

3TSF, TFS, TSS, TFF 

4 , DEL TSS , DEL TFF , DELTSF , DELTF S 
COMMON /SF/ 

2DQSF 1 ,00SF2,D0SF3,D0SF4 , VI SI , V2S1 , V3S1 ,V4S1 , V 1 S2 , V2S2 , V3S2 , 

3DV2S1 , D V3 SI , D V 1 S2 , D V3S2 , D V 1 S3 , D V2S3 , D V4S3 , D0D03 
V2S2=<OOSF2*( 1 •— AM I *CO AF ) — BET A* { l.-COAF)*( I , + AMI > ) / { AM AF* VF ) 
RETURN 
END 


SUBROUTINE SV3S2 
D IMENSION DELT A( 5, 5 ) 

COMMON /SSFFS/ 

1 CO AS, CO AF, CO A SAF, AM AS, A MAF , AMASAF, AMI, VF, GAMAS, GAM AF , BET A , DELTA , 
2 I AFAR , I AS A R , DOSF , OOFS , DOSS , DOFF , I NDXSF , I NDXFS , I NDXSS , I NDXFF , 

3TSF ,TFS , TSS, TFF 
4 , DELT S S , D EL TFF , D ELT SF , D EL TF S 
COMMON /SF/ 

2D0SF1 , D0SF2 , DO SF3 ,DQSF4,VIS1 , V2SI , V3S 1 , V4 SI , VI S2 , V2S2 , V3S2 , 
3DV2S1 , DV3S1,D VI S2,DV3S2,DV1 S3, DV2S3 ,D VAS3,D0D03 
V3S2=(D0SF3*( 1 .-AMI*CQAF)-BETA*n .-COAF)*( l.+AMI) )/{AMAF*VF) 
RETURN 
END 
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SUBROUTINE SDV2S1 
DIMENSION DELTA (5 ,5) 

COMMON /SSFFS/ 

ICO AS, COAF, COASAF, AM AS, AMAF, AMASAF, AMI ,VF, GAMAS, GAMAF, BETA , DELTA* 
21 AFAR, I ASAR, DOSF, DOFS, DOSS, DOFF, I NDXSF , I NDXF5 , INDXSS, INDXFF, 

3TSF ,TF S, TSS, TFF 
4 , DELTSS , DELTFF , DELTSF , DELTFS 
COMMON /SF/ 

2DOSFI , D0SF2 ,DGSF3,DOSF4 , VI SI , V2S1 , V3S1 , V4S 1 , V 1 S2, V2S2, V3S2 , 
3DV2S1 ,DV3S1 ,0V1S2,DV3S2 ,D V 1 S3 , D V2 S3 , 0 V4S 3 , D0D03 
ASF 1=1 .-{1 .-AMI*C0ASAF)*F2./AMASAF 

RSF1=-2.*8ETAF ( 1 .-COAF+ ( AM I FCOASAF- 1 • ) *{ 1 ,-COAF-AMI* 

1 (COASAF-COAS) ) /AM ASAF ) 

CISFI=BETA*F2.F(2 .*< 1 *-COAF )-< 1 .-COAF-AM lF<COASAF-COAS)J **2/ 

1 AM ASAF ) 

A=DCSF2*£2 .FASFl 4-D0SF2’i'eSFl-«-C l SF I 
IF (A) 10,20,20 

10 CONTINUE 

PRINT 100, I ASAR, IAFAR, A 

100 FORMAT {//, ’^SUBROUTINE SDV2S1 ,SQ* RUOTE ARGUMENT NEGATIVE F,3X, 
l*ALFAS = *, 13, 3X, *ALFAF =* , I 3, 3X , *SQRT • ARG. =*,F10.5/} 

GO TO 40 
20 CONTINUE 

DV2S1— SORT (A) 

40 CONTINUE 
RETURN 
END 


SUBROUTINE SDV3S1 
DIMENSION DELTA (5, 5) 

COMMON /SSFFS/ 

1 CO AS , CO AF , COASAF, AM AS, AMAF, AM AS AF, AMI , VF , GAMAS, GAMAF, BETA, DELTA, 
21 AFAR, I ASAR, DOSF y DOFS, DOSS , DOFF, l NDXSF , I NDXFS , INDXSS , I ND XFF , 

3TSF , TFS, TSS, TFF 

4, DELTSS, DELTFF, DELTSF, DELTFS 

COMMON /SF/ : — — 

2DOSF1 , DQSF2 , DOSF 3 ,D0SF4,V1S1 ,V2S1 ,V3S1 , V4SI , V 1 S2 , V2S2 , V3S2 , 
3DV2S1 ,DV3S1 , D V 1 S2 , 0 V3S2 ,D V 1 S3 , D V2S3 , D V4S3,00D03 
ASF 1=1 .-( 1 .-AMIFCOASAF) FF2./AMASAF 

3SFl=-2.*BETA*( 1 •-COAF + ( AM I *COA SAF- 1 . ) * { l.-COAF-AMI* 

1 (COASAF-COAS) ) /AM ASAF ) 

C 1 SF 1 = BET A**2 .*{ 2.4(1 .-COAF >- ( 1 .-COAF-AM IF I COA S AF— CO AS ) >**2/ 
l AM ASAF ) 

A=DOSF3^F2.FASFI+DOSF3FBSF1+C1SF1 
IF I A) 10,20,20 
10 CONTINUE 
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PRINT 100, IASAR, IAFAR, A 

100 FORMAT { //,*SUSROUTlNE SDV3S1 ,SQ. ROUTE ARGUMENT NEGATIVE *,3X, 
l^ALFAS — * t I3f SXj^ALFAF 13, 3X,TS0RT.ARG. =*,F10.5/1 

GO TO 4 0 __ 

20 CONTINUE 

DV3S1=SQRT( A) 

40 CONTINUE 
RETURN 
END 


SUBROUTINE SOV1S2 
DIMENSION DEL TA ( 5 , 5 ) 

COMMON /SSFFS/ 

1 CO AS , CO AF , CCA SAP, AM AS, A MAP , AM AS A F, AMI , VF , GAMAS , GAMAF, BET A , DEI.TA , 
21 AFAR, IASAR, DOSF, DOFS, DOSS, DOFF, INDXSF, I NDXFS , I NDXS S , INDXFF, 

3TSF , TFS, TSS,TFF 

4 , DEL T SS , O EL TFF , D EL T SF , DEL TFS 
COMMON /SF/ 

2D0SF 1 , DOSF 2 , DOSF 3 , DOSF 4 , V1S1 ,V2S1 , V3S1 ,V4S1 , VI S2 , V2S2 , V3S2 , 
3DV2S1 ,DV3SI , 0 V 1 S2 , D V3S2 , D V 1 S3 , DV2S3 , D V4S3 , DODO 3 
A SF 2= 1 .-{ 1 • —AM I FCOAF ) **2./AMAF 

8SF2=-2.*BETA*( 1 .-COAF-( 1 . - AM I «C0 AF ) * I l.-COAF)*!! ,+AMI I/AMAF) 

Cl SF2=3ET AXc'X2 . # ( 2 * * ( 1 . -COAF )- ( 1 ,-COAF ) =**2 .*{ 1 • fA MI) **2./AMAF } 
A=DQSF1**2.*ASF2+D0SF 1 4BSF2+C1SF1 
IF ( A ) 10,20,20 

10 CONTINUE 

PRINT i 00 , IASAR , IAFAR , A 

100 FORMAT (//,'XSUDROUT INE SDV1S2,SQ. ROUTE ARGUMENT NEGATIVE *,3X, 

I X'ALFAS , 13 , 3X, *ALFAF , 1 3 , 3X , #SQRT . ARG . =4, FI 0.5/) 

GO TO 40 
20 CONTINUE 

DV1S2=S0RT(A) 

40 CONTINUE 
RETURN 
END 


SUBROUTINE SDV3S2 
DIMENSION DEL T A { 5 * 5 ) 

COMMON /SSFFS/ 

ICO AS, COAF, COASAF, AM AS, AM AF,AMASAF, AMI, VF , GAMAS, GAMAF, BET A, DELTA, 
2 I AFAR, IASAR, DOSF , DOFS , DOSS , DOFF , I NOXSF, I NDXFS , INDXSS, INDXFF, 

3TSF , TFS, T5S,TFF 
4 , DELTSS , DELTFF , DELTSF , DELTFS 
COMMON /SF/ 

2D0SF 1 , DOSF 2 , D0SF3 , DOSF 4 , VIS! , V2S1 , V3S1 ,V4S1 , V 1 S2 , V2S2 , V3S2 , 
3DV2S1,DV3SI ,DVIS2,DV3S2 ,DV1 S3 , DV2S3 , DV4S3 ,D0D03 
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ASF 2= 1 •— ( 1 .-AMt*C0AF)**2./AMAF 

BSF2=-2.*BETA*{ 1 .-COAF-I 1 .-AMI4COAF)* C l.-COAF>*M .+AMI J/AMAF) 
C1SF2=BETA**2.*C 2*^(1 .-COAF) - ( 1 .-COAF ) ** 2.* ( 1,+AMI ) 442 . /A.MAF ) 
A=DOSF3442 .4ASF2+D0SF348SF2+C ISF2 
IF CA) 10,20,20 
10 CONTINUE 

PRINT TOO t T AS AR , I AFAR , A 

100 FORMAT {//^SUBROUTINE SDV3S2,SQ. ROUTE ARGUMENT NEGATIVE 4,3X, 
14ALFAS = 4, I3,3X f 4ALFAF =4, l 3 , 3X , 4SQRT . AR G . =4,F10.5/) 

GO TO 40 
20 CONTINUE 

OV3S2=SQRT( A) 

4C CONTINUE 
RETURN 
END 


SUBROUTINE SOV1S3 
DIMENSION DELTA (5, 5) 

COMMON /SSFFS/ 

ICO AS, COAF , CO A SAP, AM AS , AMAF, AMASAF, AMI ,VF , GAM AS, GAMAF, BETA , DELTA . 
21 AFAR, IASAR,DOSF, DOFS, DOSS, DOFF, INDXSF, I NDXFS , I NDXSS , I NDXFF , 
3TSF, TFS, TSS, TFF 
4 , DELT SS , D EL TFF , DELTSF , D EL TFS 
COMMON /SF/ 

2D0SF1 , D0SF2 , D0SF3 , DO SF4 , VI S 1 , V2S l ,V351 f V4Sl , V 1 S2, V2S2 , V3S2 , 
3DV2S1 -DV3S1 ,DV1S2 ,DV3S2 ,DV1S3, DV2S3 , D V4S 3 , DODO 3 
DV1S3=BET AP< DCSF1 -BETA ) 4AMI-VF4{ 1 .-AMI ) 4( GAM AS/ C AM 14 VF ) — 

1 ( DOSF1-BET A >/VF) 

RETURN 

END 


SUBROUTINE SDV2S3 
DIMENSION OELTACS, 5) 

COMMON /SSFFS/ 

1 COAS , CD AF , CO AS AF , AM AS , AMAF , AMASAF , AM I , VF , GAMAS, GAMAF, BETA , DELTA , 

2 I AFAR , I AS AR, DOSF , DOFS, DOSS ,D0FF, INDXSF, I NDXFS , I NDXSS , INOXFF. 

3TSF, TFS, TSS, TFF 

4, DEL TSS , DELTFF , DELTSF, DELTFS 
COMMON / SF/ 

2DOSFI ,DasF2,D0SF3,D0SF4, V1S1 , V2S 1 , V3S I , V4S 1 , V I S2 , V2S2, V 3S2 , 

3D V2S1 , D V 3 S 1 , D V 1 S2 , D V3S2 , D VI S3 , O V2 S3 , 0 V4S3 , D0DO3 
DV2S3=BETA+-( D05F2-QET A ) *AM I -VF4 < 1 .-’AM I ) *( GAMAS/ (AMI 4 VF ) — 

1 I DO SF 2-3ET A ) / VF ) 

RETURN 

END 
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SUBROUTINE SDV4S3 
DIMENSION DELTA(5,5) 

COMMON /SSFFS/ 

1 CO AS, COAF, CO ASAF, AM AS, AMAF, AM ASAF, AM I , VF , G AM AS , GAM AF, BETA , DELTA , 
21 AFAR , IAS AR, DOSF, DOFS, DOSS, DOFF, I NDXSF, I NDXFS , I NDXS5 , INDXFF , 
3TSF, TF S , TSS , TFF 
4 , DELTSS , D ELTFF , DELTSF , DELTFS 
COMMON /SF/ 

2DOSF1 , DOS F 2 , D0SF3 , DOSF4 , VI S 1 , V23 1 , V3S 1 , V4S 1 , V 1 S2, V2S2, V3S2 , 

30 V2S1 t DV3Sl ,DVIS2,0V3S2 ,DV1 S3 , D V2S3 ,0 V4S3 , 00003 
DV4S3=8ETA-H D0SF4— 3 ET A ) *AMI-VF*{ 1 .-AMI ) * { CAMAS/ (AMI ) — 

1 (D0SF4-BETA)/VF) 

RETURN 

END 


SUBROUTINE SD0D03 
DIMENSION DELTA (5, 5) 

COMMON /SSFFS/ 

1C0AS, COAF , CO ASAF, A MAS, AMAF , AMASAF .AMI , VF , GAMAS , GAMAF, BET A , DELTA , 
21 AFAR , I ASAR, DO SF , DOFS , DO SS , DOFF , I NDXSF, I NDXFS , INDXSS , INDXFF, 

3T SF , TF S , TSS , TFF 

4, DELTSS, DELTFF, DELTSF, DELTFS 
COMMON /SF/ 

2D0SF1 , D0SF2 , OOSF3 , DOSF 4 , Vi SI ,V2SI ,V3S1 ,V4S1 , V 1 S2 , V2S2 * V3S2 , 
3DV2S1 ,DV3SI , D VI S2 , D V3 S2 , D V 1 S3 , D V2S3 , D V4S3 , D0003 
D0DQ3-=D0SF3**2.-2.*D0SF3*BETA*( 1 .-COAF ) +2 .*8ETA**2 •* ( l • “COAF ) 
RETURN 
END 


SUBROUTINE 3D0FS 
DIMENSION DELT A{5,5 ) 

COMMON /SSFFS/ 

I CO AS , COAF, COASAF, AM AS , AMAF, AMASAF, AMI ,VF , GAMAS , GAMAF, BETA , DELTA , 
21 AFAR, I ASAR, DOSF, DOFS, DOSS, DOFF, INDXSF, INDXFS, I NDXSS , I NDXFF , 

3TSF, TFS,TSS t TFF 
4, DELTSS, DELTFF, DELTSF, DELTFS 
COMMON /F 5/ 

1 DOFS 1 , D0FS2 , DOF S4 , D0FS5 , T 1 S 1 , T2S 1 ,T4SI ,T 1 S2 » T2S2 , T4S2 , DT2 SI , DT 1 S2 

2 , T 5 S 1 , D T 5 S I 
T!2=“B£TA/VF 
CALL SDOFSI 
CALL ST1S1 

IF I T 1 S I “ T I 2 ) 100,100,120 

100 continue 
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CALL ST1S2 

IF (T1S2-T12) 105,105.110 

105 CONTINUE 

OOFS=OOFS1 
INDXFS-10 10 
GC TO 995 
110 CONTINUE 

CALL SDT1S2 

IF ( DTI S2-0ELTFS) 1 20 ,120, 1 15 
115 CONTINUE 

DGFS=D0FS1 
I NDXF S- 1020 
GO TO 995 
120 CONTINUE 

CALL SQOFS2 
CALL 5T2S2 

IF IT2S2-T12) 300,300,125 
125 CONTINUE 

IF ( T2S2) 130,130,200 

130 CONTINUE 

CALL ST2S1 

IF ( T2S1-T12) 1 AO, 140, 135 

135 CONTINUE 

DOF 5= DOFS 2 
I NDXF S-20 1 0 
GO TO 995 
140 CONTINUE 

CALL SDT2S1 

IF { DT2S1-DELTFS) 150,150,145 
145 CONTINUE 

D0FS=D0FS2 
I NDXF S=2020 
GO TO 995 
150 CONTINUE 
DOF S= 0 • 

I NDXF S= 1 
GO TO 995 
200 CONTINUE 

T5=G AMAS/VF 

IF (T252-T5) 205,205,240 
205 CONTINUE 

CALL ST2S1 

IF ( T 25 1 — T 1 2 ) 215 , 215,210 
210 CONTINUE 

DOFS-DOFS2 
I NDXF 5- 21 10 
GO TO 995 
215 CONTINUE 
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CALL SDT2SI 

IF { DT2S1— DELTFS ) 225,225,220 

220 CONTINUE 

DOFS=DGFS2 
INDXF 5=2120 
GO TO 995 
225 CONTINUE 
OOF S= C . 

I NDXFS=2 
GO TO 995 
240 CONTINUE 

CALL SO OF S 5 , 

CALL ST5S1 

IF (T5S1-T12) 250,250,245 
245 CONTINUE 

OOF S= DOFS 5 
INDXFS=51 10 
GO TO 995 
250 CONTINUE 

CALL SDT5S1 

IF { DT5S1-DELTFS) 260,260,255 
255 CONTINUE 

DQFS=D0FS5 
I NDXFS=5120 
GO TO 995 
260 CONTINUE 

CALL SDOFS1 
DOFS=DOFSl 
T NDXFS= 1 030 
GO TO 995 
3C0 CONTINUE 

CALL SD0FS4 
CALL ST4S1 
CALL ST4S2 

IF (T4S1-T12) 320,320,305 
305 CONTINUE 

IF (T4S2-T12) 310,310,315 
310 CONTINUE 

DGFS=D0FS4 
INDXFS=4000 
GO TO 995 
315 CONTINUE 
DOF 5=0 • 

IN0XFS=3 
GO TO 995 
320 CONTINUE 

IF (T4S2-T12) 325,325,330 
325 CONTINUE 


174 



OOF S= 0 • 
INDXFS=4 
00 TO 995 
330 CONTINUE 
DOFS = 0 • 

I NOXFS— 5 
GO TO 995 
995 CONTINUE 
RETURN 
END 


SUBROUTINE SD0FS1 
DIMENSION DELTA<5,5) 

COMMON /SSFFS/ 

1CCAS,C0AF ,COASAF,AMAS, AM AF,AMASAF, AMI ,VF, GAM AS, GAM AF, BETA, DELTA, 

2 1 AFAR , I AS AR , DOSF , DOFS , OQSS , OOFF , I NOXSF , I NDXFS , t ND XSS, I ND XFF , 

3TSF , TF S,TSS,TF r 
4 , DELTSS , DELTFF , DELT SF , DELTFS 
COMMON /FS/ 

1 DOF SI , DOFS 2 , DQF S4 , DOF S5 , T 1 S 1 ,T2S1 ,T4S1 . T I S2 , T 2S 2» T4S2 , DT2S1 , OT 1 S2 
2, T5S1 , DT55 1 

AFS1 =1 CQASAF-AMT )+*2/AMASAF 

8^S1 = 2 , *BETA*< ( CO A SAF— C 0 AS ) — I CO AS AF— AM I } $ { { 1 ♦ -CO AF ) - AM I * 

1 ( CO A S AF— CO AS) t/AMA^AF ) 

C1FS1 = BETA**2*C2.*( 1 .-COAF)-( ( 1 .-CO AF ) -A M I * ( COASAF-CO AS ) ) * *2/ 

I AMASAF) 

CFSl=CiFSl-DELTFS**5> 

DETFS 1 = BF S 1 ^ *2—4 • ^AF S 1 ^CFS I 
IF ( DETFS 1 ) 10,20,30 

10 CONTINUE 

PRINT 100,1 AS AR , I AFAR , DETFS 1 

100 FORMAT { // 1 $ SU3R0UT I NE 5D0FS1 , DETERMINANT NEGATIVE *, 3X, *ALFAS =* 
1 , 13, 3X, *AU=AF =*, I3,3X,*DETFS1 =*,F10.5/) 

GO TO 40 
20 CONTINUE 

D0FSl=-8FSl/( 2.*AFS1 ) 

GO TO 40 
30 CONTINUE 

DO IPS i~(-SFSl -SORT ( DETFS 1 > >/<2.*AFSl) 

D02FS 1 = ( — 3F S 1 +SQR T < DETF Si ) >/< 2.*AFS1 ) 

OOFS 1— D02F S 1 
40 CONTINUE 
RETURN 
END 


SUBROUTINE SDQFS2 
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DIMENSION DELTA (5,5) 

COMMON /SSFFS/ - 

I CO AS, CO AF, COASAF, AM AS , AM AF , AM ASA F, AM I , VF , GAM AS, GAMAF , BE TA ,QELTA , 
21 AFAR, I ASAR, DOSF, DOFS , DOSS , DOFF, INDXSF, INDXFS, I NOXSS, INDXFF, 
3T5F,TFS,TSS,TFF 
4,DELTSS,DELTFF, DELTSF, DELTFS 
COMMON /FS/ 

1DOFSI ,DOFS2,DOFS4,DOFS5,T1SI ,T2S1 , T4S 1 , T 1 S2 , T2S2 , T4S2 , D T2S1 , DTI S2 
2, T5SI , D T5 S 1 
A = 1 •— ( COAS— AMI ) s r s r2/AM AS 
IF (A) 10,10,20 

10 CONTI NUE - - 

PRINT 100 , I ASAR, I AFAR , A 

ICO FORM AT ( //,*SU3R0UTI NE SDOFSL jQ, ROUTE ARGUMENT NEGATIVE * # 3X, 
l^ALFAS =*, I3,3X,*ALFAF =* , I 3, 3X , *SQRT • ARG. =4 , F 1 0 • 5/ 1 
GO TO 40 

20 CONTINUE 

OOFS2=DELTF5/SQRT{ A) 

40 CONTINUE : 

RETURN 

END 


SUBROUTINE SD0FS4 
DIMENSION DEL T A ( 5 , 5 ) 

COMMON /SSFFS/ 

1 COAS, CO AF, COASAF, AM AS, AM AF, AM AS AF, AMI, VF, GAM AS, GAMAF, BETA , DELTA, 

21 AFAR, I ASAR, DOSF, DOFS, DOSS, OOFF, INDXSF, INDXFS, INDXSS, INDXFF, 
3TSF,TFS, TSS, TFF 
4, DELTSS,DELTFF, DELTSF, DELTFS 
COMMON /FS/ 

1 DOF S 1 ,D0FS2,DQFS4,D0FS5,T1SI ,T2S1 ,T4S1 ,T IS 2, T2S2 , T4S2 , DT2S l , D T 1 S2 
2 , T5S1 ,D T5S1 
AFS4 = 1 « 

BFS4=-2.*3ETA*(C.0AS-AMI ) 

C1FS4=BETA**2*AMAS 
CFS4=C1FS4-DELTFS**2 
DETFS4=8FS4**2-4. *AFS4*CFS4 
IF ( DETFS4 ) 10,20,30 
10 CONTINUE 

PRINT 100, IASAR, IAFAR,DETFS4 

100 FORMAT (//, MSUBRQUT INE SDOFS4, DETERMINANT NEGATIVE *,3X,*ALFAS s* # 
1 I 3,3X,*ALFAF =* , I 3 , 3X , *DE TFS4 ~*,F10.5/) 

GO TO 40 ’ 

20 CONTINUE 

D0FS4-— BF S4/( 2 • ♦AF S4 ) 

GO TO 40 
30 CONTINUE 
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D01FS4=(-8FS4-SQRT(DETFS4) }/( 2.*AFS4) 
D02FS4— (-3FS4+SQRT(DETFS4) }/( 2.4AFS4J 
OOFS4=D02FS4 
40 CONTINUE 
RETURN 
END 


SUBROUTINE S00FS5 
DIMENSION DELTA<5,5> 

COMMON /SSFFS/ 

1COAS , COAF ,COA$AF, AM AS, AMAF, AMASAF, AMI ,VF , GAMAS, GAM AF, BETA, DELTA, 

21 AFAR , IASAR, DO SF t DOFS r OOSS, DOFF, INDXSF, I NDXFS f I NO XSS , I NDXFF , 

3T5F , TFS , TSSf TFF 
4, DELTSS,DELTFF,DELTSF,DELTFS 
COMMON /FS/ _ 

1 OOFS l ,DQFS2,OOFS4,OOFS5,T1 SI , T2Sl,T4Sl ,T I S2 ,T 2S2 , T4S2, DT2S1 , DT 1 S2 
2 , T5 S l , D T 5 S I 
AFS5=I * 

BFS5=GAMAS4{ COAS-AMI 1 * 2 * 

C 1 FS 5= GAM AS# *2 4 AM AS 
CFS5 = C IPS 5-DEL TFS*«2 
DETFSS=BFS5**2-4.*CFS5 
IF ( DETFS5 ) 10,20,30 

10 CONTINUE 

PRINT 1 00 , IASAR, IAFAR , DETFS5 

100 FORMAT (//.^SUBROUTINE SD0FS5 , DETERMI NANT NEGATIVE *, ‘3X, *ALFAS « 4 
1 , I3,3X,*ALFAF =*, I3,3X,*DETFS5 =*,F10.5/) 

GO TO 40 
20 CONTINUE 

D0FS5=— BFS5/2 . 

GO TO 40 
30 CONTINUE 

D0FS5= (— BFS5+SQRT (DETF55) ) / 2 • 

40 CONTINUE 
RETURN 
END 


SUBROUTINE STISI 
DIMENSION DELT A ( 5 , 5 ) 

COMMON /SSFFS/ 

1 COAS, COAF ,COASAF f AM AS, AMAF, AMASAF, AMI ,VF, GAMAS, GAMAF, BET A, DELT A, 
21 AFAR, IASAR, DO SF , DOFS , DOSS , DOFF , I ND XSF , I NDXFS , I NDXSS , I NDXFF , 

3TSF, TFS,TSS, TFF 
4,DELTSS, DELTF F , DELTSF , DELTF S 
COMMON /FS/ 

1 DOF SI , DOF S2 , D0FS4 , D0FS5 , T I S V , T2S l , T4S 1 , T IS 2 , T2S2 ,T4S2 , DT2S 1 , DT l S2 
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2, T5S1 fOTSSl 

T I SI ~ — i DOJFS 1*1 COASAF— AM I ) -*-3ETAF( I 1 ,-COAF COASAF-COAS) ) )/ 

1 < AM A S AF£ VF ) 

RETURN 

END 


SU3ROUT I NE ST2S1 
DIMENSION DELTA! 5,5? 

COMMON /SSFFS/ ' 

1 CO AS , CO AF , COASAF , AM AS , AMAF, AM AS AF, AM I , VF , G AM AS , GAM AF , SET A , DELT A , 

2 1 AFAR, IASAR, DO SF, DOFS , D OSS , DOFF , INDXSF, I NDXFS , I NDXSS , I NDXFF , 

3TSF ,TFS , TSS, TFF 

4 ,DELTSS,DELTFF, DELTSF, DELTFS 
COMMON / FS / 

1 DOFS 1 , DOF S2 , DOF S4 , DOF S5 , T I S 1 ,T2S 1 , T4S 1 , T 1S2 , T2S2 , T4S2, DT2S1 , DT I S2 
2 , T5 S l ,DT5S 1 

T2Sl=-( DOFS 2^' ( COASAF— AM I ) + 3£T A# ( ( 1 .-CO AF ) - AMI *( COASAF-COAS) > )/ 

X < AM A SAF# VF ) 

RETURN 

END 


SUBROUTINE ST4S1 
DIMENSION DELT A{ 5 , 5 ) 

COMMON /SSFFS/ 

1 COAS, CO AF, COASAF, AM AS, AMAF, AM A SAP, AMI , VF,GAMAS, GAMAF , BETA t DELT A , 
21 AFAR, IASAR* OOSF , DOFS , DOSS, DOFF , INDXSF, I NDXFS , I NDXSS, I NDXFF, 

3T SF , TF S , TSS , TFF , 

4, DELT SS, DEL TFF, DELTSF , DELTFS 
COMMON /FS/ 

ID OF SI ,DOFS2, DOFS4 ,D0FS5 ,Tl SI ,T2S1 , T4S V,T 1 S2 , T2S2 , T4S2 , D T2S 1 ,DTIS2 
2, T5S1,DT5S1 

T 4S I — — ( DQFS4# ( COASAF— AM I ) +3ET A#( ( 1 .-CO AF ) - AM I * ( CO ASAF-CO AS ) ) ) / 
i i AMASAF*VF ) 

RETURN 

END 


SUBROUTINE ST5S1 
DIMENSION DELT A ( 5 , 5 ) 

COMMON /SSFFS/ 

1COAS, CO AF, COASAF, AM AS, AMAF, AM AS AF, AMI ,VF, GAMAS, GAMAF, BETA , DELTA , 
21 AFAR, IASAR, DOSF , DOFS ,DOSS ,OOFF, INDXSF, I NDXFS , I NDXSS , I NDXFF , 

3TSF ,TFS,T3S* TFF 
4, DELTSS,DELTFF, DELTSF, DELTFS 
COMMON /FS/ 

I DOFS 1 , DOFS2,DOFS4,DOFS5 f TlSl ,T2S1 , T4S 1 , T 1S2 ,T2S2 , T4 S2 , D T2S1 , DTI S2 
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2 f T5S 1 ,DT5S1 

T5S1=- C D0FS5* (COASAF- AM I ) +BETA4C < 1 ,-CO AF } -AM I 4( CQASAF-COAS) } )/ 
1J AMASAF^VF) 

RETURN 

END 


SUBROUTINE ST i S2 
DIMENSION DELTA! 5,5) 

COMMON /SSFFS/ 

I CO AS, CCAF,COASAF, AMAS,AMAF, AMASAF, AMI T VF f GAM AS, GAMAF, BETA, DELTA, 

2 I AFAR , I AS AR, DOSF, DOFS, D OSS , OOrF , I NDXSF, I NDXFS, I NDXSS • I NDXFF , 
3TSF,TFS, TSS,TFF 
4, DELTSS , DELTFF , DELTSF , DELTFS 
COMMON /FS/ 

1 DOF SI, OOFS 2 t DOF S4 , OOF S5 , T 1 S 1 ,T2SI ,T4S1 , T I S2 , T2S2 , T4S2 , D T2S 1 , DT 1 S2 
2, T5S1 , D T5 S 1 

T 1 S2=-DOF S 1* { CO AS-AMI 3 / { AMAS*VF ) 

RETURN 

END 


SUBROUTINE ST2S2 
DIMENSION DELT A ( 5 , 5 ) 

COMMON /SSFFS/ 

1COAS , COAF, COASAF, AM AS, AMAF, AM AS AF, AM I , VF , GAMAS, GAMAF » BET A, DELTA, 

2 I AFAR, I AS AR, DOSF, DOFS, DOSS , DOFF, I NDXSF , I ND XF S , I NDXSS , I NDXFF , 
3TSF,TFS, TSS,TFF 
4, DELTSS, DELTFF , DELTSF, DELTFS 
COMMON /FS/ 

I DOFS 1 , DOF S2 , DCFS 4 , DOFS 5 , T I S 1 ,T2S1 , T4S1 , T l S 2 , T 2S2 , T4S2, DT2S1 , DTI S2 
2, T5S1 ,DT5S 1 

T2S2=-DOFS2F( COAS-AMI )/{ AMAS*VF) 

RETURN 

END 


SUBROUTINE STAS2 
DIMENSION DEL TA ( 5 , 5 j 
COMMON /SSFFS/ 

I CO AS, COAF , COASAF, AM AS, AMAF, AM AS AF, AMI , VF , GAM AS , G AMAF, BETA , DELTA , 
21 AFAR, IASAR , DOSF, DOFS, DOSS , DOFF, INDXSF, I NDXFS , I NDXSS , INDXFF , 

3TSF , TF S , TSS , TFF 

4, DELTSS, DELTFF, DELTSF, DELTFS 
COMMON /FS/ 

1 DOFS l , DOFS 2 , D0FS4 , DOFS 5 , T I S 1 , T 2S 1 , T4S 1 , T 1S2, T2S2 , T4S2, DT2S 1 , DTI S2 

2 , T 5 S 1 ,DT5SI 

T4S2=-D0FS4F( CO AS-AMI ) / ( AMAS*VF ) 
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RETURN 

END 


SUBROUTINE SDT2S1 
DIMENSION DELTA (5*51 
COMMON /3SFFS/ 

lCOAS,COAF*COASAF, AM AS, A MAF, AMASAF, AMI ,VF, GAMAS, GAM AF, BET A, DELTA, 

2 I AFAR, I ASAR, DOSF, DOFS , DOSS , DOFF , INDXSF, I NDXFS , I NDXSS , I NDXFF, 
3TSF,TFS,TSS t TFF 
4, DELTSS,DELTFF,DELTSF,DELTFS 
COMMON /FS/ 

I DOFS 1 , DCFS2 , D0FS4 , DOF S5 , T I SI ,T2SI ,T4SI , T 1 S2 ,T2S2 , T4S2, OT2S1 , DTI S2 
2 , T5 SI , D T5S i 

AFS1=I • - < COASAF-AM2 ) **2/AM ASAF 

BFS1=2 **BETA*< < CO ASAF-C O AS ) - { CO AS AF-AM I )*{ ( I,-COAF)-AMI*(COASAF- 
1COAS) ) /AMASAF ) 

ClFSl = 3ETA»*2*{2.*< l«-COAF)-< ( 1 .-COAF )-AMI* ( COASAF-COAS ) ) **2/ 

1 AMASAF) 

A-=D0FS2#*2*AFSH-D0FS2*BFS1+C1FS1 
IF (A) 10,20,20 

10 CONTINUE 

PRINT 1 00* I ASAR, IAFAR, A 

100 FORMAT < / /, ^ SU3R0UT I NE SDT2S1,SQ. ROUTE ARGUMENT NEGATIVE *,3X, 
l^ALFAS , 13, 3X, *ALFAF =* , I 3 , 3X , *SQRT • AR G . =*,F10.5/> 

GO TO 40 
20 CONTINUE 

DT2S1=SQRT( A) 

40 CONTINUE 
RETURN 
END 


SUBROUTINE SDT5S1 
DIMENSION DEL TAI 5 , 5 ) 

COMMON /SSFFS/ 

ICO AS, COAF, COASAF, AM AS , AM AF , AM AS AF , AM I , VF , G AM AS , GAM AF , BETA , DELTA , 
21 AFAR, I ASAR, DOSF, DOFS* DOSS, DOFF, I NDXSF , I NDXFS , I NDXSS , I NDXFF, 
3TSF,TFS,TSS,TFF 
4,DELTSS , DELTFF , DELTSF, D EL TFS 
COMMON /FS/ 

1 DOFS 1 , DOFS2 , D0FS4 , OOFS 5 , T1S1 ,T2S1 , T4S I * T 1 S2 , T2S2 , T4S2 , DT2 SI , DTI S2 
2 , T5S 1 , DT5 S 1 

AF31=1 .-{ COASAF-AMI )**2/AMASAF 

■ B F S 1 = 2 • * BET A * ( (COASAF-COAS ) - ( CO AS AF-ANS I )*< t 1 COAF ) - AMI * { COASAF- 
1COAS ) ) /AM AS AF ) 

ClFSl=8ETA**2*(2.*( l.-COAFj-I { 1 • -CO AF ) -AMI *( COASAF-COAS ) )**2/ 

1 AMASAF > 








A = D0FS5**2*AFS14-D0FS5*8FSI+CIFS1 
IF (A) 10 , 20,20 

1 0 CONTINUE 

PRINT 100 , I ASAR , I AFAR , A 

ICC FORMAT I//,*SU3ROUT INE SDT5Sl,SQ. ROUTE ARGUMENT NEGATIVE *,3X, 
l'FALFAS = * : , 13 , 3X r^ALFAF = *, 13 , 3X , #-SQRT . ARG* =*,F10*5/> 

GO TO 40 
20 CONTINUE 

DTSSI=SORT( A) 

40 CONTINUE 
RETURN 
END 


SUBROUTINE SDT1S2 
DIMENSION DEL T A { 5 , 5 ) 

COMMON /SSFFS/ 

1 COAS, C0AF,C0ASAF, AMAS,AMAF, AMASAF, AM I , VF , GAMAS, GAMAF, BETA,OELTA, 

2 I AFAR , I ASAR , DOSF, DOFS , DOSS T DOFF , INDXSF, INDXFS, INDXSS , INDXFF, 

3TSF, TFS,TSS ,TFF 
4 , DELTSS , DELTFF , D ELTSF , DELTFS 
COMMON /FS/ 

I OOF S 1 ,DQFS2,DCFS4,DQFS5,T1 SI ,T2S1 ,T4SI , T 1 S2 , T2S2 , T4S2 , DT2 S 1 , DT 1 S2 
2,T5S1,DT5S1 
A-=I *-( COAS-AM I ) **2/AMAS 
IF (A) 10,20,20 

10 CONTINUE 

PRINT 100, I ASAR, I AFAR, A 

100 FORMAT!//, ^SUBROUTINE SDT1S2,SQ„ ROUTE ARGUMENT NEGATIVE *,3X, 
1*ALFAS =* , 13, 3X,*ALFAF , I 3 , 3X , *SQRT . ARG. =*,F10 .5/ J 

GO TO 40 
20 CONTINUE 

DT1S2-=DOFS1*SORT{ A) 

40 CONTINUE 
RETURN 
END 


SUBROUTINE SDOSS 
DIMENSION D EL T A ( 5 , 5 ) 

COMMON /SSFFS/ 

1 COAS, COAF,COASAF, AM AS , AM AF, AMASAF, AM I , VF , GAM AS, GA MAF , BETA , DELTA , 
21 AFAR , I ASAR, DOSF , DOFS , D OSS , DOFF , I NDXSF , I NOXFS , INDXSS, INDXFF, 

3TSF , TFS , TSSjTFF ’ 

A , DELTSS, DELTFF , DEL TSF, DELTFS 
DOSS 1 =DELT SS/SQRT! . 5-K* 5 *C0 AS ) 

TlSl=D0SSl/{2 . #AMI*VFJ 
TS-G AMAS/ ( AMI^VF ) 






' i ."N*#**' VW - ,T-. ■_ -** r 
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!'~ {T1S1-T5) 10,10,110 
10 CONTINUE 
DOSS=OOSSl 
INDXSS=1000 
GO TO 200 
110 CONTINUE 
A SS5= 1 • 

BSS5=-2.*GAMAS*M .-COAS) 

C 1 S S 5= 2 G A M A S* *2'-* ( 1 .-COAS ) 
CSS5=C 1 SS 5— DEl_TSS 3 !' , i c 2 
DETSS5=3SS5**2-4.=i=CSS5 
D O S S 5 = ( - B S S 5 + S Q R T ( O E T S S 5 ) ) /2 . 
00SS-00SS5 
INOXSS-AOOO 
200 CONTINUE 
RETURN 
END 


SUBROUTINE SDOFF 
DIMENSION OELT A { 5 , 5 ) 

COMMON /SSFFS/ 

1 COAS, COAF, COASAF, AM AS , AMAF, AM ASAP, AMI , VF, GAMAS, GAMAF, BETA , DELTA , 
21 AFAR, IASAR, OOSF, DOFS, DOSS, DOFF, INDXSF, I NDXFS , I NDXSS , I ND XFP , 

3T SF ,TFS,TSS, TFF 
4 , OELTSS,DELTFF ,DELTSF , DELTAS 
DOFF 1 = DEL TFF /SORT ( • 5+ • 5 *CO AF ) 

T IS I =DOFF 1/(2 »vVF) 

T5=GAMAF/VF 
IF (T1S1-T5) 10,10,110 

10 CONTINUE 
OQFF = DOFF 1 
I NDXFF= 1000 
GO TO 200 
110 CONTINUE 
AFF5=1. 

QFF5=-2.*GAMAF*{ 1 .-COAF) 

C 1 FF 5= 2 .FG AMA FFF2 y - c I 1 .-COAF ) 

CFF5=CIFF5-0ELTFFF*2 
DETFF5=3FF5#*2-4.*CFF5 
DOFFS= { — BFF5 + SQRT t DETFF5 ) )/2 • 

D0FF=D0FF5 
I NDXFF= 4 00 0 
200 CONTINUE 
RETURN 
END 
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SUBROUTINE DOT 
DIMENSION DELTA! 5,5) 

COMMON /SSFFS/ 

1 CO AS* CDAf-', COASAF, AM AS * AMAF, AMASAF, AMI , VF , GAMAS* GAMAF , BETA , DELTA . 
21 AFAR, IASAR,DOSF, DOF S , DOSS , DOFF , INOXSF, INOXFS, INDXSS. INDXFF. 

3TSF , TF S , TSS , TFF 
4 , DELTSS , DELTFF r DELTSF , DELTFS 
TSS=DOSS/( A.MI*Vr ) 

TFF=OOFF/VF 

TSF= { GAMA S-f-DO SF ) / VF— G AM AS/ ( AM I* VF ) 

TF3=< GAMAS+DOFS) / { A MI * VF ) -GAMA S/ VF 

RETURN 

END 
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PROGRAM COMPUTING CAPSF 

PROGRAM CAPSF C INPUT, OUTPUT ) 

DIMENSION I SARI 10 , 19) , I FAR! 10, 19),D00SF( 10,19), INXSF{ 10,19), 
1D00FSI 10, 19) , INXFSI 10,19) , DQ0SS! 10) , INXSS(IO) ,DQ0FF{ 1 0,19) , 
2INXFFM0, 19) , CAP AC I 10,19), 

3DELT A < 5 , 5 ) 

COMMON /SSFFS/ 

1 CO AS , CO AF:, CO AS AF , AM AS , AMAF , AM A SAF , AM I , VF , G AM AS , G AMAF, BETA , DELT A , 
21 AFAR, I AS AR , DOSF , DOFS, DOSS ,D0FF s I NDXSF , I NOXFS , I NDXSS , INDXFF, 

3TSF , TFS , TSS , TFF 
4 , DELT SS, DEL TFF , DELT SF , DELTFS 

DO LOOP TO REPEAT THE WHOLE PROGRAM WITH NEW INPUT DATA 

9999 CONTINUE 

INPUT DATA 

. PALF-1.5707? 

DALF=0 • 1T4 53 

IPALF=90 

IDALF=10 

3010 FORMAT! 2F1 0.5 ) 

3015 FORM AT (4F10.5) 

READ 30lO,VFKT,VSKT 
READ 3010 , GAMAF, GAMAS 

READ 3015, DEL TSS , OELTSF , OELTFS , DELTFF 
READ 3010, PF , PS 

INITIAL DATA CONVERSION 

VS=VSKT/3600. 

V F= VFKT/3600. 

AM 1= VS/ VF 

BET A— GAMAF— GA MAS 

PFF=PF*PF 

PSS=PS*PS 

PFS=PS*PF 

PSF=PS*PF 

PRINT INPUT DATA 

PRINT 3110, VFKT 

3110 FORMAT I 1HI , 1 0X,*VF -*,F7.2/) 

• PRINT 3120, VSKT 
3120 FORMAT I 1 IX,*VS =*,F7.2/) 
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PRINT 3130, AMI 

3130 FORMAT (11 X , *M I =*,F10.5/) 
PRINT 314 0* GAMAF 
3140 FORMAT { 8X * ❖GAMAF s*,F7.2/} 
PRINT 3150 » GAMAS 
3150 Format (8X ,*gamas =*,F7.2/> 
PRINT 3160, BETA 
3160 F0RMATC9X,*3ETA =*,F7.2/) 
PRINT 317 1 * DELTSS 
PRINT 3172, DELT SF ' 

PRINT 3173, DELTFS 
PRINT 3 1 74 * DELTFF 

3171 F0RMATI7X , *DELTSS =*,F7.2/) 

3172 F0RMAT{7X,*DELTSF =*,F7,2/) 

3173 FORMAT < 7X t^OELTFS =*,F7.2/) 

3174 FORM ATI 7 X , MQELTFF =*,F7.2/) 
PRINT 3180, PF 

PRINT 3190 , PS 

3180 FORMAT! 1 1 X,*PF =*,F7.2/> 
3190 FORMAT ( 1 1 X,4PS =*,F7,2/> 

ANGLE MATRIX GENERATION 

DO 999 1=1,10 
EI1 = FL0AT( 1-1 ) 
AS=PALF-DALF*EI 1 
l AS AR= I PALF — I DALF^I I-i ) 

DO 998 J=I , 19 
IF ( I — J ) 2120,2110,2120 
2110 CONTINUE 

AF=AS-0 •5 > i t DALF 
I AF AR= IAS AR— I DALF/ 2 
GO TO 2130 
2120 CONTINUE 

EJi=FLOAT( J-l ) 
AF=PALF-DALF*E J1 
r AFAR=IPALF-.IDALF*( J-l ) 

2130 CONTINUE 

I SAR ( I , J)=IASAR 
IFARI I , J ) = I AF AR 

PREPARATION FOR SUBROUTINES 

COAS=CQS( AS) 

CGAF=CQS(AF) 

COASAF=CQS( AS-AF) 
AMAS=AMI**2. -2.*AMI *COAS+l . 
AM AF= AM 14 * 2 • - 2 .* AMI *CO AF 4- 1 . 
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AMASAF»AMI**2*-2**AMI*C0ASAF+l 


MAIN PROGRAM 

CALL SOOSF 
OGOSF< I , J )-OOSF 
I NX SF( I , J)=INDXSF 
CALL SDOFS 
OOOFS ( I , J )=DGFS 
INXFSI I , J )=INDXFS 
CALL SDOFF 
OOOFF ( I , J >-DOFF 
I NXFF ( t • J ) = l NDXFF 
CALL SDOSS 
DOOSS ( I )=DOSS 
I N XSS ( I 1= INDXSS 
CALL DOT 

T BAR =T S5%t-P.SF + T FS^PFS+TFF^PFF+TSS 4PSS 
CAP=36O0. /T0AR 
CAP AC( l, J)=CAP 

998 CONTINUE 

999 CONTINUE 

PRINT TABLES 

PRINT 4700 
4700 FORMAT ( 1 H 1 , 20 X, 4CAPSF*//// ) 

PRINT 4805, MSARI I* J), 1 = 1 ,10) 

DO 4710 J=l,19 
I ANGF=100-1 0*J 

PRINT 4711, I ANGF, (CAPAC{ I , J > * 1 = 1 , 1C) 
4711 FORM AT I I 6 , 4X , 1 OF 1 0 • 3/ ) 

.4710 CONTINUE 

PRINT 4800 

4800 FORMAT! 1 H 1 , 20 X , *DOSF*//// ) 

PRINT 480 5 , (I SARI 1 , J> , 1 = 1 , 10 > 

DO 4310 J= 1 , 19 
I ANGF=100-10*J 

PRINT 48 11,1 ANGF , ( DOOSF (I, J), 1=1,10) 
4810 CONTINUE 

PRINT 4820 

4320 FORMAT ( IH1 , 20 X ,*I NDXSF*//// ) 

PRINT 430 5, < I SARI l, J), 1=1,10) 

DO 4830 J=1 , 1 9 
I ANGF=10Q-10*J 

PRINT 4812, I ANGF, I INXSF( I, J) , 1=1 , 10) 
4330 CONTINUE 
4805 FORM AT { 1 0 X , I 0 1 1 0 / / ) 
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4811 FORMAT ( I 6 , 4X , 1 OF 1 0*5/) 

4812 FORMATX I6,4X, 10I10/J 
PRINT 4840 

4840 FORMAT! 1 HI , 20 X , *DVFS#/SS / ) 

PRINT 4805, ( TSARfT , J) ,1=1 ,10) 

00 4850 J=l, 19 
I ANGF=100-10#J 

PRINT 4811,1 AN GF , IDQQFSI I ,4), 1=1 ,10) 
4850 CONTINUE 

PRINT 4660 

4860 F 0 RM AT ! 1 H 1 » 20 X , * I NO XF S * //// ) 

PRINT 4805, X I SARI I , J ) , 1 = 1 , 10 ) 

DO 4870 J=l,19 
I ANGF= 1 00-1 Q*J 

PRINT 4812 , I ANGF, ( INXFS( I, J), 1=1, 1C) 
4370 CONTINUE 

PRINT 4880 

4880 FORMAT < IH1 ,20X ,*DOSS, INDXSS*//// ) 
PRINT 4305, ( I SAR( I , J) , 1=1 ,10) 

PRINT 4831 .XDOOSS! I ), 1=1 , 101 

4881 FORMAT! 1 OX, 1 0F1 0 *5/ > 

PRINT 4882, X INXSSX I ),I=1 ,10) 

4382 FORMAT ( VOX, 10 I 10) 

PRINT 4890 

4890 FORMAT <1H1 , 20 X , *DOFF*//// J 

PRINT 4805, { I SAR! I, J), 1=1, 10) 

DO 4900 J=l , 19 
I ANGF=100-10«J 

PRINT 4811 , I ANGF, ( DOOFF I I , J ) , 2 =1 , 1 0 5 
4900 CONTINUE 

PRINT 4910 

4910 FORMAT { 1H 1 , 20 X , * I NO XFF *//✓/ ) 

PRINT 4805, USAR! I, J) ,1=1, 10) 

DO 4920 J=1 , 19 
I ANGF= 100-10# J 

PRINT 481 2, IAN GF , { I NXFF C I , J ) , 1=1 , 10) 
4920 CONTINUE 

GO TO 9999 

STOP 

END 


B.4. Program CAP 3 (Case with three air craft types). 



C PROGRAM COMPUTING CAP3 

C . . . 

PROGRAM CAP3< INPUT, OUTPUT) 

DIMENSION VKT<5) , G AM A { 5 ) , P ( 5 ) , PP (5, 5 > , DELTA ( 5,5), TT <5,5), 

1 I ALF 1(50) , I ALF2( 50 > , I Al>3< 50 ) , CAPAC( 50 ) , CAPMI N< 50 ) , 

2 1 ALFM 1 { 50 ) , I ALFM2 ( 50 ) , I AUFM3 < 50 ) 

common /ssffs/ 

1C0AS,C0AF,C0ASAF, A MAS, AMAF, AMASAF, AMI , VF, GAMAS,GAMAF, BETA, DELTA, 
21 AFAR, I ASAR, DOSF, DOFS , OOSS , DOFF, I NDXSF , I NDXFS , I NDXSS , I ND XFF , 
3TSF,TFS,TSS,TFF 
4, DELTSS, DEL IFF , DELT SF, DELTFS 

C DO LOOP TO REPEAT THE WHOLE PROGRAM WITH NEW INPUT DATA 

C ' > . , 

9999 CONTINUE 
C 

C INPUT DATA 

C ■.■!■! ' - ; . ■;:.■■■■ 

M — 3 
- NM-40 

READ 30 10, ( VKT ( I ) ,1=1 , M ) 

3010 FORMAT (3F10.5) 

READ 3010, { GAM A( I ) , 1=1 ,M) 

READ 3010, (P( I ) , 1 = 1 ,M> 

READ 3010, {{DELTA! I, J) , J— 1 , M ) ,1— 1 » M) 

C 

C PRINT INPUT DATA 
C 

PRINT 3111 

3111 FORMAT ( 1H1 , 1 OX,* INPUT DATA*,///) 

DO 3110 1=1, M 

PRINT 31 12, I , VKTC I ) , I, GAMA ( I ) , I ,P( I ) 

3112 FORM AT( 11X,*VKT*, I 1 , * =*, F7. 2 , 3X , *GAMA*, II,* =*, F7. 2, 3X, *P*, 1 1 , 
1* =*,F4.2/) 

3110 CONTINUE 

PRINT 3120 

3120 FORMAT!//, 1 IX, *DELTA MATRIX*//) 

PRINT 3130 - 1 -- — : ■ 

3130 F0RMAt(20X,*TRAILING A/C*/) 

PRINT 3131 , ( J, J=l , M) 

3131 FORMAT ( 19X ,3 ( II ,9X) ,/) 

, DO 3140 I — 1 , M 

PRINT 3132, I, ( OELTAC I , J ) , 3=1 , M) 

3132 FORMAT ( 9X , I 2, 3F 10*2,/) 

3140 CONTINUE 

C . . 

C INITIAL DATA PREPARATION 
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DO 4010 I — 1 * M 
DO 4009 J = 1 ,M 
PP( I , J )=P( I >*P< J) 

4009 CONTINUE 

4010 -CONTINUE 

DO 4020 N=1,NM 
C AP ACT N } = 0 • 

C AP M I N { N ) — 99 99* 
t ALF1 <N) = 0 r-T 
I ALF2 ( N} = 0 
I ALF3(N)=0 
4020 CONTINUE 

PALF=1 .57077 


DALF=0. 1 


7453 


I PALF=90 


I DALF= 1 0 



MAIN PROGRAM 


CASE A 


ANGLES GENERATION 


DO 2495 t 1-1 , 3 

Ell -FLOAT < 11-1 ) 

Al=PALF-DALF4Ei 1 

I A1 AR= IPALF-IDALF*{ 1 1-1 ) 

I IPLSI = I l + l 

DO 2490 t 2= 1 1 PLS1 « 9 

E21— FLOAT ? 12-1 ) 

A2=PALF-DALF*E21 

I A2AR=IPALF-T0ALF*{ 12-1 ) 

I 2PL51— 12+ 1 

DO 2485 I 3= I 2 PL SI * 19 

E31=FLOAT{ 13-1 J 

A3=PALF-DALF*E31 

I A3AR=IPALF-IDALF*( 13-1) 

AIRCRAFT 1 AND 2 

VS=VKT{ 1 ) /3600* 

VF=VKT( 21/3600. 

AM I = V S/VF 
GAMAF=GAMA(2) 

GAMAS=GAMA{ 1 ) 

BETA=GAMAF— GAMAS 
DEL TSS= DELTA! 1,1) 
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DELTFF-0ELTAI2 ,2) 

DELTSF=DELTA< 1 ,2 ) 
DELTFS=DELTA<2,1 ) 

AS-Ai 

AF=A2 

IAFAR=IA2AR 
I AS AR= I A1 A!R 
CQAS=COS( AS) 

COAF=COS< AF) 

C0ASAF=:CQ5 ( AS-AF) 
AMAS=AMI**2-2.*AMI*C0AS+1 • 
AMAF=A.M IF*2-2 .*AMI*C0AF*1 • 

AM AS AF= AM I **2-2.* AMI *‘CQ ASAF + I • 
CALL SDOSF 
CALL SDOFS 
CALL SDOFF 
CALL 500SS 
CALL DOT 
TT ( 1*1) =TSS 
TT( i f 2 )=T SF 
• TT ( 2 « 1 ) =TFS 
TT<2,2)=TFF 

AIRCRAFT 1 AND 3 

VS=VKT( 1 ) /36C 0 « 

VF=VKT{ 3 ) /3600. 

AMI=VS/VF 

GAMAF=GAMA(3) 

G AM AS=GAMA< 1 > 

B ET A= G A M AF -G AM AS 
DELTSS— DELT A { 1,1) 

0ELTFF = DELTA(3 ,3 ) 

DELTSF— DELT A { 1 ,3) 

DEL TF S=OELT A ( 3 » i ) 

AS=A1 

AF=A3 

I AFAR- I A3 AR 
I ASAR-I ATAR 
CQAS-COSI AS) 

COAF=CDS( AF) 

COASAF=COS< AS-AF) 

AMAS=AMI**2-2 » #AMI*COAS-l-l • 

A M AF= A M I * * 2-2 • * AM 1 4 CO AF + 1 • 
AMASAF=AMI **2-2.*AMI*C0A$AF+l • 
CALL SDOSF 
CALL SDOFS 
CALL SDOFF 
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CALL SOOSS 
CALL DOT 
TT ( 1 » 3 ) =T SF 
TY£ 3 f 1 1— TFS 
TT(3,3)=TFF 

AIRCRAFT : 2 AND 3 

VS=VKTt 2J/3600. 

VF— VXT t 3 } /360 0 • 

AM I = VS/ VF 
GAMAF=GAMA (3 ) 

G AM AS=G AMA ( 2 ) 

BET A = GAMAF — GAMAS 
DELTSS= DELTA! 2,2) 

DELTFF=DELTA< 3 f 3 > 

DELTSF=DELTA( 2, 3 ) 

DELTF S^DELT A { 3 , 2 ) 

AS- A 2 
AF=A3 

- I AS AR= I A2AR 
I AF AR = I A3 AR 
CaAS=COS(AS) 

CQAF=COS( AF) 

COASAF=CaS C AS-AF) 

AMAS=AMI**2-2 .*AM| SCOAS+l • 
AMAF^AM I**2-2.*AMI*CaAF+l • 
AMASAF=AMI**2-2.*AMI*C0ASAF+1 • 
CALL SDOSF 
CALL SDOFS 
CALL SDOFF 
CALL SDOSS 
CALL DOT 
TTt 2,3) ~T SF 
7T< 3, 2 ) = TFS 

MEAN INTERARRIVAL TIME AND CAPACITY 

T3AR=0. 

DO 2420 I — 1 , M 
00 2410 J=1,M 
TBARsTBARFTTC I, J) *PP ( 1 , J ) 

2410 CONTINUE 
2420 CONTINUE 

CAP=3600./TBAR 
N = NM 

IF ICAP.LT. CAPAC( N) ) GO TO 2470 
DO 2450 Nl-ItN ■ 
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IF (CAP-CAPAC<N1) } 2450,2450,2440 
2440 CONTINUE 
NN=N1 

GO TO 2455 
2450 CONTINUE 
2455 CONTINUE 

N MI NNN= N—NN 

DO 2460 N0P= 1 , NMI NNN 

N2-N—NOP 

N2Pl=N2+l 

CAP AC ( N2P I ) — CAPACt N2 ) 

I ALF1 (N2P1 ) = IALFHN2J 
I ALF2< N2P 1 ) = I AL.F2 ( N2 ) 

I ALF3(N2Pl> -IALF3IN2) 

2460 CONTINUE 

CAPACI NN)=CAP 
I ALF 1 { NN J = I A I AR 
I ALF2C NN)=IA2AR 
I ALF3(NN>=I A3AR 
2470 CONTINUE 
• N = NM 

IF <CAP.GT.CAPMIN£N> > GO TO 2478 
DO 2474 Nl=l f M 

IF ( CAP-CAPMININl ) ) 2472,2474,2474 
2472 NN=N 1 

GO TO 2475 

2474 CONTINUE 

2475 CONTINUE 
NMI NNN=N—NN 

DO 2476 NOP= 1 , NMINNN 

N2= N— NOP 

N2PI=N2+1 

CAPMIN(N2P1 )-CAPMIN(N2) 

I ALF Ml { N2P1 >=I ALFM1 ( N2 ) 

I AL* r M2{ N2P1 ) = I ALF, M2 ( N2 ) 

I ALFM3 ( N2P1 >= I ALFM3 £ N2 ) 

2476 CONTINUE 
CAPMIN(NN)=CAP 

I ALFMi ( NN)= I A1 AR 
I ALF M2< NN ) - 1 A2AR 
IALFM3(NN)=I A3AR 
2478 CONTINUE 
2435 CONTINUE 
2490 CONTINUE 
2495 CONTINUE 

MAIN PROGRAM 
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CASE 8 

ANGLE GENERATION 

00 2595 11=1,9 
Etl=FL0AT< 11-t) 

A 1 =PALF-0 ALF*E1 1 

1 A 1 AR = f PALF— lOALF^C 11— l 1 
I 1PLS1=I 1+1 

DO 2590 12=11,19 
E2l=FL0AT< 12-1 1 
A2=PALF-DALF*E21 
I A2AR=I PALF— IOALF*< 12-1 1 
I 2MINI= £2—1 

00 2585 I 3= 1 1 PLS1 , 1 2M INI 
E31=FL0AT ( 1 3— 1 1 
A3=PALF-DALF*E31 
I A3AR=tPALF-IDALF*< 13-1 1 

AIRCRAFT 1 ANO 2 

VS=VKT< U/3600* 

VF=VKT 12) /3600* 

AMI= VS/VF 
GAMAF=GAM A1 2 1 
GAMAS=GAMA< 1 ) 

BETA=GAMAF-GAMAS 
DELTSS=DELTA< 1,1) 

0ELTFF=0ELTAC 2 ,2 1 
DELTSF=DELTA< 1 ,21 
0ELTFS=0ELTA<2,1 1 
AS*A 1 
AF=A2 

SAFAR=IA2AR 
IASAR=1A1AR 
C0AS=C0S( AS) 

C0AF=C0SC AF) 

COASAFssCOSC AS— AF) 
AMAS=AMI**2-2«*AMI*C0AS+1 • 
AMAF=AM1**2— 2 • *AMI*C0AF+1« 

AMAS AF= AM I«*2-2**AM I ♦COASAF+l . 

CALL SDOSF 

CALL SOOFS 

CALL SDOFF 

CALL SOOSS 

CALL DOT 

TTC 1 , i )=TSS 

TT( 1 ,2)=TSF 
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TT (2,1 1— TFS 
TTC 2 *2 >— TFF 

AIRCRAFT 1 AND 3 

VS=VKT{ 1 ) /3600 • 

VF=VKT(3)y3600* 

AM1=VS/VF 
GAMAF=GAMAC 3) 

GAMAS=GAM AI 1 ) 

BET A=GAMAF— GAMAS 
DELTSS=OELTA( 1 , 1 ) 
DEL7FF=DELTA(3 ,3) 

OELTSF=DELTA< 1,31 
DELTF S-=DELT A! 3 , 1 > 

AS— A 1 
AFss A3 

I AFAR= IA3AR 
I ASAR-1A1AR 
CO AS- COS! AS ) 

-COAF=COS! AF) 

COASAF=COS( AS-AF) 

AMAS=AMI**2-2 .*AMI*COAS+l . 
AMAF=AMI**2-2 •❖AM^COAFM • 

AMAS AF= AM I **2-2 **AM I FCQASAF-H • 

CALL SDOSF 

CALL SDOFS 

CALL SDOFF 

CALL SDOSS 

CALL DOT 

TT! 1 ,3)-TSF 

TT( 3, l)=TFS 

TT( 3 ,3) =TFF 

AIRCRAFT 2 AND 3 

VS=VKT( 2J/3600* 

VF=VKT!3>/3600* 

AMI — VS/VF 
GAMAF=G AM A (3) 

GAMAS=GAMA!2) 

8ETA=GAMAF— GAMAS 
DELTSS= DELTA! 2,2) 

OELTFF— DELTA! 3,3) 

DELTSF— DELTA! 2 ,3 ) 

OELTFS^ DELTA! 3,2) 

AS=A2 

AF=A3 
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I ASAR=I A2AR 
I AFAR-IA3AR - 
COAS=COS( AS) 

COAF=COS( AF) 

COASAF=COSI AS-AF) 
AMAS=AMI*42-2«4AMI*CQAS+1* 
AMAF=AMI**2-2**AMI*C0AF«-1 . 
AMASAF=AM I 4*2-2 . * AM I *C0 ASAF+ l * 
CALL SDOSF 
CALL SDOFS 
CALL SODFF 
CALL SDOSS 
CALL DOT 
TT(2,3)=TSF 
TT{3,2)=TFS 

MEAN INTERARRIVAL TIME 

TBAR-O* 

DO 2520 1 = 1 ,M 
-DO 2510 J= 1 » M 
T8AR =TBAR+TT{ I,J)*PP< I , J) 

2510 CONTINUE 
2520 CONTINUE 

C AP=360 0 • /TBAR 
N=NM 

IF < C AP «LT« CAP AC{ N ) ) GO TO 2570 
00 2550 Nl= 1 , N 

IF T CAP— CAPACI N1 ) ) 2550,2550,2540 
2540 CONTINUE 
NN=N1 

GO TO 2555 
2550 CONTINUE 
2555 CONTINUE 

NMINNN=N-NN 

DO 2560 NOP=l , NMINNN 

N2=N— NOP 

N2PI=N2+1 

C AP ACIN2P1 )=CAPAC<N2> 

I ALF 1 { N2P 1 ) = I ALF 1 ( N2 ) 

I ALF2< N2P 1 )=I ALF2CN2) 

I ALF3CN2P1 )=IALF3(N2> 

2560 CONTINUE 

CAPACI NN)=CAP 
I ALF I < NN)=I A1 AR 
I ALF2(NN)=I A2AR 
I ALF3( NN)=I A3AR 
2570 CONTINUE 
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N=NM 

IF (CAP.GT.CAPMININ)) GO TO 2578 
DO 2574 NT— 1 , N ■ 

IF ( CAP-CAPMIN(NI )S 2572,2574,2574 
2572 NNsNl 

GO TO 2575 

2574 CONTINUE 

2575 CONTINUE ! - = L - 

NMINNN=N-NN 

00 2576 NOP— 1 , NMINNN 
N2=N— NOP 
N2PI=N2+1 

CAPMI N CN2P1 )— CAPMI N (N2) 

1 ALFM1I N2P1 )=IALFMKN2) 

I ALFM2IN2P1 )=I ALFM2IN2) 

IALFM3(N2P1)=IALFM3(N2> 

2576 CONTINUE 
CAPMI N ( NN ) =C AP 

I ALFM I (NN )— I Al AR 
I ALF M2 ( NN ) = I A2 AR 
- I ALF M3 1 NN ) =1 A3 AR 
2578 CONTINUE 
2585 CONTINUE 
2590 CONTINUE 
2595 CONTINUE 

PRINTING RESULTS 

PRINT 4050 •' 

4050 FORM ATI 1H1 ,///, 13X,*MAXIMAL CAPACITY-OPTIMAL ANGLES*,////) 
PRINT 4100 

4100 FORMAT M4 X , *CAP ACI TY4 , 5X , * ALF A 1 * , 5X » * ALF A2* , 5X , *ALFA3*, ///) 
DO 4150 NP= 1 , NM 

PRINT 4151 ,NP, CAPACCNP) , IALF1 <NP1 , I ALF2( NP>, IALF3INP) 

4151 FORMAT < 5X , I5,Fl 0*3,31 10) 

4150 CONTINUE 

PRINT 4250 

4250 FORMAT < 1H 1 ,/// ,13X,*MINI MAL CAPACITY-THE WORSE ANGLES*,////) 
PRINT 4100 
DO 4260 NP— 1 , NM 

PRINT 4151, NP, CAPMI N( NP ) , I ALFM1 ( NP I , IALFM2IMP), I ALF M3 (NP) 
4260 CONTINUE 

GO TO 9999 ! 

STOP ■ - -—i. 

ENO L_ 
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PROGRAM COMPUTING CAP4 

PROGRAM CAP4( INPUT, OUTPUT) 

DIMENSION VKT (5 ) , GAMA { 5 ) ,P! 5> ,PP{5,5) , DELTA (5,5) ,TT( 5,5) , 

1 I ALFl ( 50) , I ALF2 ( 50 ) , IALF3! 50) , I ALF4( 50 ) , CAP AC ( 50 ) , CAPMIN! 50) , 

21 ALFM1 ( 50) , I ALFM2( 50 ) , I ALFM3! 50 ) , I A LF M 4 { 50 ) 

COMMON /SSFFS/ 

1C0AS, COAF , CO AS AF, AM AS , A MAF , AM AS AF, AMI , VF , GAMAS, GAMAF , BETA , DcLTA , 
2IAFAR,IASAR,DQSF, DOFS, DOSS, DOFF, INDXSF, INDXFS, INDXSS , INDXFF , 
3TSF, TFS,TSS,TFF 
4 , DELTSS , DELTFF , DELTSF , DELTFS 

DO LOOP TO REPEAT THE WHOLE PROGRAM WITH NEW INPUT DATA 

>99 CONTINUE 

INPUT DATA 

M=4 
NM=4Q 

READ 30i0,(VKT( I), 1 = 1 ,M) 

>10 FORMAT! 4F 10*5) 

READ 3010, (GAMAC I ) , 1=1 ,M) 

READ 3010 , (PI I ) ,1-1 ,M) 

READ 3010, ( (DELTA! I ,J> ,J=1 ,M) ,1=1, M) 

PRINT INPUT DATA 

PRINT 3111 

3111 FORMAT! 1 H 1 , 1 0 X , * I NPUT DATA*,///) 

DO 3110 1=1, M 

PRINT 3112, I ,VKT( I ) , I , GAMA ( I > , I,P( I) 

3112 FORMAT! 11 X, 4 VKT#, tl ,* =*,F7*2 , 3X,*GAMA*, ( 1 , * s*,F7,2 f 3X,*P*,U f 
1* =*,F4.2/> 

3110 CONTINUE 

PRINT 3120 

3120 FORMAT <//, 1 IX, *DELTA MATRIX#//) 

PRINT 3130 

3130 F QRM AT (21 X,# TRAILING A/C*//) 

PRINT 3131 , ( J, J=1,M) 

3131 FORMAT! 13X, 4! II ,9X) ,/) 

DO 3140 1=1, M 

PRINT 3132, I, (DELTA! I, J), J=1 ,M) 

3132 FORMAT (7X, I2,2X,4(FI0*2),/) 

3140 CONTINUE 

INITIAL DATA PREPARATION 
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DO 4010 1=1 ,M — 
DO 4009- J=1 ,M 
pp( i ,j)=p( i )*p(ji 

4009 CONTINUE 

4010 CONTINUE 

DO 4020 N-l t NM 
CAPAC(N)=0. 
CAPMIN{N)=9999* 
IALF1 CN)=0 
IALF2(N)=0 
I ALF3 ( N } = 0 
I ALF4 ( N) = 0 
4020 CONTINUE 

PALF=1 ,57077 
. DALF= 0.1 7453 
V IPALF=90 
IDALF=10 


MAIN PROGRAM CASE A 

ANGLES GENERATION 

DO 2195 11=1,7 
Ell =FLOAT ( 1 1-1 ) 

A1 = PALF-DALF*EU 
I A1 AR=I PALF— I DALF* { Il-l 1 
I 1PLS1=I 1+1 
DO 2190 I2=11PLS1,3 
E21=FLOAT( I 2— I ) 

A2=PALF— DALF4E21 
I A2 AR= I PALF— I DALF4( 12-1 ) 
I 2PL S 1= 1 2+ 1 

00 2185 I 3= I 2PLS1 , 9 
E 31 =FLQ AT ( 13-1 ) 
A3=PALF-DALF*E31 

1 A3AR=IPALF— IDALF*< 13-1 ) 
I 3PLS 1 =13+1 

DO 2180 I 4= 1 3PLS1 ,19 . 

E41=FLOAT( 14-1 > 

A4=PALF— DALF^EA l 
I A4AR = 1 PALF-IOALF*( 14-1 ) 

AIRCRAFT 1 AND 2 

VS=VKTI 1 >73600. 

VF=VKT< 2J/3600. 
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AMI=VS/VF 

GAMAF-=GAMA(2) 

GAMA'5=GAMA{ 1 ) 

BET A=GAMAF— GAMAS 
DELTSSsDELTAC 1,1 ) 
DELTFF=DELTA(2,2) 

OELTSF=DELTA( 1 ,2) 

DEL TFS= DELTA (2,1 ) 

AS= A 1 
AF-A2 

I AFAR=IA2AR 
I AS AR= I A 1 AR 
CO AS=COS( AS} 

COAF=CQS(AF) 

CO A SAF = COS C AS— AF ) 

AMAS=AM I**2-2 . *AMI *COAS+l . 

AM AF= AM 1**2— 2 • *AM I *COAF + 1 • 

AMAS AF = AM I *>1=2-2 • * AM I *COASAF+ 1 . 
CALL SDOSF 
CALL SDOFS 
' CALL SOOFF 
CALL SDOSS 
CALL DOT 
TTC 1 t 1 ) — TSS 
T T { 1 , 2 > =T SF 
TT{2,1)=TFS 
TT{ 2 ,2) =TFF 

AIRCRAFT 1 AND 3 

VS=VKT{ 1J/3600. 

VF= VKT (3)/36Q0* 

AMI.— VS/VF 
GAMAF=GAMA(3) 

GAMAS=GAMA< 1 ) 

BET A = GAMAF — GAMAS 
DELTSS= DELTA { 1 ,1 ) 

DELTFF = DELTA(3 ,3) 

DELTSF=sDELTA< 1 ,3) 
DELTFS=0ELTA(3 , 1 ) 

AS=A1 

AF=A3 

l AFAR=IA3AR 
I AS AR= I A 1 AR 
CO AS— COS { AS ) 

CO AF = COS ( AF ) 

COASAF=COS( AS-AF) 

AMAS=AM I**2-2 • *AMI *COAS+l • 
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AMAF=AMI**2-'2.*AMI*CaAF+l . 

AMASAF= AM I #*2— 2* * AM I =5*00 A SAF 1 • 

CALL SOOSF 

CALL SOOFS 

CALL SQGFF 

CALL SOOSS 

CALL DOT 

TTC 1 , 3)=TSF 

TT ( 3 » 1 ) =TFS 

TTi3,3)-TFF 

AIRCRAFT 1 AND 4 

V S= V K T { 1 ) / 3 6 0 0 • 

VF=VKT{ 4) /3600. 

AMI-VS/VF 
G AM A F— GAM A { 4 ) 

GAM AS=GAM A { 1 ) 

BETA=GAMAF— GAMAS 
DELTSS = DELTAI 1 , 1 ) 

* DELTFF-DELTAI4 ,-4 ) 

DELTSF=0£LT A ( 1 ,4) 
DELTFS=DELTA<4,1 ) 

AS= A 1 
AF=:A4 

I AFAR-IA4AR ! 

IASAR=IA1AR 

COAS=COS(AS) 

CQAF=COS{ AF> 

COASAF=COS( AS-AF) 
AMAS=AMI**2-2**AMI*C0AS4-1# 

A M A F = A M I * * 2 - 2 * * A M I * C 0 A F ♦ 1 . 

AMASAF=AKI*#2-2.*AMI*C0ASAF*1 . 

CALL SDOSF 

CALL SDGFS 

CALL SDOFF 

CALL SDOSS 

CALL DOT .. 

TT{ l ,4)=TSF 
TT(4, 1 )=TFS 
TT( 4 , 4 > = TFF 

AIRCRAFT 2 AND 3 

VS=VKT( 2) /3600 • 

VF=VKT< 3)/3600. 

AMT=VS/VF 
GAMAF=GAMA( 3) 
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GAMAS*GAMA<2J 
BETAs GAM AF— GAMAS 
0ELTSS=0ELTA(2,2> 

DELTFFsOELTA{ 3,3) 
0ELTSF=DELTA(2,3) 
DELFSsOELTA(3,2) 

AS=A2 

AF=A3 

I AS ARs I A2AR 

I AFAR = T A3 AR . — t- ■ 

CQAS-COS(AS) 

COAF=COS( AF) 

COASAF=COS( AS-AF) 

AM A S— AM I **2—2 • * AM I *CQA SFl « 
AMAF=AMI**2-2 .*AMI*COAF+l . — - 
AMASAF = AMI**2-2.*AMI*C0ASAF-»-l • 
CALL SDOSF 
CALL SDQFS 
CALL SDOFF 

CALL SDOSS - — 

• CALL DOT 
TT( 2 » 3 J =T SF 
TT { 3,2) =TFS 
C 

C AIRCRAFT 2 ANO 4 
C 

VSsVKT { 2J/3600* 

VF=VXT ( 4 ) y 360 0 • 

AMIsVS/VF 
GAMAF=GAMA{ 4 ) 

GAMAS=GAMA<2) 

BETA= GAMAF-GAM AS 
DELTSSsDELTAf 2,2 ) 
0ELTFF=DELTA(4,4 ) 

D EL TSF= DELTA (2,4) 
DELTFS=DELTA(4,2) 

AS= A 2 
AF = A4 

I ASAR= I A2 AR 
I AFARs I A4AR 
COASsCOSCAS) 

COAF=COS(AF) 

COASAF=COS( AS-AF) 
AMAS=AMl**2-2.*AMI*CCAS+I . 
AMAFs AM I **2—2 • *AM I *C0AF *■ 1 . 
AMASAFsAM I **2-2. *AM I*COASAF FI • 
CALL SDOSF 
CALL SDQFS 
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CALL SDOFF 
CALL SDOSS 
CALL DOT 
TTC2 t 4)=TSF 
TT<4,2>=TFS 
C 

C AIRCRAFT 3 AND 4 
C 

VS=VKT< 3) /360C* 

VF= VKT { 4) /3600* 

AM I— VS/ VF 
GAM AF=G AM A ( 4 ) 

GAMAS=GAMA(3) 

BET A = GAMAF-GAMAS 
DELTSS— DELTA (3*3) 
DELTFF=DELTA(4,4) 

DELT SF= DELTA <3*4) 

DELTF S— OELT A < 4 *3 ) 

AS=A3 

AF=A4 

- I ASAR= I A3 AR 
I AF AR= I A4 AR 
COAS=COS( AS) 

COAF=COSCAF) 

CCAS AF=CO S { AS- AF) 

AMAS=AM I*»2-2 #*AMI 4C0AS + 1 • 
AMAF-AMI**2-2«*AMI*C0AF+1 . 

AM AS AF = AM I #42—2 • £ AM I ^CO AS AF + 1 • 

CALL SDOSF 

CALL SDCFS 

CALL SDOFF 

CALL SDOSS 

CALL DOT 

TT ( 3,4)=TSF '■ 

TT(4 # 3)=TFS 
C ■ 

C MEAN INTER ARRIVAL TIME 
C ■ 

TBAR=0 • 

DO 2120 I*1,M 
DO 2110 J=1 ,M 
TBAR=TQAR+TT( I , J)TPP( I , J) 

2110 CONTINUE 
2120 CONTINUE 

CAP=3600*/TBAR 

N=NM 

IF < CAP.LT.CAPAC(N) ) GO TO 2170 
DO 2150 N 1 — 1 |N 
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IF (CAP-CAP ACINI) ) 2150,2150*2140. 
2140 CONTINUE — - 

NN=N1 

GO TO 2155 ' -■ ■ ' . 

2150 CONTINUE . _ i 

2155 CONTINUE 

NMINNN^N-NN 

DO 2160 NQP=1,NMINNN 

N2-N-N0P 

N2P1=N2+1 

CAPAC (N2P1 J-CAPACC N2) 

CAP AC (N2P 1 ) =C APAC ( N2 ) 

I ALFKN2PI ) = I ALFKN2) 

IALF2(N2P1 )sIALF2(N2) 

I ALF3IN2P1 >=I Ai_F3{ N2 ) 

I ALF4 ( N2P 1 )= I ALF4 { N2 ) 

2160 CONTINUE 

CAPAC ( NN) =CAP 
I ALF1 ( NN ) = I A 1 AR 
I ALF2(NN)=rA2AR 
I AL.F 3 ( NN)=IA3AR 
1 ALF4< NN} = I A4AR 
2170 CONTINUE 
N=NM 

I F ( CAP »GT«CAPMIN(N) ) GO TO 2173 

00 2174 N1=1,N 

IF (CAP-CAPMIN(NIl) 2172,2174,2174 
2172 NN=N1 

GO TO 2175 

2174 CONTINUE 

2175 CONTINUE - 

N’MI NNN=N— NN 

DO 2176 NOP= 1 , NMT NNN 

N2=N-N0P 

N2P1 — N2+ 1 

CAP-MINI N2P1 ) =C APM I N ( N2 ) 

IALFM1 (N2P1 >=I ALFM1 (N2) 

I ALFM2 ( N2P1 ) = l ALFM2 ( N2 ) 
IALFM3(N2P1)=I ALFM3TN2) 

I ALFM4 { N2P 1 ) = I ALFM4 ( N2 ) 

2176 CONTINUE 
CAPMIN (NN )=CAP 

I ALFM1 { NN ) = I A1 AR 
I ALFM2 { NN ) = I A2 AR 
I ALFM3 ( NN ) = I A3 AR 

1 AL.FM4 ( NN ) = I A4 AR 
2178 CONTINUE 

2180 CONTINUE 
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2185 CONTINUE 
2190 CONTINUE 
2195 CONTINUE 


MAIN PROGRAM CASE 8 

ANGLES GENERATION 

DO 2295 I 1 = 1 ,8 
Ell=FLOAT( II — 1 ) 

A1 = PALF**DALF#E11 
I Al AR=IPALF-IOALF*< I l — l > 
1 1 PL S 1 = 11 + 1 
DO 2290 I2=I1PLS1»9 
E21—FL0AT ( 12—1 ) 
A2=PALF-DALF*E21 
I A2AR=IPALF-I0ALF*{ 12-1 ) 
I 2PLS 1=12+1 
DO 2285 13=11 i 19 
. E3i=FL0AT( 13-1 ) 
A3=PALF-DALF*E31 
I A3AR=IPALF-IDALF«( 13-1 ) 
I 3MIN1= 13—1 

DO 2280 I 4= 1 2PLS1 * 1 3MI N 1 
E4.1=FL0AT( 14-1 > 
A4-PALF-DALF*F41 
I A4AR = I PALF— I DALF’l' ( 14-1 ) 

AIRCRAFT 1 AND 2 

VS=VKTM J/360 0* 

VF=VKT( 2) /36C0* 

AMI =VS/VF 

gamaf=gama ( 2 y 
gamasAgav.a ( 1 ) 
beta=gamaf-gamas 

DELTSS=DELTA( 1 fl ) 

OELTFF= DELTA < 2 ,2 ) 
DELTSF=DELTA( 1 ,2) 

OELTF S— DELTA! 2 » 1 ) 

AS=Al 

AF=A2 

I AFAR= IA2AR ' 

I ASAR=I A1AR 
CQAS=COS( AS) 

C0AF=C0S( AF) 

C0ASAF=C0S ( AS-AF) 
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A.MAS=AMl**2-2* *AMI*C0AS+1 
AMAF= AM I #*2-2 • * AM I*C0AF+1 . 

AM AS AF= AM I **2 — 2 • * AM I *C0 AS AF+l . 

CALL SDOSF 

CALL SDCFS 

CALL SOOFF 

CALL SDOSS 

CALL DOT 

TT < i , 1 ) = TSS 

TT<l*2)=TSF 

TT (2,1) =TFS 

TT { 2 # 2 ) =TFF 

A I RCft AFT 1 AND 3 

VS=VKT< 1 )/3600. 

VF~VKT{ 31/3600 . 

AM I — VS/VF 
GAMAF=GAMA(3) 

GAMAS=GAMA{ 1 ) 

• 3ET A=GAMAF — G AM AS 

deltss=delta{ i,n 

DELTFF=D£LTA{ 3, 3) 

DELTSF=DELTA< 1 ,3) 
DELTFS=DELTA{3,1 ) 

AS=Al 
A Pis A 3 

I AFAR= I A3 AR 
I AS AR = I A 1 AR 
CUAS=COS( AS) 

CCAF=COS( AF) 

COASAF=COS{ AS-AF) 

AM AS= A Ml **2 — 2 . #AM I *COA S + 1 . 

A M AP = A M I # # 2—2 • * AM I * CO AF + 1 • 

AM AS AF— AM l#*2-2«*AMI*COASAF+l . 

CALL SDOSF 

CALL SDOFS 

CALL SDOFF 

CALL SDOSS 

CALL DOT 

TT< 1 *3)=T5F 

TT{ 3 , 1 )=TFS 

TT<3,3)=TFF 

AIRCRAFT l AND 4 

VS=VKT( 1 ) /360C. 

VF=VKT( 4)/ 3600 • 


AMI = VS/VF 

■■■: .• GAMAF=GAM A( A ) ; 

GAMAS=GAMA{ 1 ) 

BETAS GAM AF-G AM AS 
DELTSSsDELT A( 1,1) 
DELTFF=0ELTA(4,4> 

DELTSFsDELTAC 1 ,4) 
DELTFS=0ELTA(4,1 ) 

A S- A 1 
AF=A4 

I AFAR= I A4AR 
I AS AR= I A 1 AR 
COAS=CQS{ AS) 

COAF=COS<AF) 

COAS AF=COS ( AS— AF) 
AMASsAMI**2-2«*AM!*C0AS+l • 
AMAF=AMI**2-2 .*AMI 4C0AF-U . 
AMASAFs AM I#<=2-2 • # AM I *COASAF+l . 
CALL SDOSF 
CALL SDOFS 
- CALL SDOFF 
CALL SOOSS 

CALL OOT : ! ' ; 

TT< I » 4 ) — T SF 
TT<4,1)=TFS 
TT(4,4)=TFF 
C 

C AIRCRAFT 2 AND 3 
C 

VS=V.KT(2) /3600. 

VF=VKT( 3)/3600. 

AM I = VS/VF 
GAMAF=GAMA 13 ) 

G AM AS=G AM A { 2 ) 

BET A = GAM AF— GAM AS. 

OELTSS=DELTA{ 2 ,2 ) 

OELTFF-OELTAI 3, 3) 
OELTSF=DELTAt2 ,3) 

(5ELFS = DELTA(3 t 2) 

AS=A2 
A F— A 3 

I ASAR=I A2AR 

I AFARs IA3AR 

COAS=COS{ AS) 

COAF=COS{ AF) 

CQASAF=COS( AS-AF) 
AMAS=AMI**2-2. *AMI*COAS+l . 

A MAP = AM I $£2— 2 • *AMI 4COAF + 
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AMASAF=AMI**2-2.*AMI*COASAF+l . 
CALL SDOSF 
CALL SDQFS 
CALL SDOFF 

CALL SOOSS . 

CALL DOT 
TT(2,3)=TSF 
T T { 3 , 2 ) =TFS 

AIRCRAFT 2 AND 4 

VS=V<T( 2J/3600. 

VF— VKT ( 4 ) /360 0 • 

AMI — VS/VF 
GAMAF=G.AMA<4) 

GAMAS=GAMA(2I 
BET A=GAMAF— GAMAS 
DELTSS^DELTAI 2,2) 

DELTFF=DELTA( 4 f 4 ) 

DEL T SF = DEL T A ( 2 , 4 ) 
DELTFS=DELTA(4,2 > 

AS=A2 

AF=A4 

I AS AR = I A 2 AS 
I AFAR=IA4AR 
COAS=CQS{ AS) 

COAF=COS<AF) 

COASAF=COS ( AS— AF) 
AMAS=AMI**2-2.*AMI*CQA S+l • — 
AMAF=AMI**2-2 .*AMI*COAF+l • 

AMAS AF= AM l **2-2 . * AM I *CO AS AF F 1 . 

CALL SDOSF 

CALL SDGFS 

CALL SDOFF 

CALL SDOSS 

CALL DOT 

TT ( 2 » 4) =T SF 

TT( 4,2)=TFS 

AIRCRAFT 3 AND 4 

VS=VKT( 3) /3600. 

VF-VKT( 4)/3600. 

AM I = VS/VF 
G AM AF— GAMA ( 4 ) 

G A M A S = G A M A ( 3 ) 

BET A = GAMAF— GAM AS 
DELTSS=DELTA( 3,3 ) 
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DELTFF=0ELTA<4,4) 

DELTSF=DELTA( 3,4) 
DElTFS=0ELTA<4,3) 

AS=A3 

AF=A4 

I ASAR=I A3AR 
I AFAR= I A4 AR 
COAS=COS( AS) 

COAF=COS( AF> 

COASAF=CO S{ AS- AF ) 

AM AS=AM 1 4 *2—2 • 4AMI 4CQAS4- 1 • 

AMAF=AM 1*42-2. 4AMI4COAF+1 . 

AMASAF=AM 1442-2. 4AMI4C0ASAF 4-1 • 

CALL SOOSF 

CALL SDOFS 

CALL SDOFF 

CALL SOOSS 

CALL OOT 

TT(3,4)=TSF 

TT<4,3)=TFS 

MEAN INTERARRIVAL TIME 

TBAR= 0 . 

00 2220 I = 1,M 
DO 2210 J=1,M 
TBAR-TBAR4-TT( I , J)4PP( I , J) 

2210 CONTINUE 
2220 CONTINUE 

CAP-3600. /TSAR 
N = NM 

IF ( CAP .L T » CAP AC I N ) 1 SO TO 2270 
DO 2250 N 1 = l * N 

IF ( CAP-CAP AC( N1 J ) 2250,2250,2240 
2240 CONTINUE 
NN=N 1 

GO TO 2255 
2250 CONTINUE 
2255 CONTINUE 

NMINNN=N— NN 

DO 2260 N0P=1,NMINNN 

N2=N-NOP 

N2P1=N2+1 

CAP AC ( N2P 1 )=CAPAC(N2) 

IALFKN2P1 )- I ALF1 { N2 ) 

I ALF2( N2P1 )=I ALF2IN2) 

I ALF3( N2P 1 ) = ! ALr 3 ( N2 > 

IALF4IN2P1 >=IALF4(N2> 
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2260 CONTINUE 

CAP AC C NN ) = CAP 
I ALF i { NN ) = I A 1 AR 
I ALF2{ NN)=I A2AR 
I ALF 3 C NN ) = I A3 AR 
I ALF 4 ( NN ) = I A4 AR 
2270 CONTINUE 
N=NM 

IF ( C AP • GT * C APM IN ( N ) ) GO TO 2278 
DO 2274 N 1= 1 ,N 

IF { CAP— CAPM I N ( N1 ) ) 2272,2274,2274 

2272 NN=N1 

GO TO 2275 

2274 CONTINUE 

2275 CONTINUE 
NMINNN=N-NN 

00 2276 NOP— 1 , NMI NNN 
N2=N— NOP 

N2PI— N2+ 1 

CAPMININ2P1 )=CAPMIN{N2) 

• IALFM1IN2P1 )=IALFM1(N2) 

I ALFM2IN2PI )= I ALFM2IN2 J 
I ALF M3 ( N2P 1 ) = I ALF M3 ( N2 ) 

I ALFM4 { N2P1 ) - I ALFM4CN2 ) 

2276 CONTINUE 

C APM INI NN ) =C AP 
I ALF M 1 ( NN ) = I A 1 AR 

1 ALF M2< NN ) = I A2 AR 
I ALFM3( NN )= I A3 AR 
I ALF M4 ( NN ) = I A4 AR 

2273 CONTINUE 
2230 CONTINUE 
2235 CONTINUE 
2290 CONTINUE 
2295 CONTINUE 


MAIN PROGRAM CASE C 
ANGLES GENERATION 

, DO 2395 11-1,8 
El l=FLOAT< Il-l ) 

A 1 = P ALF— D ALF* El 1 
I A1AR=IPALF-IDALF*{ Il-l > 
I 1PLS1=I 1+1 
DO 2390 12=11,19 
E21=FLDAT ( 12-1 ) 
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A2=PALF-DALF4E21 
I A2AR=I PALF-IDALF4< 12-1 ) 

I 2M INI =12-1 
00 2385 I 3= I 1 PLS1 »9 
E3i=FL0AT< 13-1 ) 
A3-PALF-DALF4E31 
I A3AR-IPALF-I0ALF»( 13-1 ) 
I3PLS1-T3+1 

DO 2380 I4=I3PLS1, X2MIN1 
£41 -FLOAT < 14-1 ) 

A4=PALF— DALF4E41 
r A4AR=IPALF-IDALF4{ 14-1 ) 

AIRCRAFT l AND 2 

V S= VKT { 1 J/3600. 

VF=VKT( 2) /360C. 

_ AM 1= VS/VF 
GAMAF=GAMAI2) 

G AMAS=GAMA ( 1 ) 

• 8ETA=GAMAF— GAM AS 
OELTS3=DELTA{ 1 ,11 

0 EL TFF = D£ L T A ( 2 , 2 ) 

QELTSF=DELTA{ 1 ,2) 
0ELTFS=DELTA(2,1 ) 

AS=A1 

AF=A2 

1 AFAR- I A2 AR 
I ASAR= I A 1 AR 
COAS=CQS( AS) 

COAF=COS<AF) 

C0ASAF=C0S( AS-AF) 
AMAS=AMI442-2. 4AM I 4C0AS+1 • 
AMAF=AM 1442-2. 4AM I4C0AF+1 • 

AMA SAF= AM I 442—2 • 4 AM 1 4C0ASAF +1 • 
CALL SDOSF 

CALL 5D0FS 
CALL SOOFF 

CALL SDOSS 

CALL DOT 
TT ( 1 , 1 >=TSS 
TTI I, 2 ) — TSF 
T T ( 2 » 1 )=TFS 
TT( 2 » 2 ) =T FF 

AIRCRAFT 1 AND 3 

VS=VKT< 1 ) /3600 . 
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VF=VKT( 3) /3600. 

AM I=VS/VF 
G AM A F=6 AM A < 3 ) 

GAMAS=GAMA{ 1 ) 

BET A=GAMAF— GAMAS 
OELTSS=OELTA( i , 1 ) 
DELTFF=0ELTA<3,3) 

DELTSF=DELTA{ l ,3) 
DELTFS=DELTA(3, 1 ) 

A 3= A 1 
AF=A3 

I AF AR = I A3 AR 
I ASAR= I A 1 AR 
COAS— COSt AS) 

CQAF = COSC AF) 

COASAF=COS( AS-AF) 

AMAS = AMI**2-2.#AMl*COAS-i-I • 
AMAF=AMI**2-2.*AMI*C0AFFl ♦ 
AMASAF-AM I**2-2-.*AMI*C0ASAF-t-l • 
CALL SDOSF 
* CALL SOQFS 
CALL SDOFF 
CALL SDOSS 
CALL DOT 
TT< 1 ,3)=TSF 
TT ( 3 « 1 } =TFS 
TT { 3 f 3 ) —TFF 

AIRCRAFT 1 AND 4 

VS= i/KT ( 1 ) /3600 • 

VF=VKT< 4)/3600. 

AMI — VS/V/F 
GAMAF=GAMA(4 ) 

GAMAS=GAMA{ 1 > 

BET A=GAM AF—GA MAS 
DELTSS=DELTA< 1 ,1) 

DELTFF- DELTA { 4 ,4 ) 

DELTSF=DELTA( 1 ,4) 

DELTFS=DELTA< 4 , 1 ) 

AS— A l 
AF=A4 

I AF AR=IA4AR 
I ASAR=I A l AR 
CO AS=COS { AS ) 

COAF=COSt AF) 

COASAF=COS( AS-AF) 

A M A S- AM I * * 2-2 . * AM I # COA S 1 • 
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AMAF=AM I#' #'2^2 «$AMI $COAF+i . 

A M A S AF = A M I ** 2- 2 . * AMI *C 0 A S AF + 1 . 

CALL 500SF I 

CALL SOOFS 

CALL SOOFF 

■CALL SOOSS 

CALL DOT - - 

TT( 1 ,4>=TSF 
TT ( 4 » 1 )— TFS 
TT{ 4 »,4 ) “TFF 

AIRCRAFT 2 AND 3 

V S — VKT ( 2 ) /36C 0 • 

VF= VK.T (3)73600. 

AMI=VS/VF 
GAMAF=GAMA( 3) 

G AMAS= GAM A ( 2 ) 

3tT A=GAMAF— GAMAS 
DEL TSS= DELTA ( 2 » 2 ) 
DELTFF=DELTA(3,3> 

OELTSF= DELTA (2*3) 
DELFS=DELTA(3,2) .1. 

A S= A 2 
AF=A3 

IASAR=IA2AR 

IAFAR=IA3AR 

COAS=COS( AS) 

COAF=COS( AF) 

COASAF=COS( AS-AF) 
AMAS=AMI#*2-2.*AMI*CQASFi • 
AMAF=AM I**2-2 .SAMItfCOAF+l . . 

A.MASAF=AMI**2~2.*AMI*C0ASAF + 1 . 
CALL SQOSF 
CALL SDOFS 
CALL SDOFF 

CALL SDDSS 

CALL DOT 
TT( 2,3)=TSF 
TT( 3,2)=TFS 

AIRCRAFT 2 AND A 

VS=VKT( 2)73600. 

VF=VKT(4) 73600. 

AMI=VS/VF 
G AM AF=GAM A ( A ) 

GAMAS=GAMA{ 2 ) 
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beta=gamaf-gamas 

DELTSS=DELTA{2,2> 
DELTFF=OELTA(4 t 4) ~ 
DELTSF=DELTA(2,4) 

DELTFS=DELTA (4,2) 

AS=A2 

AF=A4 

I ASAR=IA2AR 
I AFAR= I A4AR 
COAS=COS( AS) 

GQAF=COS(AF) 

CQASAF=COS( AS— AF) 

AM AS" AM I ❖ + 2—2 • *AMI 4C0AS-*-l • 
AMAF=AMI442-2 .4 AM I*COAF+l . 

AMAS AF= AM 1442— 2*4AMI 4CQ ASAF +■ t * 
CALL SDOSF 
CALL SDOFS 
CALL SDOFF 
CALL SDOSS 
CALL DOT 
• Tt(2,4) =TSF 

TT(4,2)=TFS 

AIRCRAFT 3 AND 4 

VS=VKT(3>/3600* 

VF=VKT( 4)/3600. 

AMI— VS/VF 
GAMAF=GAMA(4) 

GAMAS=GAMA<3) 

BET A=6AMAF— GAMAS 
DELTSS=DELTA( 3 ,3 ) 

DELTFF— DELT A ( 4 , 4 ) 

OELTSF- DELT A { 3 , 4 ) 

DELTFS=DELTA (4,3) 

AS=A3 

AF=A4 

I ASAR= I A3 AR 
IAFAR=I A4AR 
COAS=COS { AS ) 

COAF=COS(AF) 

COASAF=COS< AS-AF) 

AM A S= A M 1 4 4 2 -2 • * AM 1 4 CO A S +1 • - 
AM AF = AM 1 4 42—2 • 4 AM I4C0AF+1 . 

AM ASAF= AM 1442— 2*4AMI 4C0ASAF-M . 
CALL SDOSF 
CALL SDOFS 
CALL SDOFF 
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CALL SOQSS 
CALL DOT 
TT<3,4)=TSF 
TT( 4, 3)=TFS 

MEAN IN TER ARRIVAL TIME 

T'3AR=0 . 

DO 2320 1=1 

DQ 2310 J=1,M 
TQAR = T3AR + TT ( I ,J)*PP( I , J)„ 

2310 CONTINUE 
2320 CONTINUE 

CAP=3600*/T3AR 

N=NM 

IF ( CAP • LT • CAPAC C N ) ) GO TO 2370 
DO 2350 N1=1,N 

IF ( CAP— CAP AC INI)) 2350,2350,2340 
2340 CONTINUE 
NN=N 1 

• GO TO 2355 
2350 CONTINUE 
2355 CONTINUE 

NMINNN=N-NN 

DO 2360 NOP=l, NMINNN 

N2—N— NOP 

N2P1=N2+1 

CAPAC (N2P1 )=CAPAC{N2> 

IALF1 (N2P1)— I ALF 1 ( N2 ) 

I ALF2I N2P 11 = 1 ALF 2 ( N2 ) 

I ALF3IM2P1 )=IALF3IN2) 

I ALF 4 ( N2P 1 )-= 1 ALF4 C N2) 

2360 CONTINUE 

CAPAC(NN)— CAP 
I ALF 1 ( NN 1 = I A 1 A R 
I ALF 2 ( NN ) = I A2 AR 
I ALF 3 ( NN) * I A3 AR 
I ALF 4 (NN)=I A4AR 
2370 CONTINUE 
N=NM 

IF (CAP.GT.CAPMIN(N) > GO TO 2378 
DO 2374 Nisi, N 

IF ( CAP-C APMIN IN 1 ) ) 2372,2374,2374 
2372 NN=sNl 

GO TO 2375 

2374 CONTINUE 

2375 CONTINUE ' 

NMI NNN=N— NN 
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00 2376 N0P=1,NMINNN 
N2=N— NOP 
N2P1-N2+1 

CAPMININ2P1 )=CAPMININ2) 
I ALFM1CN2P1 )=IALFM1(N2) 
I ALFM2 ( N2P1 ) ~ I ALFM2I N2) 
I ALFM3IN2P1 >=IALFM3tN2) 

1 ALFM4( N2P1 ) = I Al_FM4 < N2 ) 
2376 CONTINUE 

C APMl NI NN J =CAP 
I ALFM1 { NN )= I A1 AR 
I ALFM2INN )=I A2AR 
I ALFM3I NN ) - 1 A3 AR 
TALFM4< NN)=I A4AR 
2378 CONTINUE 
2380 CONTINUE 
2385 CONTINUE 
2390 CONTINUE 
2395 CONTINUE 


MAIN PROGRAM CASE D 

ANGLES GENERATION 

DO 2495 11-1,9 

El l “FLOAT! M — 1 ) 

A1=PALF-DALF*E1 1 

I A1AR=IPALF-IDALF*< 1.1-1 ) 

I 1 PL S 1 = 1 1 +• l 

DO 2490 1 2— 12,19 

E21=FL0AT{ 12-1 ) 

A2=PALF-DALF*E21 

I A2AR=I PALF— IDALF* ( 12—1 ) 

I2MINl!=I2-l 

DO 2485 13=11,1 2M IN 1 

E31 -FLOAT ( 13—1 1 

A3=PALF-DALF*E31 

I A3AR=IPALF-IOALF«{ 3-1) 

I 3M IN 1 = 1 3— 1 

DO 2480 I 4= I 1PLS1 , 13MIN1 
£41 =FLO AT ( 14—1) 
A4=PALF-DALF4E41 
I A4AR = I PALF-I DALF*< 14-1 ) 

AIRCRAFT 1 AND 2 

V S= VKT { 1 ) /3600 • 
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VF-VKT( 2) /3600. 

: AMI ssVS/VF 
GAMAF-=GAMA(2 ) 

GAMAS=GAM A( 1 ) 

8ETA=GAM AF—GAM AS 
DELTSS=DELTA< 1 ,1 ) 

OELTFF^OEL T A < 2 , 2 > 
DELTSF=DELTA< 1 ,2) 
DELTFS=D£LTA( 2, 1 1 
AS=A1 
AF=A2 

I AFAR= IA2AR 
IASAR=IA1AR 
COAS=COS( AS) 

COAF=COS(AF) 

COASAF=COS( AS-AF) 

AM AS= AM I ♦^'2—2 « * AMI #COAS+ 1 * 

AMAF=AMl**2-2.*AMI*CGAF+l . 

AMASAF= AM I**2-2.*AMI *COASAF-H 

CALL SDOSF 

CALL SO OF S 

CALL SDOFF 

CALL SOOSS 

CALL DOT ------ -- 

TT{1, U =TSS 
TTI1,2)=TSF 
TT<2, 1 )*TfS 
TT(2,2)=TFF 

AIRCRAFT l AND 3 

VS-VKTI 1) /3600. 

VF=V KTt 31/3600. 

AM 1= VS/VF 
G AMAF—GAM A C 3 1 
GAMAS=GAMA( 1 ) 

8ET A= GAM AF— GAMAS 
OELTSS=DELTA{ 1 ,1 ) 

OELTFF =DELT A ( 3 f 3 ) 
DELTSF=OELTA< 1 ,3) 

DELTFS=DELTA (3,1) 

A S= A 1 
AF= A3 

I AFAR= I A3 AR 
I ASAR=l A1 AR 
COAS=COS< AS) 

COAF=COS( AF) 

CGASAF=COS( AS-AF) 
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AMAS=AM I **2-2 • * AM I *COAS+ l * 

AM AF— AM 1**2— 2 «#AM I *C0 AF+1 • 
AMASAF= AM t **2-2.*AM I *COASAF+ 1 • 
CALL SDOSF 
CALL SDOFS 
CALL 5 DOFF 
CALL SDOSS 
CALL DOT 
TTi 1 3} — T SF 
TT £ 3 » I ) — TFS 
TTI 3,3)=TFF 

AIRCRAFT 1 AND 4 

V S= VKT £ 1 ) /360 0 • 

VF=VKT( 4 ) / 3 60 0 • 

AMI— VS/VF 
G AM AF=GAM Ai 4 ) 

GAMAS=GAMA{ 1 ) 

BET A = GAMAF— GAM AS 
• DELTSS=OELTA( 1 ,1 ) 

DELTFF=DELTA< 4 ,4> 

DELTSF»DELTA( 1 ,4} 

OELTFS=DELTA< 4, 1 ) 

AS= A 1 
AF = A4 

I AFAR= I A4 AR 
I ASAR=I A1 AR 
COAS=COS(AS) 

COAF=COS<AF) 

C. O A S AF = CO S ( AS-AF) 

AMAS=AMI**2-2 .*AMI*C0AS+1 • 

AMAF- AM 1**2— 2 • *AM I*CQAF+1 • 

AMASAF=AMI**2-2.*AMI*C0ASAF+1 * 

CALL SDOSF 

CALL SDOFS 

CALL SDOFF 

CALL SDOSS 

CALL DOT 

TT<1,4)=TSF 

TT £ 4 } 1 ) — TF S 

TTt4, 4)=TFF 

AIRCRAFT 2 AND 3 

VS=VKT( 2>/3600. 

VF=VK.T{ 3)/36CG. 

AM I = V S/VF 
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GAMAF=GAMA(3) 

G AM AS=GAM A (2 ) 

BET A = GAM AF— GAMA S 
DELTSS-DELTAt 2,2) 
DELTFF=DELTA<3,3> 

DELTSF=DELTA( 2,3) 

O ELF S= DELTA (3,2) 

AS=A 2 
AF=A3 

I ASAR-I A2AR 
I AFAP.= I A3AR — - 

COAS=CQS(AS) 

COAF==CQS< AF) 

COASAF=COS( AS-AF) 
AMAS=AMI**2-2.*AKI*C0ASFl • 
AMAF=AMI« ; ! t 2-2.*AMI^C0AF+i . 
AMASAF=AM I**2-2 .*AMI *COASAF+l * 
GALL SDOSF 
CALL SDOFS 
CALL SDOFF 
• CALL SDOSS 
CALL DOT 
TT(2,3)=TSF 
TT{ 3,2>=TFS 

AIRCRAFT 2 AND 4 

V$=VKT(2) /3600. 

VF=VKT{ 4>/3600. 

AM I = VS/ VF 
GAMAF=GAMA<4) 

GAM AS=GAM A ( 2 ) 

BET A^GAMAF— GAMAS 
DEL TSS= DELTA (2,2) 
DELTFF=DELTA{4,4) 

DELTSF =DELTA (2,4) 

DELTF S = DcLT A < 4 , 2 ) 

AS- A 2 
AF = A4 

I ASAR=I A2AR 
I AFAR = I A4 AR 
COAS=COS( AS) 

CO Ar — C OS * AF ) 

COASAF-CGS{ AS-AF) 

AM A S= A M I F *2—2 . #AMI *COAS+ I • 
AMAF=AMl**2-2. *AMI#COAF+l . 
AMASAF=AMI**2~2* 3 S : AMI*COASAF*l . 
CALL SDOSF 


221 



uuu oou 


CALL SDOFS 
CALL SDOFF 
CALL SDOSS 
CALL DOt 
TT{ 2,4>=TSF 
TT(4,2>=TFS 

AIRCRAFT 3 AND 4 

VS=VKT( 3) /36C0* 

VF=VKT< 41/3600. 

* AMI=VS/VF 
GAMAF=GAMA{ 4) 

G AM AS=G AM A ( 3 ) 

BETA=GAMAF— GAMAS 
DELTSS=DELTA(3 ,3 ) 
DELTFF=DELTA<4,4) 

DELT SF = DELTA (3*4) 

DELTFS=DELTAC 4,3) 

AS=A3 
- AF=A4 

IASAR=IA3AR 
I AFAR=I A4AW -T- 

COAS-COS(AS) 

C OAF— COS ( AF) 

COASAF = COS< AS-AF) 
AMAS=AMI**2-2.*AMI4C0AS+1 • 
AMAF-AMI**2-2.*AM1*C0AF>1 • 
AMASAF=AMI**2-2.*AMl*C0ASAF+l • 
CALL SDOSF 
CALL SDOFS 
CALL SDOFF 
CALL SDOSS 
CALL DOT 
TT f 3,4 } — T SF 
TT ( 4,3) =TF S 

MEAN INTERARRIVAL TIME 

TBAR=0. 

DO 2420 1-1, M 
DO 2410 J= 1 , M 
TBAR=T3AR*TT( I , J)4PPC I , J) 

2410 CONTINUE 
24'20 CONTINUE 

CAP=3600./TBAR 

N=NM 

IF { CAP »LT • CAP AC ( N ) ) GO TO 2470 
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00 2450 N1=1,N 

IF (CAP-CAP AC (NIT) 2450,2450,2440 

continue 

NN=N1 

GO TO 2455 
CONTINUE 

CONTINUE 

NMI N 1 =N— 1 
NMINNN=N-NN 

00 2460 N0P=1,NMINNN 
N2=N-NQP 

CAP AC { N2P 1 >=CAPAC(N2) 

I ALF1 (N2P1 ) = I ALF 1 { N2 ) 

I ALF 2 ( N2P 1) = I ALF2 ( N2 ) 

IALF3(N2P1 )=IALF3(N2) 

I ALF4( N2P1 )— I ALF4< N2) 

CONTINUE 
CAP AC ( NN ) =CAP 

1 ALF 1 { NN i = I A 1 AR 
I ALF2C NN)=I A2AR 
I ALF 3 ( NN ) = I A3 AR 
t ALF4( NN)=I A4AR 
CONTINUE 

N-NM 

IF (CAP.GT.CAPMIN(N) ) GO TO 247S 
DO 2474 Nl=l ,N 

IF ( CAP-CAPMIN(N1 ) ) 2472,2474,2474 
NN=Nl 

GO TO 2475 

CONTINUE 

CONTINUE 

NMI NNN=N— NN - _ 

DO 2476 N0P=1 , NMINNN 
N2=N— NOP 
N2P1-N2+1 

C AP M I N (iN 2 P 1) = C APM I N ( N 2 ) 

I ALFM1 (N2Pl)=I ALFM 1 ( N2 ) 

I ALF M2 ( N2P1 )- I ALFM2(N2) 

I ALF M 3 ( N 2P 1 ) = I ALF M 3 ( N 2 ) 

I ALFM4< N2P 1 ) = I ALFM4( N2 ) 

CONTINUE 
C APM I N ( NN ) =C AP 
I ALF M l ( NN ) = I A 1 AR 
I ALF M2( NN > = I A2AR 
I ALF M3 ( NN )— I A3AR 
I ALFM4{ NN)=I A4AR 
CONTINUE 
CONTINUE 
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2485 CONTINUE 
2490 CONTINUE 
2495 CONTINUE 


PRINTING RESULTS 


PRINT 4050 

4050 FORMATdHl ,///,13X, ^MAXIMAL C APACI TY-QPT I MAL ANGLES*,////) 

PRINT 4100 

4100 FORMAT ( 1 4 X, * CAPACITY* , 5X , *ALFA1 *, 5X f * ALFA 2* ,5X,*ALFA3*, 

1 5X , * ALFA4* , /// ) 

DO 4150 NP=1,NM 

PRI NT 4151 , NP , CAP AC ( NP ) , I ALF 1 {NP) , IALF2C NP) , IALF3INP) , I ALF 4 ( NP ) 
4151 F ORM AT <5X» 15, FI 0*3,41 10) 

4150 CONTINUE 

PRINT 4250 

4250 FORMAT ( 1H1 ,///, 13X,*MINI MAL CAP AC I TY-THE WORSE ANGLES*,////) 
PRINT 4100 
DO 4260 NP= 1 , NM 

• PRINT 4151 ,NP,CAPMIN<NP) , I ALFM 1 I NP ) , I ALF M2INP) , I ALF M3 (NP), 

1 1 ALF M4 { NP ) 

4260 CONTINUE 

GO TO 9999 
STOP 

END V ' ~ — — .. , 
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DOSS * INDXSS 


90 

sc 

70 

60 

50 

40 

30 

20 

10 

0 

24253 

3.91618 

3.66230 

3.46408 

3.31012 

3.19253 

3. 10582 

3.04628 

3.01146 

3 *0 00 CO 

1000 

1000 

1000 

10C0 

1000 

1000 

1000 

1000 

1000 

1000 


NS 



DOFF 



V/ 90 

80 

70 

60 

50 

40 

30 

20 

10 

0 

90 

4.06898 

-1 

-I 

— I 

-! 

-I 

-I 

-I 

-I 

— I 

80 

3.91 618 

3.78139 

-I 

- 1 

-I 

-I 

‘ - 1 

-I 

-I 

-I 

70 

3.66230 

3.6623C 

3.55705 

- I 

-1 

— I 

-I 


-1 

-1 

60 

3.46403 

3.46408 

3.46408 

3.38213 

-I 

-I 

-I 

-i 

-t 

-I 

sc ; , 

3.31012 

3 . 31 C 1 2 

3.31012 

3.31012 

3.24717 

-I 

-l 

-i 

-! 

-l 

40 

3. 1 9253 

3. 1 9253 

3. 1 9253 

3.19253 

3.19253 

3.14559 

-r 

-i 

- I 

-I 

3C 

3.10582 

3. 1 0532 

3.10532 

3. 10582 

3. 10582 

3.10582 

3.07284 

-i 

-I 

- 1 

20 

3.04623 

3.04628 

3.04628 

3 . 0462 8 

3.04628 

3.04628 

3.C4623 

3.02589 

-I 

-1 

10 

3.01146 

3.01146 

3.01146 

3.0 1 146 

3.01146 

3.01146 

3.01146 

3.01146 

3.0C236 

-I 

0 

3.00000 

3 « COOOO 

3.00000 

3.00000 

3.00000 

3.00000 

3. OOOCC 

3.00000 

3.00000 

3.00286 

10 

3.01146 

3.01146 

3.01146 

3.01 146 

3.01146 

3.01146 

3.01146 

3 ,01 146 

3.01146 

3.01146 

20 

3.04628 

3.04628 

3.04628 

3.04623 

3. 04623 

3.04628 

3.04628 

3.04628 

3.04628 

3.04623 

3 0 

3. 10582 

3.1 058 2 

3.10582 

3. 1 0582 

3. 10502 

3.10582 

3. IC582 

3.10 582 ' 

3. 10582 

3.10582 

40 

3.19253 

3.19253 

3. 1 9253 

3. 19253 

3. 19 253 

3.19253 

3.19253 

3. 19253 

3. 19253 

3.1 9253 

50 

3.31012 

3.31012 

3.31012 

3.31012 

3.31012 

3.31012 

3.31012 

3.31012 

3.31012 

3,31012 

60 

3.46408 

3.4 64 08 

3.46408 

3.46403 

3.46408 

3.46408 

3.46408 

3.464 08 

3.46403 

3.46403 

70 

3.66230 

3.66230 

3.66230 

3.66230 

3.66230 

3.66230 

3.66230 

3.6623C 

3.66230 

3.66230 

80 

3.91618 

3.91613 

3.91618 

3.91613 

3.91618 

3.91618 

3.91618 

3.91618 

3.91618 

3.91613 

90 

4.24258 

4.24258 

4.24258 

4.24258 

4.24258 

4.24253 

4.24258 

4.24259 

4. 24253 

4.24253 
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tNDXFF 



90 
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70 

60 

so 
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0 
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R 
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R 
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R 

n 

SO 

1000 

1000 

R 

R 

R 

R 

• R 

R 

R 

R 

70 

1000 

1C0C 

100 0 

R 

R 

R 

R 
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60 

1000 

1000 

1000 

1000 

R 

R 

R 

R 
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R 

50 

o 
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o 

H 
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R 

R 

40 
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R 

R 

R 

R 
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100 0 

ICOO 

l ooo 

1000 

100 0 
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R 

R 

R 

20 
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100 0 

100 0 

1 000 

1 ooo 

1 OOO 

1000 

R 


10 

10C0 

100 0 

1000 

1000 
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1 ooo 

1 ooo 
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1000 

r-v 

r\ 

0 

1000 

10CC 

100C 

1000 

1000 

1 coc 

icoo 

1000 

1000 

1000 

-10 

1000 

1000 

1000 

1 000 

1 ooo 

l ooo 

1000 
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1000 

ICOO 

-2 0 

10 CO 

1 oco 

1000 

ICOO 
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ICOO 
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100 3 

-30 
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1 000 

1 ooo 

1000 

1 ooo 

100C 
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1000 

-4 0 

1000 

1000 

1000 

1000 

I oco 

ICOO 
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1000 

1000 

1000 

-50 

10 00 

I 000 _ 

1000 

1 00 0 

1 ooo 

l coo 

1030 

1000 

100 0 . 

icco 

1 

o 

icoo 

icoo 

1 C3C 

1000 

1 ooo 

1000 

1000 

1000 

1000 

1000 

-70 

1000 

1000 

1000 

1 000 

1 ooo 

ICOO 

o 

o 

o 

13 00 

1 ooo 

1 C C 3 

-83 

ICOO 

ICOO 

1000 

1000 

1 ooo 

1000 

1000 

1000 

1000 

1000 

-90 

100 3 

1000 

1000 

1000 

1000 

1000 

1000 

1C0C 

1003 

100 0 
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C.3. Example of CAP 4 output (MLS case) 


INPUT DATA 


VKTl 

= 

100.00 

gamai = 

2.00 

PI 

= 

.20 

VKT2 

- 

120.00 

GAMA2 = 

2.00 

P2 

V 

.20 

VKT3 

= 

140.00 

GAM A3 = 

4.00 

P3 

= 

.40 

VKT4 

— 

150.00 

GAM A4 = 

4.00 

P4 

5S 

.20 


DELTA MATRIX 

TRAILING A/C 



1 

2 

3 

4 

I' 

3.00 

3.00 

3.00 

3.0 0 

2' .7.x.' 

3.00 

3. CO 

3.00 

3.00 

3 /: 

3.00 

3.00 

3.00 

3.00 


4 


3.00 


3.0 0 


3.00 


3.00 




MAXIMAL CAPACITY-OPTIMAL ANGLES 



CAPACITY 

ALFA! 

ALFA2 

ALFA3 

ALFA 4 

1 

39. 887 

40 

-30 

10 

-10 

2 

3 9. 887 

40 

-30 

-10 

1 0 

3 

39. 824 

40 

-30 

20 

-10 

A 

39.822 

40 

-30 

-20 

c 

5 

39.782 

50 

-30 

-10 

1C 

5 

39.775 

,40 

-3 0 

-10 

0 

7 

39.750 

50 

30 

-10 

10 

6 

39.743 

50 

-30 

— 10 

30 

9 

39.724 

5C 

-30 

20 

0 

to 

39.719 

50 

-30 

20 

-10 

l 1 

39. 7*.S 

5 0 

-30 

-2 0 

0 

12 

39.7C7 

4 0 

— 30 

20 

10 

13 

39. 705 

40 

-40 

10 

-1C 

1 A 

39.6G4 

4 0 

-20 

10 

0 

15 

39.701 

4 0 

-30 

-20 

-10 

16 

39.591 

40 

-40 

-10 

20 

17 

39.6 86 

5C 

30 

20 

0 

18 

39.682 

50 

30 

20 

-1C 

19 

39.671 

50 

-30 

10 

0 

20 

39.671 

5 C 

-30 

-10 

0 

21 

39.666 

4C 

-40 

10 

-30 

22 

39.656 

4 0 

-20 

2 0 

0 

23 

39.649 

40 

-3 0 

10 

— 2C 

24 

39. 647 

40 

-40 

2 0 

n 

25 

39.64 4 

50 

30 

10 

-40 

2 6 

39.643 

40 

-40 

20 

- 10 

27 

39.642 

50 

30 

20 

-30 

2ft 

39.639 

50 

30 

1 0 

0 

29 

39.639 

50 

30 

-10 

0 

33 

39. 629 

40 

-40 

20 

-20 

31 

39.607 

30 

-30 

-20 

0 

32 

3 9. 6 04 

40 

-40 

20 

-30 

33 

39.603 

50 

-30 

20 

1 0 

34 

3 9.631 

5 C 

-40 

10 

-10 

3 5 

39.593 

5C 

-30 

-20 

-10 

36 

39.595 

40 

-40 

-10 

0 

37 

39.538 

50 

-30 

30 

10 

38 

39.538 

50 

-30 

30 

-10 

39 

3 9. 5 67 

50 

-40 

10 

— 2C 

40 

3 9.5 f-7 

50 

-40 

-10 

20 
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MINIMAL CAPACITY-THE WORST. ANGLES 




CAPACITY 

ALFA l 

ALFA 2 

ALFA3 

ALFA4 


i 

32.284 

90 

80 

-90 

-80 


2 

32.5 74 

90 

-90 

80 

-80 


3 

32.574 

90 

-90 

-80 

80 

Q & 

4 

32.583 

90 

80 

-90 

70 


5 

32.588 

90 

80 

-90 

-7 0 

*xj cFj 

6 

32.719 

90 

10 

-9.0 

-80 

■ o a 

7 

32.741 

90 

70 

-90 

-80 

o S 

8 

32. 81 G 

90 

30 

-90 

60 

& >■ 

g 

32.310 

90 

80 

-90 

-60 

& 

1 0 

22. 833 

90 

-90 

8 0 

70 

40 w 

1 1 

32.383 

90 

-90 

SO 

-70 

a & 

I 2 

32.833 

9C 

-90 

—so 

7 C 

}>■ 

I 3 

32.333 

9 C 

-90 

— QC 

-70 

c* 

I* 

32.944 

90 

80 

-90 

50 

taf ra 

15 

32. 944 

90 

80 

-90 

-50 


16 

33.031 

90 

10 

-90 

-70 

1 7 

33.046 

90 

80 

-90 

4Q 


1 3 

33.046 

90 

80 

-90 

-4Q 


19 

33.053 

90 

70 

-90 

-70 


20 

33. C 72 

90 

60 

-90 

-SO 


21 

33 .0 32 

90 

30 

-90 

30 


22 

33.092 

90 

dO 

-90 

-30 


23 

33.109 

90 

— 30 

SO 

60 


24 

33.109 

90 

-90 

so 

— 60 


2 5 

33. 109 

90 

-90 

-80 

60 


26 

33.1 09 

90 

-90 

-80 

-5 0 


27 

33.109 

90 

3 0 

-9 0 

20 


25 

33. 1 09 

90 

BO 

-90 

-20 


29 

33.120 

9C 

80 

—90 

1 0 


30 

33.120 

90 

SO 

-90 

-10 


31 

3 8. 1 23 

90 

80 

-9 0 

0 


32 

33.156 

90 

-90 

70 

-3 0 


33 

33. 1 56 

90 

-90 

-70 

80 


34 

33.1 85 

80 

70 

-90 

-SC 


3 5 

33. 197 

90 

i 0 

-30 

C 


36 

33. 246 

90 

-90 

80 

50 


37 

33.246 

90 

-90 

-80 

50 


33 

33.246 

90 

—90 

60 

. -50 


3 9 

33.246 

90 

-90 

-00 

-50 


40 

33.259 

90 

1 0 

-90 

-6 0 


N3 

-P- 

VC 


